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Quing  Zhu  (PI) 

INTRODUCTION: 

Diffuse  optical  tomography  in  the  near  infrared  (NIR)  is  an  emerging  modality  with  potential 
applications  in  radiology  and  oncology  [1-8].  Optical  tomography  with  NIR  light  is  made  possible 
in  a  spectrum  window  that  exists  within  tissues  in  the  700-900  nm  NIR  region,  in  which  photon 
transport  is  dominated  by  scattering  rather  than  absorption.  Optical  tomography  offers  insight  into 
functional  parameters,  such  as  tumor  angiogenesis,  required  for  tumor  growth  and  metastasis,  tumor 
hypoxia,  indicator  of  tumor  response  to  various  forms  of  therapy,  and  tumor  metabolism.  If  a  single 
optical  wavelength  is  used,  optical  absorption  related  to  tumor  angiogenesis  and  other  normal  blood 
vessels  can  be  measured.  If  two  or  more  optical  wavelengths  are  used,  both  oxy-hemoglobin  and 
deoxy-hemoglobin  concentrations  can  be  measured  simultaneously.  However,  optical  tomography 
alone  has  not  been  widely  used  in  clinical  studies.  The  fundamental  problem  remains  the  intense 
light  scattering  in  tissue,  which  makes  the  lesion  localization  difficult. 

We  have  developed  a  novel  near-infrared  optical  technique  with  ultrasound  localization  [9-16]. 
During  the  last  year,  we  have  made  significant  progress  in  clinical  studies.  Preliminary  results  with 
a  group  of  biopsied  patients  have  shown  that  early  stage  invasive  cancers  present  two-fold  greater 
total  hemoglobin  concentration  than  fibroadenomas  and  benign  lesions  [Appendices  A  and  B]. 
Preliminary  results  of  advanced  cancers  have  shown  good  correlation  between  cancer  angiogenesis 
distributions  imaged  by  the  hybrid  imaging  technique  and  histological  microvessel  density  counts 
[Appendices  C  and  D].  Our  unique  approach  uses  ultrasound  localization  to  overcome  the  poor 
localization  problem  that  the  diffused  light  has  and  thus  significantly  improves  the  sensitivity  and 
specificity  of  optical  tomography. 

The  original  tasks  proposed  in  this  proposal  were 

Statement  of  Work 
Taskl: 

To  simultaneously  deploy  an  ultrasound  array  and  NIR  sensors  on  a  hand-held  probe  and  to 
perform  joint  acoustic  and  optical  imaging  of  the  target  (months  1-8). 

a.  Complete  NIR  imaging  system  by  adding  two  more  detector  channels  (months  1-3). 

b.  Simultaneously  deploy  an  ultrasound  commercial  probe  and  NIR  sensors  on  a  hand-held 
probe  and  perform  calibration  of  the  combined  probe  with  phantoms  (months  4-6). 

c.  Implementing  Rytov  and  Bom  optical  imaging  reconstmction  algorithms  based  on  the  NIR 

sensor  distribution  shown  in  Figure  7  (months  1-6). 

d.  Perform  joint  acoustic  and  optical  imaging  of  testing  targets  (months  4-8). 

Task  2: 

To  acquire  ultrasound  and  optical  image  data  from  90  patients.  To  correlate  optical  deoxy,  blood 
volume,  absorption,  scattering  functional  parameters  with  ultrasound/mammography  findings 
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and  to  assess  the  value  of  optical  imaging  as  an  adjunct  tool  to  improve  breast  cancer  diagnosis 
and  to  reduce  unnecessary  biopsies. 

a.  Refine  optical  imaging  reconstruction  algorithms  using  clinical  data  (9-24). 

b.  Perform  clinical  studies  (months  9-24). 

Extract  ultrasound  properties  of  the  lesions  and  correlate  the  optical  parameters  with  the  ultrasound 
findings.  Biopsy  results  will  be  used  as  diagnosis  standard  (months  9-24). 

PREGRESS: 

Task  1. 

The  proposed  subtasks  were  all  completed.  The  NIR  imager  shown  in  Fig.  1  was  built,  calibrated 
(10,  12,13),  and  used  at  University  of  Connecticut  Health  Center  and  Hartford  Hospital  to  collect 
data  from  patients  who  were  scheduled  for  breast  biopsy. 


(a) 
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(b) 

Fig.  1.  (a)  Picture  of  the  combined  probe  and  the  frequency  domain  NIR  imager,  (b)  Sensor 
distribution  of  the  combined  probe.  Smaller  circles  in  (b)  are  optical  source  fibers  and  big  circles  are 
detector  fibers.  A  commercial  ultrasound  probe  is  located  at  the  center  of  the  combined  probe  and 
the  optical  source  and  detector  fibers  are  distributed  at  the  periphery  of  the  ultrasound  probe. 


The  optical  imaging  reconstruction  algorithms  have  been  optimized  for  imaging  breast  lesions 
(Appendices  A  and  C).  Briefly,  the  NIR  reconstruction  takes  advantages  of  ultrasound  localization 
of  lesions  and  segments  the  imaging  volume  into  finer  grid  in  lesion  region  and  coarser  grid  in  non¬ 
lesion  regions.  A  modified  Bom  approximation  is  used  to  relate  the  scattered  field  measured  at  the 
source-detector  pairs  to  absorption  variations  in  each  volume  element  of  two  regions  within  the 
sample.  With  this  dual-mesh  scheme,  the  inverse  optical  reconstmction  is  well-conditioned  and 
converges  in  few  iterations. 

Task  2: 

To  acquire  ultrasound  and  optical  image  data  from  90  patients.  To  correlate  optical  deoxy,  blood 
volume,  absorption,  scattering  functional  parameters  with  ultrasound/mammography  findings 
and  to  assess  the  value  of  optical  imaging  as  an  adjunct  tool  to  improve  breast  cancer  diagnosis 
and  to  reduce  unnecessary  biopsies. 


a.  Refine  optical  imaging  reconstmction  algorithms  using  clinical  data  (9-24). 
c.  Perform  clinical  studies  (months  9-24). 

Extract  ultrasound  properties  of  the  lesions  and  correlate  the  optical  parameters  with  the  ultrasound 
findings.  Biopsy  results  will  be  used  as  diagnosis  standard  (months  9-24). 
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The  proposed  subtasks  are  nearly  completed.  We  had  initial  delay  on  clinical  study  at  Hartford 
Hospital.  With  the  help  of  Hartford  Hospital  Co-PI  Dr.  Edward  Cronin  and  Integrated  Breast  Care 
Program  Coordinator,  Roxanne  Rotondaro,  we  have  successfully  recruited  70  biopsy  patients  and  are 
continuously  recruiting  more  patients  to  our  study.  We  expect  to  complete  the  enrollment  of  90 
patients  by  the  end  of  the  year.  We  have  another  non-DOD  funded  clinical  site  at  the  University  of 
Connecticut  Health  Center  on  combined  breast  imaging  study  and  57  patients  have  been  recruited  to 
our  study. 

Initial  results  obtained  from  both  clinical  sites  are  very  encouraging  and  have  shown  that  early  stage 
invasive  cancers  present  two-fold  greater  total  hemoglobin  concentration  than  fibroadenomas  and 
other  benign  lesions  [Appendice  A  and  B].  Initial  results  of  advanced  cancers  have  shown  that  the 
total  hemoglobin  distributions  are  quite  heterogeneous  and  the  distributions  correlate  with 
histological  micro  vessel  density  counts  [Appendices  C  and  D]. 

Examples  are  given  here. 

Examples  of  small  invasive  cancers 

The  first  example  is  obtained  from  a  55-year-old  woman.  Figure  2  (a)  shows  a  gray  scale  ultrasound 
image  of  a  non-palpable  lesion.  The  lesion  was  located  at  the  4  o’clock  position  of  the  right  breast. 
The  ultrasound  showed  a  nodular  mass  with  internal  echoes  and  the  lesion  was  considered 
suspicious.  The  estimated  lesion  diameters  measured  from  two  orthogonal  ultrasound  images  was  8 
mm.  An  ultrasound  guided  core  needle  biopsy  was  recommended  and  biopsy  result  revealed  that  the 
lesion  was  intraductal  and  infiltrating  ductal  carcinoma  (nuclear  grade  II,  histological  grade  III).  The 
cancer  once  removed  from  the  breast  measured  1  cm  in  greater  diameter  and  was  composed 
predominantly  of  invasive  carcinoma  (>95%)  extending  to  surgical  margins. 

The  optical  absorption  maps  at  780nm  and  830  nm  are  shown  in  Fig.  2(b)  and  (c),  respectively.  In 
both  (b)  and  (c),  the  first  slice  is  0.7  cm  deep  into  the  breast  tissue  from  the  skin  surface  and  the  last 
slice  is  closer  to  the  chest  wall.  The  spacing  of  the  slices  is  0.5  cm.  The  horizontal  and  vertical  axes 
of  each  slice  are  spatial  x  and  y  dimensions  of  9  cm  in  size.  The  lesion  is  well  resolved  in  slice  #5 
and  has  shown  much  larger  spatial  extension  at  830  nm  than  that  at  780  nm.  The  measured 
maximum  absorption  coefficients  are  pf°=0.24  cm"'  and  pf“=0.28  cm"',  respectively,  and  the 

absorption  maximums  at  both  wavelengths  are  located  at  (0  1 .0  2.25)  cm,  which  is  close  to  the  lesion 
center  measured  by  ultrasound. 

The  total  hemoglobin  distribution  of  the  lesion  is  shown  in  Fig.2(d).  The  measured  maximum  total 
hemoglobin  concentration  for  lesion  is  122.03  p  moles,  and  the  measured  average  background 
hemoglobin  concentration  is  13.73  p moles.  The  spatial  extent  of  lesion’s  angiogenesis  is  measured 
from  the  full  width  at  half  maximum  (FWHM)  of  the  total  hemoglobin  map  and  it  is  estimated  to  be 
2  cm.  This  number  is  about  two  times  larger  than  the  8  mm  diameter  measured  by  ultrasound  and 
suggests  that  optical  contrasts  extend  well  beyond  the  cancer  periphery  because  of  angiogenesis. 
This  example  suggests  that  the  1cm  sample  removed  from  the  breast  with  invasive  cells  extending  to 
surgical  margins  is  not  large  enough  and  2  cm  cancer  size  imaged  with  NIR  may  be  a  good 
estimation. 
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Another  example  was  obtained  from  a  56-year-old  woman  who  had  a  non-palpable  lesion  located  at 
the  10  o’clock  position  of  the  left  breast.  Ultrasound  showed  a  solid  mass  with  internal  echoes 
measuring  9  mm  in  size  and  the  lesion  was  considered  suspicious  (see  Fig.3(a)).  An  ultrasound 
guided  core  needle  biopsy  was  recommended  and  biopsy  results  revealed  that  the  lesion  was  in-situ 
and  invasive  ductal  carcinoma  with  ductal  and  lobular  features  (nuclear  grade  II,  histological  grade 
II).  The  tumor  once  removed  from  the  breast  measured  1.5  cm  in  greater  diameter  and  was 
composed  predominantly  of  invasive  carcinoma  (>80%)  extending  to  inferior/anterior  surgical 
margin. 

The  optical  absorption  maps  at  both  780  nm  and  830  nm  wavelengths  are  shown  in  Fig.3  (b)  and  (c), 
respectively.  In  (b)  and  (c),  the  first  slice  is  0.4  cm  deep  into  the  breast  tissue  from  the  skin  surface 
and  the  last  slice  is  closer  to  the  chest  wall.  The  spacing  between  the  slices  is  0.5  cm.  This  lesion  is 
well  resolved  in  slice  #4  and  has  shown  much  larger  spatial  extension  at  780  nm  than  that  at  830  nm. 
The  measured  maximum  absorption  coefficients  are  pf°=0.29  cm'*  and  =  0.25  cm'*, 

respectively.  The  absorption  maximums  at  both  wavelengths  are  located  at  (0  1.0  1.9000)  cm, 

which  is  very  close  to  the  lesion  center  measured  by  ultrasound.  The  calculated  maximum  total 
hemoglobin  concentration  of  the  lesion  is  127.60  p.  moles,  and  the  background  concentration  is 
24.16  p  moles.  The  measured  FWHM  of  the  total  hemoglobin  map  is  3.5  cm,  which  is  more  than 
three  times  larger  than  9  mm  diameter  measured  by  ultrasound  and  again  suggests  that  optical 
contrasts  extend  well  beyond  the  cancer  periphery  due  to  angiogenesis. 


780nm 


830  nm 


(c)  (d) 

Fig.  2,  (a)  shows  a  gray  scale  ultrasound  image  of  a  non-palpable  lesion  of  a  55-year-old  woman.  The  lesi 
pointed  by  the  arrow  was  located  at  the  4  o’clock  position  of  the  right  breast  measuring  8  mm  in  diamet 
Ultrasound  showed  nodular  mass  with  internal  echoes  and  the  lesion  was  considered  suspicious,  (b)  and  (v 
are  optical  absorption  maps  at  780  nm  and  830  nm,  respectively,  x  and  y  axes  are  spatial  dimensions  of  9  ci 
by  9  cm.  The  slices  beginning  from  left  to  right  correspond  to  absorption  maps  from  0.7  cm  underneath  sk’ 
surface  to  chest  wall  in  0.5  cm  spacing,  (d)  is  the  total  hemoglobin  concentration  and  the  vertical  scale  it 
micro  moles. 
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Fig.  3.  This  example  was  obtained  from  a  56-year-old  woman  who  had  non-palpable  lesion  located  at  the 
o’clock  position  of  the  left  breast.  Ultrasound  showed  a  solid  mass  (a)  with  internal  echoes  measuring  9  n 
in  size  and  the  lesion  was  considered  suspicious,  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  83 
nm,  respectively.  The  slices  beginning  from  left  to  right  correspond  to  absorption  maps  from  0.4  ci 
underneath  skin  surface  to  chest  wall  in  0.5  cm  spacing,  (d)  is  the  total  hemoglobin  concentration  and  tf 
vertical  scale  is  in  micro  moles. 


Example  of  fibroadenoma  and  another  benign  lesions 

Figure  4  (a)  shows  an  ultrasound  image  of  a  hypoechoic  mass  of  a  37-year-old  woman.  TL. 
diagnosis  was  that  the  lesion  likely  was  a  fibroadenoma;  however,  there  was  concern  that  the  lesion 


could  be  a  carcinoma  because  of  the  irregular  shape  and  the  posterior  shadowing  seen  by  ultrasound 
An  ultrasound  guided  core  breast  biopsy  revealed  that  the  lesion  was  simply  a  fibroadenoma. 

Optical  absorption  maps  are  shown  in  Fig.  4(b)  and  (c)  as  well  as  the  total  hemoglobin  distribution  ii 
(d).  Compared  with  the  invasive  cancer  cases,  the  spatial  distributions  of  the  absorption  coefficient  ■ 
as  well  as  the  total  hemoglobin  concentration  are  quite  diffused  and  higher  absorption  region  at  S'" 
nm  corresponds  to  lesion  region  seen  by  ultrasound.  The  measured  maximum  lesion  absorpti  i 
coefficients  at  780  nm  and  830  nm  are  |if'’  =  0.10  cm''  and  =  0.12  cm'',  respectively,  and  the 

values  are  less  than  half  of  those  small  invasive  cancer  cases.  The  calculated  maximum  tot.  ^ 
hemoglobin  concentration  is  52  ji  moles  and  the  background  hemoglobin  concentration  is  U 
fi  moles.  The  calculated  average  maximum  total  hemoglobin  concentration  from  four  NIR  images 
acquired  at  lesion  region  is  59  moles  (±5.1  moles)  and  the  average  background  is  16  moles 
(±1.0  (I  moles). 
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Fig.4.  (a)  shows  an  ultrasound  image  of  a  solid  hypoechogenic  lesion  located  at  9  to  10  o’clo(  . 
position  of  a  37-year-old  woman.  The  diagnosis  was  that  the  lesion  could  represent  fibroadenonr.  , 
however,  there  was  a  chance  that  the  lesion  could  be  a  carcinoma  because  of  the  in'egular  shape  an  I 
posterior  shadow  seen  by  ultrasound,  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  830  nir , 
respectively.  The  slices  starting  from  0.5  cm  underneath  skin  surface  to  chest  wall  in  0.5  cr.i 
spacing,  (d)  is  the  total  hemoglobin  concentration. 


The  statistics  of  maximum  total  hemoglobin  concentration  of  three  groups  obtained  from  the  1' 
cases  are  given  in  Ref  15  (Appendix  A).  The  benign  group  of  fibroadenoma  (15  cases)  an  ’ 
hyperplasia  (1  case)  has  an  average  of  67  ji  moles  (±  17.0  p  moles),  the  combined  fibrocystic  chan 
and  non-invasive  neoplasia/carcinoma  in  situ  case  has  a  maximum  of  48  p  mole,  and  the  invas 
cancer  group  of  2  cases  have  shown  about  twofold  greater  average  of  1 19  p  moles  (±1.6  pmole; 

If  average  total  hemoglobin  concentration  is  use,  the  values  are  46  p  moles  (±  11.3  p  moles),  3 
p  moles  and  86  p  moles  (+2.12  p  moles),  for  the  three  corresponding  groups,  respectively.  The 
malignant  group  also  presents  about  two-fold  greater  average  hemoglobin  concentration  than  that  o*' 
the  benign  group.  The  average  sizes  of  lesions  of  the  three  groups  measured  by  ultrasound  are  1.05 
cm  (+0.3  cm),  1.1  cm,  0.9  cm  (±  0.07  cm),  respectively.  The  lesion  size  is  the  geometric  mean  of 
diameters  measured  in  x  and  z  dimensions. 

Recently,  we  have  completed  data  analysis  of  15  patients  (5  malignant  cancers  and  16  beni;  + 
lesions).  These  15  patients  were  representative  examples  of  the  70  patients  recruited  fror  i 
Hartford  Hospital.  The  results  of  this  group  of  biopsied  patients  have  shown  that  sma  1 
invasive  cancers  can  be  well-resolved  in  3D  optical  absorption  maps.  The  average  maximui  t 
total  hemoglobin  concentrations  for  malignant  and  benign  groups  are  110.3 /mmoles  and  47 
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^  moles,  respectively.  The  malignant  group  presents  more  than  two-fold  greater  total 
hemoglobin  concentration  than  that  of  the  benign  group.  These  results  agree  with  earlier  data 
obtained  from  University  of  Connecticut  Health  Center.  The  manuscript  that  reports  these 
flndings  is  attached  (Appendix  B). 

Evidence  that  hybrid  imaging  technique  can  be  used  for  assessing  chemotherapy  response  and 
tumor  perfusion 

During  the  diagnostic  imaging  study,  we  scanned  a  patient  who  was  undergoing  chemotherapy.  This 
44-year-old  woman  had  a  large  4  cm  x  4  cm  x  1.5  cm  palpable  mass  (Fig.5(a))  located  at  the  6  to  8 
o’clock  position  of  the  left  breast.  The  lesion  was  considered  as  highly  suspicious  for  malignancy 
and  an  ultrasound  guided  needle  biopsy  revealed  that  the  lesion  was  a  high-grade  invasive  carcinoma 
with  necrosis.  Optical  absorption  maps  of  both  wavelengths  are  shown  in  Fig.  5(b)  and  (c)  and  the 
distributions  are  highly  heterogeneous  with  high  absorption  at  cancer  periphery.  Slice  1  is  the  spatial 
x-y  image  of  9  cm  x  9  cm  obtained  at  0.5  cm  deep  from  the  skin  surface.  Slice  7  is  3.5  cm  deep 
toward  the  chest  wall  and  the  spacing  between  slices  is  0.5  cm.  The  total  hemoglobin  concentration 
map  is  shown  in  Fig.5  (d)  and  the  maximum  and  average  hemoglobin  concentrations  are  92.12 
p  moles  and  26.20  p  moles.  Since  this  cancer  was  too  large  for  breast  conserving  surgery,  the 
patient  was  treated  with  chemotherapy  in  the  neo-adjuvant  setting  for  three  months.  At  the  time  the 
patient  completed  the  chemotherapy,  we  imaged  her  again  with  the  combined  probe.  Fig.6  (a)  is  the 
ultrasound  image  of  the  caner  three  month  later  and  the  cancer  contrast  was  poor  and  cancer 
boundaries  were  completely  unclear  due  to  treatment.  Fig.6  (b)  and  (c)  are  optical  absorption  maps 
at  both  wavelengths  and  (d)  is  the  total  hemoglobin  distribution.  The  maximum  and  average 
hemoglobin  concentrations  of  the  lesion  are  83.03  p  moles  and  26.07  p moles,  respectively. 
Compared  with  the  images  acquired  before  treatment,  the  spatial  extent  of  the  angeogenesis  pattern 
is  much  smaller  and  more  confined  to  the  core  area.  The  maximum  total  hemoglobin  concentration 
is  reduced  by  about  10  p  moles  and  the  average  is  about  the  same  as  before.  This  example  clearly 
demonstrates  the  feasibility  of  monitoring  the  treatment  using  the  combined  technique. 

To  correlate  the  NIR  images  with  vessel  densities.  Dr.  Poomima  Hegde,  a  pathologist  at  the 
University  of  Connecticut  Health  Center,  performed  micro-vessel  density  counts.  Two  samples 
obtained  at  breast-conservation  surgery  marked  with  lateral  and  anterior  positions  were  used  for 
counting.  Sections  3  to  5  micrometer  thick  were  stained  on  an  immunohistochemistry  slide  staining 
system  (DAKO  autostainer)  with  factor  8/86  mouse  monoclonal  antibody  (anti-human  von 
Willebrand  factor,  DAKO  Corp,  Carpinteria,  Calif)  at  1:100  dilution  digested  by  proteinase  K  for  3 
minutes,  by  labeled  polymers  (DAKO  EnVision  plus)  immimoperoxidase  method.  The  microvessel 
density  counts  were  performed  in  ten  consecutive  fields  with  the  use  of  an  ocular  grid  at  X200 
magnification.  The  first  field  chosen  was  a  hotspot  (area  of  maximum  vascular  density  either  within 
the  infiltrating  tumor  mass  or  at  the  tumor-stromal  interface).  The  total  number  of  microvessels 
were  196  (lateral)  and  114  (anterior)  per  10  consecutive  fields,  respectively.  These  high  counts 
correlate  well  with  the  high  optical  absorption  and  total  hemoglobin  concentration  shown  in 
Fig.6(b)-(d)  NIR  images. 


Figure  5.  (a)  Ultrasound  image  of  a  large  cancer  obtained  from  a  44-year-old  woman  before  chemotherapj 
(b)  Optical  absorption  map  at  780  nm.  (c)  Optical  absorption  map  at  830  nm.  (d)  Total  hemoglob-' 
concentration  of  the  lesion.  Slice  #1  is  0.5  deep  from  the  skin  surface.  Slice  #7  is  3.5  cm  deep  from  the  sk 
surface.  The  spacing  between  slices  is  0.5  cm.  [from  Ref.  14]. 
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(c)  (d) 

Figure  6(ID19P10).  (a)  Ultrasound  image  of  the  cancer  treated  with  chemotherapy  for  three  month; 
(b)  Optical  absorption  map  at  780  nm.  (c)  Optical  absorption  map  at  830  nm.  (d)  Total  hemoglob>' 
concentration  of  the  lesion. 

The  second  imaging  example  was  obtained  from  a  47-year-old  woman  who  had  a  lump  at  2  o’cloc 
in  her  left  breast.  Ultrasound  showed  hypoechoic  mass  with  irregular  margins  (see  Fig.  7(a))  and  th  i 
lesion  was  considered  as  highly  suspicious  for  malignance.  Fig.7  (b)  and  (c)  are  optical  absorpti 
maps  obtained  at  780  nm  and  830  nm  and  the  vertical  scale  of  (b)  and  (c)  is  half  of  that  used  in  Fi'  i 
(b)  and  (c).  Fig.  7(d)  is  the  total  hemoglobin  concentration  distribution  and  the  vertical  scale  is  , 
same  as  that  used  in  Fig.  6(d).  The  light  absorption  at  both  wavelengths  is  much  lower  than  f 
previous  example  but  the  distributions  are  highly  heterogeneous.  The  measured  maximu 
absorption  coefficients  at  780  nm  and  830  nm  are  0.08  cm'‘  and  0.10  cm  ',  respectively.  Tl]e 
measured  maximum  total  hemoglobin  concentration  of  the  tumor  and  average  background  are  40.'' 
|j,  moles  and  10.6  p.  moles,  respectively.  Surgical  pathology  report  revealed  that  the  mass  was 
infiltrating  carcinoma  (histological  grade  II,  nuclear  grade  II)  with  low  mitotic  activity.  The  total 
counts  of  microvessels  obtained  from  anterior  and  posterior  core  biopsy  samples  were  61  and  40  p'' 
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10  consecutive  fields,  respectively.  The  total  counts  measured  from  anterior  and  posterior  surgical- 
samples  were  52  and  29,  respectively.  These  low  counts  coiTelate  well  with  the  low  optica! 
absorption  shown  in  Fig.  7  (b)-(c)  and  indicate  that  the  tumor  was  poorly  perfused.  : 


These  two  examples  demonstrate  that  NIR  images  correlate  with  microscopic  angiogenesi' 
distributions  and  the  NIR  imaging  technique  has  a  great  potential  for  assessing  tumor  perfusion  ar  . 
therefore  predicting  treatment  responses. _ 
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(c)  (d) 

Figure  7  (IDS),  (a)  Ultrasound  image  of  a  47-year-old  woman  with  a  highly  suspicious  lesion  located  at 
2’oclock  position  of  her  left  breast,  (b)  Optical  absorption  map  at  780  nm.  (c)  Optical  absorption  map  at  830 
nm.  Note  the  scales  for  Fig.  7(b)  and  (c)  is  half  of  that  used  in  Fig.  6(b)  and(c).  (d)  Total  hemoglob* 
concentration  of  the  lesion.  The  scale  is  the  same  as  the  Fig.6  (d). 
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To  establish  the  statistic  power  of  this  novel  technique  on  cancer  diagnosis,  we  are  continuously 
recruiting  more  patients  to  this  study. 

KEY  RESEARCH  ACCOMPLISHMENT: 

•  Initial  results  with  a  group  of  biopsied  patients  have  shown  that  early  stage  invasive  cancers 
present  two-fold  greater  total  hemoglobin  concentration  than  fibroadenomas  and  benign  lesions 
[Appendices  A  and  B]. 

•  Preliminary  results  of  advanced  cancers  have  shown  good  correlation  between  cancer 
angiogenesis  distributions  imaged  by  the  hybrid  imaging  technique  and  histological  microvessel 
density  counts  [Appendices  C  and  D]. 
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with  ultrasound  localization,”  submitted  to  Proc.  of  Nat.  Ac.  of  Sci.  (July  2003). 

Chen  NG,  Huang  MM,  Xia  HJ,  Piao  DQ,  and  Zhu  Q,  “Potable  near  infrared  diffusive  light  imager 
for  breast  cancer  detection,”  submitted  to  Journal  of  Biomedical  Optics. 

Yuan  BH,  Chen  NG,  Zhu  Q.,  “Absorption  and  Emission  Properties  of  Indocyanine  Green  in 
Intralipid  Solution,”  submitted  to  Journal  of  Biomedical  Optics,  July  2003. 

Huang  MM  and  Zhu  Q,  “Dual-mesh  Optical  Imaging  Reconstruction  with  Depth  Correction,” 
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To  be  submitted: 

Zhu,  Q,  Cronin  E,  Currie  A,  Huang  MM,  Chen  NG,  ChenX,  Benign  versus  Malignant  Breast 
Masses:  Optical  Differentiation  with  US  Localization  (work  in  progress),  to  be  submitted  to 
Radiology. 

Conference  Proceedings  (Since  August  2002) 

Zhu  Q,  Kurtzman,  S,  Chen,  NG,  Zarfos,  K,  Huang,  MM,  Kane  M,  “Breast  lesion  diagnosis  using 
combined  near  infrared  and  ultrasound,”  SPIE  Photonics  West,  San  Jose,  CA,  January,  2003.  in 
SPIE  Proceeding  of  Optical  Tomography  and  Spectroscopy  of  Tissue,  in  press. 

Chen  NG  and  Zhu  Q,  “Spread  spectrum  time-resolved  photon  migration  imaging  system:  the 
principle  and  simulation  results”,  SPIE  Photonics  West,  San  Jose,  CA,  January,  2003,  in  SPIE 
Proceeding  of  Optical  Tomography  and  Spectroscopy  of  Tissue,  in  press. 

Chen  NG,  Piao  DQ,  Xia  HJ  and  Zhu  Q,”  Portable  multi-channel,  multi-wavelengh  near  infrared 
imager,”  SPIE  Photonics  West,  San  Jose,  CA,  January,  2003,  in  SPIE  Proceeding  of  Optical 
Tomography  and  Spectroscopy  of  Tissue,  in  press. 

Chen  NG  and  Zhu  Q,  “Fluorescence  decay  profile  measured  by  spread  spectrum  excitation  method,” 
SPIE  Photonics  West,  San  Jose,  CA,  January,  2003,  in  SPIE  Proceeding  of  Optical  Tomography  and 
Spectroscopy  of  Tissue,  in  press. 

Huang  MM,  Chen  NG,  Yuan,  B,  Zhu  Q,”  Three-dimensional  near-infrared  diffusive  wave  imaging 
with  ultrasound  localization,”  SPIE  Photonics  West,  San  Jose,  CA,  January,  2003.  in  SPIE 
Proceeding  of  Optical  Tomography  and  Spectroscopy  of  Tissue,  in  press. 

Guo  PY,  Yan  SK  and  Zhu  Q,  “Three-Dimensional  simultaneous  ultrasound  and  near-infrared 
imaging,”  SPIE  Photonics  West,  San  Jose,  CA,  January,  2003.  in  SPIE  Proceeding  of  Optical 
Tomography  and  Spectroscopy  of  Tissue,  in  press. 

CONCLUSIONS: 

We  have  developed  a  novel  near-infrared  optical  technique  with  ultrasound  localization.  During  the 
last  year,  we  have  made  significant  progress  on  using  our  technique  to  improve  breast  cancer 
diagnosis.  Initial  results  with  a  group  of  biopsied  patients  have  shown  that  early  stage  invasive 
cancers  present  two-fold  greater  total  hemoglobin  concentration  than  fibroadenomas  and  benign 
lesions.  Initial  results  of  advanced  cancers  have  shown  good  correlation  between  cancer  angiogenesis 
distributions  imaged  by  the  hybrid  imaging  technique  and  histological  microvessel  density  counts. 

We  believe  our  unique  has  a  great  potential  to  improve  current  practice  on  breast  cancer  diagnosis  and 
monitor  treatment  responses  of  advanced  cancers. 
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Abstract 

The  diagnosis  of  soiid  benign  and  malignant  tumors 
presents  a  unique  chalienge  to  all  noninvasive  imaging 
modalities.  Ultrasound  is  used  in  conjunction  with 
mammography  to  differentiate  simple  cysts  from  solid 
lesions.  However,  the  overlapping  appearances  of 
benign  and  malignant  lesions  make  ultrasound  less 
useful  in  differentiating  solid  lesions,  resulting  in  a 
large  number  of  benign  biopsies.  Optical  tomography 
using  near-infrared  diffused  light  has  great  potential 
for  imaging  functional  parameters  of  1)  tumor  hemo¬ 
globin  concentration,  2)  oxygen  saturation,  and  3) 
metabolism,  as  well  as  other  tumor  distinguishing 
characteristics.  These  parameters  can  differentiate 
benign  from  malignant  lesions.  However,  optical 
tomography,  when  used  alone,  suffers  from  low  spatial 
resolution  and  target  localization  uncertainty  due  to 
intensive  light  scattering.  Our  aim  is  to  combine 
diffused  light  imaging  with  ultrasound  in  a  hovel  way 
for  the  detection  and  diagnosis  of  solid  lesions,  ihitial 
findings  of  two  early-stage  invasive  carcinomas,  one 
combined  fibroadenoma  and  fibrocystic  change  with 
scattered  foci  of  lobular  neoplasia/lobular  carcinoma 
in  situ,  and  16  benign  lesions  are  reported  in  this 
paper.  The  invasive  cancer  cases  reveal  about  two-fold 
greater  total  hemoglobin  concentration  (mean  119  mol) 
than  benign  cases  (mean  67  mol),  and  suggest  that  the 
discrimination  of  benign  and  malignant  breast  lesions 
might  be  enhanced  by  this  type  of  achievable  optical 
quantification  with  ultrasound  localization.  Further¬ 
more,  the  small  invasive  cancers  are  well  localized  and 
have  wavelength-dependent  appearance  in  optical 
absorption  maps,  whereas  the  benign  lesions  appear 
diffused  and  relatively  wavelength-independent. 
Neoplasia  (2003)  5 
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Introduction 

Ultrasound  imaging  Is  a  well-developed  medical  diagnostic 
tool  that  is  used  in  conjunction  with  mammography  for  the 
differentiation  of  simple  cysts  from  solid  lesions  in  breast 


examinations  [1].  When  the  criteria  for  a  simple  cyst  are  strictly 
adhered,  to,  the  accuracy  of  ultrasound  is  96%  to  100%  [2]. 
Ultrasound  specificity  in  breast  cancer  detection,  however,  is 
not  high  as  a  result  of  the  overlapping  characteristics  of  solid 
benign  and  malignant  lesions  [3,4].  The  insufficient  predictive 
value  of  ultrasound  has  prompted  radiologists  to  recommend 
biopsies  on  rhost  solid  nodules,  which  result  in  a  large  number 
of  biopsies  yielding  benign  breast  lesions.  Currently,  70%  to 
80%  percent  of  such  biopsies  yield  benign  findings  [5].  In 
addition,  the  diagnostic  accuracy  of  ultrasound  depends  largely 
on  the  experience  of  physicians. 

Diffuse  optical  tomography  in  the  near-infrared  (NIR)  is  an 
emerging  modality  with  potential  applications  in  radiology  and 
oncology  [6-13]!  Optical  tomography  with  NIR  light  is  made 
possible  in  a  spectrum  window  that  exists  within  tissues  in  the 
700-  to  900-nm  NIR  region,  in  which  photon  transport  is 
:  dominated  by  scattering  rather  than  absorption.  Optical  tomog- 
5  rapihy  offers  several  functional  parameters  to  differentiate 
malignant  cancers  from  benign  lesions.  It  has  been  shown  that 
breast  cancers  have  higher  blood  volumes  than  nonmalignant 
tissues  because  of  angiogenesis  [1 4].  Tumor  blood  volume  and 
microvascular  density  are  parameters  anatomically  and  func¬ 
tionally  associated  with  tumor  angiogenesis.  If  a  single  optical 
wavelength  is  used,  optical  absorption  related  to  tumor  angio¬ 
genesis  and  other  normal  blood  vessels  can  be  measured.  If 
two  or  more  optical  wavelengths  are  used,  both  oxyhemoglobin 
(oxyHb)  and  deoxyhemoglobin  (deoxyHb)  concentrations  can 
be  measured  simultaneously.  However,  optical  tomography 
alone  has  not  been  widely  used  in  clinical  studies.  The  funda¬ 
mental  problem  of  intense  light  scattering  remains,  which 
makes  the  target  localization  difficult,  especially  for  small 
lesions.  In  general,  localization  or  imaging  based  on  tomo¬ 
graphic  inverse  scattering  approaches  suffers  from  low  spatial 
resolution  and  location  uncertainty,  and  the  inversion  problem 
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is,  in  general,  underdetermined  and  ill  posed.  The  image 
reconstruction  results  depend  on  many  parameters,  such  as 
the  system  signal-to-noise  ratio,  regularization  schemes 
used  in  inversion,  and  so  on. 

Due  to  the  abovedescribed  deficiencies,  literature  data  on 
opticai  tomography  have  been  limited  to  feasibility  studies  or 
case  reports  [7-13].  However,  these  deficiencies  can  be 
overcome  by  a  novel  technique  using  a  combination  of  NIR 
diffusive  light  and  ultrasound  [15-18].  The  combined  imager 
consists  of  a  hand-held  probe  housing  a  commercial  ultra¬ 
sound  probe  and  NIR  source  and  detectors  for  coregistra¬ 
tion.  Ultrasound  is  used  to  scan  suspicious  lesions,  whereas 
coregistered  opticai  measurements  are  used  to  image  and 
characterize  the  functional  parameters  of  the  lesions.  With 
the  a  pr/on  knowledge  of  lesion  structure  initially  obtained  by 
ultrasound,  the  tissue  volume  can  be  approximately  seg¬ 
mented  into  a  lesion  region  and  a  background  region;  optical 
reconstruction  from  the  optical  measurements  can  be  per¬ 
formed  by  using  finer  imaging  voxel  size  corresponding  to 
the  lesion  region  and  coarse  voxel  size  corresponding  to  the 
background  region  [18].  Consequently,  the  image  recon¬ 
struction  is  well  defined  because  the  total  number  of  imaging 
voxels  with  unknown  optical  properties  is  reduced  signifi¬ 
cantly.  In  addition,  the  reconstruction  is  less  sensitive  to 
noise  because  the  weight  matrix  is  appropriately  scaled  for 
inversion,  and  convergence  can  be  achieved  within  a  small 
number  of  iterations. 

In  this  paper,  we  report  our  initial  clinical  results  of  using  " 
the  combined  technique  to  image  and  characterize  ultrasdn- 
ically  detected  breast  lesions.  A  total  of  18  patients  (19 
cases)  of  invasive  breast  carcinomas'(2  cases),  combined 
fibroadenoma  and  fibrocystic  changes  with  scattered  foci  of; 
lobular  neoplasia/lobular  carcinoma  in  situ  (noninvasive  car-  ^; 
cinoma;  1  case),  juvenile  atypical  ductal  hyperplasia  (1 
case),  and  fibroadenomas  and  fibrocystic  changes  (15 
cases)  are  reported  here.  The  diameters  of  these  lesions 
are  on  the  order  of  1  cm.  For  the  first  time,  we  report  high 
sensitivity  of  optical  tomography  in  diagnosing  early-stage 
invasive  cancers  and  in  differentiating  ttiem  from  small 
benign  lesions.  We  also  show  a  highly  localized  pattern  of 
angiogenesis  distribution  of  small  invasive  cancers.  Our 
initial  results  are  very  encouraging  and  suggest  that  the 
discrimination  of  benign  and  malignant  breast  masses  could 
be  enhanced  by  this  type  of  achievable  optical  quantification 
coupled  with  ultrasound  localization. 


Materials  and  Methods 

The  major  technical  aspects  of  our  combined  imaging  tech¬ 
nique  have  been  described  in  detail  previously  [16].  Briefly, 
the  NIR  system  consists  of  12  pairs  of  dual-wavelength  (780 
and  830  nm)  laser  diodes,  which  are  used  as  light  sources, 
and  their  outputs  are  amplitude-modulated  at  140  MHz.  For 
each  source  pair,  the  outputs  from  the  two  laser  diodes  at 
both  wavelengths  were  coupled  to  the  tissues  through  a  200- 
pm  optical  fiber.  On  the  receiver  side,  eight  photomultiplier 
tubes  (PMTs)  were  used  to  detect  diffusely  scattered  light 
from  the  tissues,  and  eight  optical  fibers  (3  mm  in  diameter) 


were  used  to  couple  detected  diffusive  light  from  the  tissues 
to  the  PMTs.  The  high-frequency  radiofrequency  (RF)  sig¬ 
nals  were  amplified  and  heterodyned  to  20  kHz.  The  hetero¬ 
dyned  signals  were  further  amplified  and  band  pass-filtered 
at  20  kHz.  A  reference  signal  of  20  kHz  was  also  generated 
by  directly  mixing  the  detected  RF  signals  with  the  RF  signal 
generated  from  the  oscillator.  The  reference  signal  was 
necessary  for  retrieving  phase  shifts.  Eight  detection  signals 
and  one  reference  were  sampled  and  acquired  into  the  PC 
simultaneously.  The  entire  data  acquisition  took  about  3  to  4 
seconds,  which  was  fast  enough  for  acquiring  data  from 
patients. 

A  picture  of  our  combined  probe  and  imaging  system 
used  in  the  reported  clinical  studies  is  shown  in  Figure  la, 
and  the  probe  configuration  and  optical  sensor  distribution 
are  shown  in  Figure  1  b.  Both  amplitude  and  phase  at  each 
source-detector  pair  are  obtained  and  the  resulting  total 
number  of  measurements  is  12  x  8  x  2  =  192.  The  com¬ 
bined  probe  is  made  of  a  black  plastic  plate  10  cm  in 
diameter:  therefore,  a  semi-infinite  boundary  condition  can 
be  used  for  NIR  measurement  geometry.  The  amplitudes 
and  phases  acquired  at  all  source -detector  pairs  at  the 
normal  contralateral  breast  were  used  to  calculate  the 


(b) 


Figure  1.  (a)  Picture  of  the  combined  probe  and  frequency  domain  imager, 
(b)  Sensor  distribution  of  the  combined  probe.  Smaller  circles  in  (b)  are 
optical  source  fibers  and  big  circles  are  detector  fibers.  A  commercial 
ultrasound  probe  is  located  at  the  center  of  the  combined  probe,  end  the 
optical  source  and  detector  fibers  are  distributed  at  the  periphery  of  the 
ultrasound  probe. 
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background  absorption  coefficient  and  the  reduced 
scattering  coefficient  fis' 

The  details  of  our  dual-mesh  optical  imaging  reconstruc¬ 
tion  algorithm  have  been  described  in  Ref.  [18].  Briefly,  the 
entire  tissue  volume  Is  segmented  based  on  initial  coregis¬ 
tered  ultrasound  measurements  into  a  lesion  region  L  and  a 
background  region  B.  A  modified  Born  approximation  is  used 
to  relate  the  scattered  field  Usd(rsij,rdt,co)  measured  at  the 
source  (s)  and  detector  (d)  pair  /  to  absorption  variations 
(r')  in  each  volume  element  of  two  regions  within  the  sample. 
The  matrix  form  of  image  reconstruction  is  given  by: 

[t^Sd]A,x/  =  [WL,WB]MxNlML,/WB]Ll 

where  Wl  and  Wb  are  weight  matrices  for  lesion  and 
background  regions,  respectively;  and  [Ml]  =  A/ia(r') 

dV', . A/ia(r')d®r'  and  [Mb]  =  A//a(f')d®r', . . . 

. . .  /^  A/ifl(r')d®r')  are  total  absorption  distributions  of  lesion 
and  background  regions,  respectively.  The  weight  matrices 
are  calculated  based  on  the  background  pa  and  ps'  measure¬ 
ments  obtained  frorh  the  normal  contralateral  breast. 

Instead  of  reconstructing  the  Apa  distribution  directly,  as  is 
done  in  the  standard  Born  approximation,  the  total  absorp¬ 
tion  distribution  M\s  reconstructed  and  the  total  is  divided  by 
different  voxel  sizes  of  lesions  and  background  tissues  to 
obtain  the  A|Ja  distribution.  By  choosing  a  finer  grid  for  lesion; 
and  a  coarse  grid  for  background  tissues,  we  can  maintaih 
the  total  number  of  voxels  with  unknown  absorption  on  the  r 
same  scale  of  the  total  measurements.  a  result,  the 
inverse  problem  is  less  underdetermined  ■  In  addition,  be¬ 
cause  the  lesion  absorption  coefficient  Is  higher  than  that  of 
background  tissues,  in  general,  the  total  absorption  of  the. 
lesion  over  a  smaller  voxel  is  on  the  same  scale  as  thb  total  ■ 
absorption  of  the  background  over  a  bigger  voxel.  Therefore, 
the  matrix  [Mu, Mb]  is  appropriately  scaled  for  inversion.  Iri  all 
experiments,  a  finer  grid  of  0.5  x  0.5  x  0.5  (cm®)  and  a 
coarse  grid  of  1.5  x  1.5  x  1.5  (cm®)  were  chosen  for  the 
lesion  region  and  the  background  tissues,  respectively.  The 
total  reconstruction  volume  was  chosen  to  be  9  x  9  x  4 
cm®.  The  reconstruction  is  formulated  as  a  least  square 
problem  and  the  unknown  distribution  M  is  iteratively  calcu¬ 
lated  using  the  standard  conjugate  gradient  method. 

The  perturbations  for  both  wavelengths  used  to  calculate 
absorption  maps  were  normalized  as  L/sd(/'si,/di,<»)  = 
t4(rsi,  fdi,  co),  where  UL(rsi,rdi,ffl)  and  LM 
(fsu  were  optical  measurements  obtained  at  the  lesion 
Mregion  and  the  normal  region  of  the  contralateral  breast, 
and  UB{rs\MM  was  the  calculated  incident  field  using  the 
measured  background.  This  procedure  cancels  the  un¬ 
known  system  gains  associated  with  different  sources  and 
detectors  as  well  as  electronic  channels. 

The  commercial  one-dimensional  (1-D)  ultrasound  probe 
acquires  two-dimensional  (2-D)  ultrasound  images  in  the  x- 
zplane  (zis  the  propagation  direction)  and  the  2-D  NIR  probe 
provides  three-dimensional  (3-D)  optical  measurements  for 
3-D  Image  reconstruction.  Therefore,  at  each  location,  a  2-D 
ultrasound  image  is  coregistered  with  a  corresponding  set  of 
3-D  optical  measurements  in  the  x-z  plane.  However,  if  we 


approximate  a  lesion  as  an  ellipsoid,  we  are  able  to  estimate 
its  diameters  from  two  orthogonal  ultrasound  images.  The 
3-D  lesion  center  can  be  approximated  from  two  orthogonal 
2-D  ultrasound  images.  However,  two  sources  of  error  can 
affect  accurate  lesion  location  and  volume  estimation.  First, 
the  diameter  measurements  of  large  irregular  lesions  are 
inaccurate  because  lesion  boundaries  may  not  be  well 
defined  in  ultrasound  images.  Second,  the  target  boundaries 
seen  by  different  modalities  may  be  different  due  to  different 
contrast  mechanisms.  To  account  for  these  errors,  we  use  a 
lesion  center  estimated  from  the  2-D  coregistered  ultrasound 
image  and  use  much  larger  diameters  in  both  spatial  dimen¬ 
sions  than  ultrasound-measured  ones  to  highlight  the  lesion 
region  for  fine-mesh  optica!  reconstruction.  We  have  found 
from  experiments  that  the  abovementioned  measurement 
inaccuracies  have  little  effects  on  reconstructed  optical  prop¬ 
erties  as- long  as  the  lesion  depth  is  measured  correctly  and 
the  total  number  of  unknown  voxels  is  of  the  same  order 
as  the  total  number  of  measurements.  The  lesion  depth  z 
and  the  lesion  boundaries  in  the  z  direction  can  be  estimated 
reasonably  well  from  2-D  coregistered  ultrasound  using 
normal  tissue  structure  patterns. 

Clinical  studies  were  performed  at  the  University  of 
Connecticut  Health  Center  (UCHC;  Farmington,  CT).  The 
UCHC  IRB  committee  approved  the  human  subject  proto- 
V  col.  Patients  with  palpable  and  nonpalpable  masses  that 
were  visible  on  clinical  ultrasound  and  who  were  scheduled 
for  biopsy  were  enrolled  as  research  subjects.  These 
subjects  were  scanned  with  the  combined  probe,  and 
ultrasound  images  and  optical  measurements  were  ac¬ 
quired  at  multiple  locations  including  the  lesion  region 
..  scanned  at  two  orthogonal  positions,  a  normal  region  of 
'4the  same  breast  if  the  breast  was  large,  and  a  normal 
symmetric  region  of  the  contralateral  breast  also  scanned 
at  two  orthogonal  positions. 


Results 

Examples  of  Invasive  Cancers 

Figure  2a  shows  a  gray  scale  ultrasound  image  of  a 
nonpalpable  lesion  of  a  55-year-old  woman.  The  ultrasound 
showed  a  nodular  mass  with  internal  echoes  and  the  lesion 
was  considered  suspicious.  The  estimated  lesion  diameter 
measured  from  two  orthogonal  ultrasound  Images  was  8  mm. 
An  ultrasound-guided  core  needle  biopsy  was  recommended 
and  the  biopsy  result  revealed  that  the  lesion  was  Intraductal 
and  infiltrating  ductal  carcinoma  (nuclear  grade  11,  histologic 
grade  111).  The  cancer,  once  removed  from  the  breast, 
measured  1  cm  in  greatest  diameter,  and  was  composed 
predominantly  of  invasive  carcinoma  (>95%),  extending  to 
surgical  margins.  Negative  margins  were  achieved  upon 
reexcision. 

The  average  tissue  background  absorption  coefficientsp^ 
and  the  reduced  scattering  coefficients  pis'  at  780  and 
830  nm  were  measured  as  =  0.02  cm"'\  = 

8.2  cm"\  ixf®®  =  0.04  cm"\  and  =  6.0  cm"\  respec¬ 
tively.  The  initial  estimates  of  the  lesion  center  and  the 
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Figure  2.  {a)  shows  a  gray  scale  ultrasound  image  of  a  nonpalpable  lesion  of  a  55’year-old  woman.  The  lesion  pointed  by  the  arrow  was  located  at  the  4  o’clock 
position  of  the  right  breast  measuring  8  mm  in  diameter.  Ultrasound  showed  nodular  mass  with  interna!  echoes  and  the  lesion  was  considered  suspicious,  (b)  and 
(c)  are  optical  absorption  maps  at  780  and  830  nm,  respectively,  x  and y  axes  are  spatial  dimensions  of  9  9  cm.  The  slices  starting  from  left  to  right  correspond  to 

absorption  maps  from  0.7  cm  underneath  the  skin  surface  to  the  chest  wall  in  0.5-cm  spacing,  (d)  is  the  total  hemoglobin  concentration  and  the  vertical  scale  is 
in  micromoles. 


diameter  from  coregistered  ultrasound  were  (0,  0.3,  2.3)  cm 
and  8  mm.  To  account  for  errors  in  lesion  geometry  estimate 
and  for  differences  in  optical  contrast,  we  used  a  6-cm 
diameter  in  both  x  and  y  spatial  dimensions  at  the  center  of 
(0,  0.3,  2.3)  cm  for  fine-mesh  optical  reconstruction.  The 
lesion  diameter  in  depth  is  estimated  as  1  cm  by  comparing 
the  top  and  the  bottom  of  lesion  boundaries  from  normal 
tissue  patterns.  The  white  arrow  arrays  on  the  top  and  bottom 
of  the  lesion  In  Figure  2a  point  to  normal  tissue  boundaries. 

The  optical  absorption  maps  at  both  wavelengths  are 
shown  in  Figure  2,  b  and  c,  respectively.  In  both  Figure  2, 
b  and  c,  the  first  slice  is  0.7  cm  deep  Into  the  breast  tissue 
from  the  skin  surface  and  the  last  slice  is  closer  to  the  chest 


wall.  The  spacing  of  the  slices  is  0.5  cm.  The  horizontal  and 
vertical  axes  of  each  slice  are  spatial  x  and  y  dimensions  of  9 
cm  in  size.  The  lesion  is  well  resolved  in  slice  5  and  has 
shown  much  larger  spatial  extension  at  830  nm  than  that  at 
780  nm.  The  measured  maximum  absorption  coefficients  are 
=  0.24  cm"''  and  =  0.28  cm"'* ,  respectively,  and  the 

absorption  maxima  at  both  wavelengths  are  located  at  (0, 
1 .0, 2.3)  cm,  which  is  close  to  the  lesion  center  measured  by 
ultrasound. 

By  assuming  that  the  major  chromophores  are  oxygen¬ 
ated  (oxyHb)  and  deoxygenated  (deoxyHb)  hemoglobin 
molecules  in  the  wavelength  range  studied,  we  can  esti¬ 
mate  the  distribution  of  total  hemoglobin  concentration  as 
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shown  in  Figure  2d  The  extinction  coefficients  used  for  cal¬ 
culating  oxyHb  and  deoxyHb  concentrations  were  ^Hb®°  = 
2.54,  fiHbof®  =  1.70,  =  1.80,  eHbof*^  =  2.42,  obtained 

from  Ref.  [19]  in  a  natural  logarithm  scale  with  units  of 
inverse  millimoles  times  inverse  centimeters.  The  measured 
maximum  total  hemoglobin  concentration  for  lesions  is 
122  pumol,  and  the  measured  average  background  hemo¬ 
globin  concentration  is  14  \irr\o\.  The  spatial  extent  of  the 
lesion’s  angiogenesis  is  measured  from  the  full  width  at  half 
maximum  (FWHM)  of  the  total  hemoglobin  map  and  it  is 
estimated  to  be  2.0  cm.  This  number  is  about  two  and  half 
times  larger  than  the  8  mm  diameter  measured  by  ultra¬ 
sound  and  suggests  that  optical  contrasts  extend  well 
beyond  the  cancer  periphery  due  to  angiogenesis. 

Because  our  hand-held  probe  can  be  easily  rotated  or 
translated,  we  have  acquired  at  least  three  coregistered 


ultrasound  and  NIR  data  sets  at  the  lesion  location  for  each 
patient  and  have  reconstructed  corresponding  optical  ab¬ 
sorption  maps  as  well  as  the  total  hemoglobin  concentration 
distribution  under  the  coregistered  ultrasound  guidance.  The 
average  maximum  total  hemoglobin  concentration  at  the 
cancer  region  is  115  ^imol  (±27.6  \imo\)  and  the  average 
background  total  hemoglobin  concentration  Is  13  ^.mol  (±1.5 
^jimol).  The  large  standard  deviation  at  the  cancer  region  is 
likely  related  to  the  probe  compression  at  different  spatial 
locations  and  elastic  properties  of  the  blood  vessels  in  the 
cancer  region. 

Another  example  was  obtained  from  a  56-year-old 
woman  who  had  a  solid  mass  with  internal  ultrasound 
echoes  measuring  9  mm  in  size  (Figure  3a).  The  lesion 
was  considered  suspicious,  and  ultrasound-guided  core 
needle  biopsy  revealed  that  the  lesion  was  in  situ  and 
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Figure  3.  This  example  was  obtained  from  a  56-year-old  woman  who  had  nonpalpable  lesion  located  at  the  10  o’clock  position  of  the  left  breast.  Ultrasound 
showed  a  solid  mass  (a)  with  internal  echoes  measuring  9  mm  in  size  and  the  lesion  was  considered  suspicious,  (b)  and  (c)  are  optical  absorption  maps  at  780  and 
830  nm,  respectively.  The  slices  start  from  0.4  cm  underneath  the  skin  surface  to  the  chest  wall  in  0.5-cm  spacing,  (d)  is  the  total  hemoglobin  concentration. 
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Figure  4.  (a)  shoves  a  solid  lesion  at  the  2  o’clock  position  measuring  1.1  cm  in  diameter  in  ultrasound  image.  The  lesion  was  diagnosed  as  having  benign 
appearance  because  of  well-defined  boundaries.  Ultrasound-guided  core  biopsy  revealed  scattered  foci  of  lobular  neoplasia,  carcinoma  in  situ  but  not  invasive,  and 
fibroadenoma  and  fibrocystic  changes  with  associated  microcalcifications,  (b)  and  (c)  are  optical  absorption  maps  at  780  and  830  nm,  respectively.  The  slices  start 
from  0.5  cm  underneath  the  skin  surface  to  the  chest  wall  in  0.5-cm  spacing,  (d)  is  the  total  hemoglobin  concentration. 


Invasive  ductal  carcinoma  (nuclear  grade  il,  histologic 
grade  II).  The  tumor,  once  removed  from  the  breast, 
measured  1 .5  cm  in  greatest  diameter,  and  was  composed 
predominantly  of  invasive  carcinoma  (>80%),  extending  to 
the  inferior/anterior  surgical  margin.  Negative  inferior  mar¬ 
gin  was  achieved  upon  reexcIsion. 

The  average  tissue  background  absorption  coefficient  ^ 
and  the  reduced  scattering  coefficient  \i's  at  780  and  830  nm 
were  measured  as  |.La®°  =  0.04  cm"\  =  4.7  cm"\  = 
0.03  cm“\  and  =  5.8  cm“\  respectively.  The  initial 


estimate  of  the  lesion  center  and  diameters  from  coregis¬ 
tered  ultrasound  were  (0,  0.6,  1.9)  cm  and  9  mm,  respec¬ 
tively.  Again,  we  used  a  6-cm  diameter  in  both  x  and  y  spatial 
dimensions  at  the  center  of  (0,  0.6, 1 .9)  for  fine-mesh  optical 
reconstruction.  The  optical  absorption  maps  at  both  wave¬ 
lengths  are  shown  in  Figure  3,  b  and  c,  respectively.  In 
Figure  3,  b  and  c,  the  first  slice  is  0.4  cm  deep  into  the  breast 
tissues  from  the  skin  surface  and  the  last  slice  is  closer  to  the 
chest  wall.  The  spacing  between  the  slices  is  0.5  cm.  This 
lesion  is  well  resolved  in  slice  #4  and  has  shown  much  larger 
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spatial  extension  at  780  nm  than  that  at  830  nm.  The 
measured  maximum  absorption  coefficients  are  =  0 
.29  cm"''  and  =  0.25  cm“\  respectively.  The  absorption 
maxima  at  both  wavelengths  are  located  at  (0  1.0  1 .9000) 
cm,  which  is  very  close  to  the  lesion  center  measured  by 
ultrasound.  The  calculated  maximum  total  hemoglobin  con¬ 
centration  of  the  lesion  Is  128  p,mol,  and  the  background 
concentration  is  24  p,mol.  The  measured  FWHM  of  the  total 
hemoglobin  map  is  3.5  cm,  which  is  more  than  three  times 
larger  than  the  9  mm  diameter  measured  by  ultrasound  and 
again  suggests  that  optical  contrasts  extend  well  beyond  the 
cancer  periphery  due  to  angiogenesis.  The  average  maxi¬ 


(a) 

830  nm 


(c) 


mum  total  hemoglobin  concentration  calculated  from  three 
NIR  data  sets  is  123  fxmol  (±6.2  p,mol)  and  the  average  back¬ 
ground  hemoglobin  concentration  is  24  [xmol  (±0.5  jimoi). 
The  variation  of  total  hemoglobin  concentrations  acquired  at 
slightly  different  spatial  probe  positions  is  much  smaller  than 
that  In  the  previous  case. 


An  Example  of  Combined  Fibroadenoma  and  Fibrocystic 
Changes  with  Scattered  Foci  of  Lobular  Neoplasia/Lobular 
Carcinoma  In  Situ  (Noninvasive  Carcinoma) 

This  patient  had  a  solid  lesion  measuring  1.1  cm  In 
diameter  as  shown  in  Figure  4a  ultrasound  image.  The  lesion 


780  nm 


(b) 


(d) 


Figure  5.  (a)  shows  an  ultrasound  image  of  a  solid  hypoechogenic  lesion  located  at  9  to  10  o’clock  position  of  a  SV-year-old  woman's  breast  The  diagnosis  was 
that  the  lesion  could  represent  fibroadenoma;  however,  there  was  a  chance  that  the  lesion  could  be  a  carcinoma  because  of  the  irregular  shape  and  posterior 
shadow  seen  by  ultrasound,  (b)  and  (c)  are  optical  absorption  maps  at  780  and  830  nm,  respectively.  The  slices  start  from  0.5  cm  underneath  the  skin  surface  to 
the  chest  wall  in  0.5-cm  spacing,  (d)  is  the  total  hemoglobin  concentration. 


Neoplasia  •  Vo!.  5,  No.  5,  2003 


BEST  AVAILABLE  COPY 


8  Ultrasound-Guided  Imaging  of  Breast  Lesions  Zhu  et  al. 


Table  1.  Measured  Parameters  of  19  Cases. 


Reference  Biopsy  Results 
Number 


Lesion  Size  ROI 

(Ultrasound)  [cm] 


13 

Invasive  cancer 

0.9  X  0.6 

6x6x1 

10  X  10  X  1 

25 

Invasive  cancer 

0.8  X  0.5 

6x6x1 

10  X  10  X  1 

32 

Fibrocystic  and  lobular 
neoplasia/carcinoma 

1.1  X  1.1 

10  X  10  X  1.2 

11 

Fibroadenoma 

1.2  X  0.8 

8x8x1 

8 

Fibroadenoma 

2.2  X  1.3 

8  X  8  X  1.6 

9 

Fibroadenoma 

0.9  X  0.7 

8x8x1 

7 

Fibroadenoma  and 
sclerosing  adenosis 
with  extensive 
calcification 

1.0  X  0.6 

10  X  10  X  1 

17 

Fibrocystic  changes 

0.6  X  0.7 

10  X  10  X  1 

20 

Fibroadenoma 

1.2  X  0,6 

10  X  10  X  1 

30 

1:  Fibrocystic  change 

1.7  X  1.1 

8  X  8  X  1.3 

30 

2:  Sclerosing  adenosis 
and  epithelial 
hyperplasia 
without  atypical 

1.4  X  1.2 

8  X  8  X  1.3 

22 

Fibroadenoma 

1.0* 

10  X  10  X  1.5 

36 

Fibroadenoma 

1.9  X  0.9 

9x9x1 

28 

Fibroadenoma 

1.9  X  1.1 

9  X  9  X  1.3 

33 

Fibroadenoma 

1,2  X  0.7 

10  X  10  X  1 

37 

Fibroadenoma 

0.8  X  0.6 

8x8x1 

35 

Fibroadenoma 

1.3  X  1.3 

9x9x2 

38 

Hyalin  ized  fibroadenoma 

0.8  X  0.4 

6x6x1 

29 

Atypical  ductal  hyperplasia 

1.1  X  0.8 

9  X  9  X  t 

*See  by  X-ray. 


was  diagnosed  as  having  a  benign  appearance;  however, 
the  biopsy  result  revealed  combined  fibrocystic  changes 
associated  with  microcalcifications  and  scattered  foci  of 
lobular  neoplasia/lobular  carcinoma  In  situ  with  pagetoid ' 
extension  along  ducts,  but  no  invasive  carcinoma.  Alsb 
present  was  a  fibroadenoma  involved  by  lobular  neoplasia. 


140 


120-- 

benign 

fibroadenomas 
100--  Ahypeiplasia 


meanD  1.05  cm 


combined 
fibroadenoma 
and  fibrocystic 
changes  with 
&neoplasia/carci 
noma  in  situ 


Dl.l  cm 


invasive 

cancers 


mean  D  0.85  cm 


FWHM 
(NIR)  [cm] 

Max  pe  [cm“^ 
780  nm 

],  Max  Me  [cm 
830  nm 

\  Max  Total 
Hb[M.M] 

Average 

Total 

HbbiM] 

3.5 

0.28 

0.24 

123 

87 

3.5 

0.27 

0.24 

122 

88 

2.0 

0.22 

0.26 

115 

84 

2.1 

0.22 

0.26 

115 

84 

6.1 

0.12 

0.08 

48 

30 

4.9 

0.11 

0.10 

49 

37 

2.7 

0.07 

0.03 

24 

16 

4.9 

0.15 

0.18 

77 

52 

3.8 

0.16 

0.11 

59 

41 

5.8 

0.16 

0.16 

76 

53 

5.7 

0.14 

0.14 

67 

45 

4.0 

..  0.10 

0.15 

59 

39 

4.0 

0.15 

0.23 

88 

57 

3.7 

0.15 

0.15 

72 

46 

4.7 

0.13 

■  j  0.18 

73 

50 

.  0.14 

0.i6 

'  71 

45 

7.0 

0.13- 

0.13 

61 

44 

6.0 

0.16 

0.18 

80 

53 

4.3 

^  0.14  ‘ 

0.15 

59 

41 

2:5 

0.17 

0.22 

90 

67 

3:7 

0.10 

0.17 

64 

44 

The  average  background  tissue  absorption  coefficients 
and  reduced  scattering  coefficients  at  both  wavelengths 
were  measured  as  =  0.04  cm“\  =  3-3  cm“^, 
^  cm"\  and  ns'“°  =  3.5  cm-\  respectively.  The 
initial  estimate  of  the  lesion  center  and  the  diameter  from 
ultrasound  were  (0, 0.9, 1.4)  and  1.1  cm.  The  reconstructed 
absorption  maps  as  well  as  total  hemoglobin  concentration 
distribution  are  shown  in  Figure  4,  b-d,  respectively.  The 
measured  maximum  absorption  coefficients  are  =  0.13 
cm"^  and  =  0.09  cm"\  respectively,  which  are  about 
half  of  those  measured  in  invasive  cancer  cases.  In  addition, 
the  absorption  coefficient  distributions  are  diffused  at  both 
wavelengths  and  no  resolvable  or  localized  lesions  are  seen. 
The  calculated  maximum  total  hemoglobin  concentration  is 
53  M,mol,  which  is  less  than  half  of  those  in  invasive  cancer 
cases.  The  calculated  background  value  is  13  fimol.  The 
average  maximum  hemoglobin  concentration  calculated 
frorn  four  reconstructed  NIR  images  at  lesion  location  is 
48  pmol  (±2.8  ^.mol)  and  the  value  for  the  background  is 
12  M,mol  (±0.4  jimol).  This  example  suggests  that  tumor 
vascularization  may  not  be  developed  at  the  early  noninva- 
sive  stage  of  combined  fibrocystic  change  and  neoplasia/ 
carcinoma  In  situ. 


Figure  6.  The  average  maximum  total  hemoglobin  concentration  obtained 
from  benign  fibroadenomas  and  hyperplasia,  and  combined  fibroadenoma 
and  fibrocystic  change  with  naoplasla/carcinoma  in  situ  (noninvasive)  and 
Invasive  cancer  groups  are  shown  in  bars.  The  standard  deviations  of  three 
groups  are  also  provided.  The  average  sizes  of  lesions  of  the  three  groups 
are  1.05  (0.3),  1. 1,  and  0.9  (0.07)  cm,  respectively. 


Example  of  Fibroadenoma  and  Other  Benign  Lesions 
Figure  5a  shows  an  ultrasound  image  of  a  hypoechoic 
mass  of  a  37-year-old  woman.  The  diagnosis  was  that  the 
lesion  likely  was  a  fibroadenoma;  however,  there  was 
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concern  that  the  lesion  could  be  a  carcinoma  because  of 
the  irregular  shape  and  the  posterior  shadowing  seen  by 
ultrasound.  An  ultrasound-guided  core  breast  biopsy 
revealed  that  the  lesion  was  simply  a  fibroadenoma. 

Optical  absorption  maps  are  shown  in  Figure  5,  b  and  c, 
as  well  as  the  total  hemoglobin  distribution  in  Figure  5d. 
Compared  with  the  invasive  cancer  cases,  the  spatial  dis¬ 
tributions  of  the  absorption  coefficients  as  well  as  the  total 
hemoglobin  concentration  are  quite  diffused,  and  a  higher 
absorption  region  at  830  nm  corresponds  to  lesion  regions 
seen  by  ultrasound.  The  measured  maximum  lesion  ab¬ 
sorption  coefficients  at  780  and  830  nm  are  |4®°=0.10  cm“^ 
and  p|®°=0.12  cm''\  respectively,  and  these  values  are  less 
than  half  of  those  small  invasive  cancer  cases.  The  calcu¬ 
lated  maximum  total  hemoglobin  concentration  is  52  pmol 
and  the  background  hemoglobin  concentration  is  14  p,moI. 
The  calculated  average  maximum  total  hemoglobin  concen¬ 
tration  from  four  NIR  images  acquired  at  lesion  region  is  59 
^imol  (±5.1  p<mol)  and  the  average  background  is  16  p-mol 
(±1 .0  pjnol). 

Table  1  lists  all  the  measured  parameters  of  the  19  cases 
[from  left  to  right:  biopsy  result,  lesion  size  measured  by 
ultrasound  in  xand  zdimensions  (because  lesions  are  small, 
the  lesion  size  in  y  Is  similar  to  that  measured  in  the  x 
dimension),  region  of  interest  (ROI)  used  for  fine-mesh  NIR 
imaging,  FWHM  measured  from  NIR  imaging,  maximum 
absorption  coefficients  measured  at  both  wavelengths, 
maximum  and  average  total  hemoglobin  concentrations 
(the  average  is  computed  within  FWHM)].  For  the  two 
invasive  cancer  cases,  NIR  parameters  are  given  for  two 
different  ROIs  and  the  results  show  that  the  choice  of  ROI 
has  negligible  effects  on  the  absorption  and  hemoglobin 
measurements. 

The  statistics  of  maximum  total  hemoglobin  concentration 
of  three  groups  obtained  from  the  19  cases  are  shown  in 
Figure  6.  The  benign  group  of  fibroadenoma  (15  cases)  and 
hyperplasia  (1  case)  has  an  average  of  67  iirnol  (±17.0 
p,moI),  the  combined  fibroadenoma  and  fibrocystic  change 
with  noninvasive  neoplasia/carcinoma  in  situ  case  has  a 
maximum  of  48  pmol,  and  the  invasive  cancer  group  of  two 
cases  has  shown  about  two-fold  greater  average  of  1 1 9  \ixr\o\ 
(±1.6  p^mol).  if  average  total  hemoglobin  concentration  is 
used,  the  values  are  46  (±11.3),  30,  and  86  (±2.12)  pmol  for 
the  three  corresponding  groups,  respectively.  The  malignant 
group  also  presents  about  two-fold  greater  average  hemo¬ 
globin  concentration  than  that  of  the  benign  group.  The 
average  sizes  of  lesions  of  the  three  groups  measured  by 
ultrasound  are  1.05  (±0.3),  1.1,  and  0.9  (±0.07)  cm,  respec¬ 
tively.  The  lesion  size  is  the  geometric  mean  of  diameters 
measured  in  x  and  z  dimensions. 

Because  our  study  is  In  its  early  clinical  trial  stage,  we 
do  not  Intend  to  provide  the  sensitivity  and  specificity  of  the 
combined  method  due  to  the  limited  sample  size  available. 
However,  our  initial  findings  are  very  encouraging  and 
suggest  that  early-stage  small  invasive  cancers  have  much 
higher  optical  contrast  than  benign  lesions  and  could  be 
diagnosed  with  high  specificity  possibly  due  to  tumor 
angiogenesis. 


Discussion 

As  described  earlier,  due  to  intense  light  scattering,  optical 
tomography  alone  has  not  been  widely  used  In  clinical 
studies.  Data  in  the  published  literature  have  been  limited 
to  feasibility  studies  or  case  reports.  Pogue  et  al.  [9]  reported 
pilot  results  of  one  invasive  ductal  carcinoma  1  cm  in  size 
and  one  benign  fibroadenoma  3  cm  In  size.  The  reported 
maximum  total  hemoglobin  concentrations  were  68  iimol  for 
the  cancer  case  and  55  pmol  for  the  benign  case.  Although 
the  invasive  cancer  size  reported  by  this  group  is  comparable 
to  ours,  the  system  they  used  can  only  acquire  NIR  data  from 
a  ring  area  with  optical  sources  and  detectors  deployed 
around  the  breast  (2-D  Imager).  Therefore,  incomplete  infor¬ 
mation  could  lead  to  smaller  reconstructed  absorption  coef¬ 
ficients  and  total  hemoglobin  concentrations  than  those 
reported  here  using  our  3-D  NIR  imager.  In  addition,  an 
investigation  with  phantoms  has  shown  that  NIR  alone.  In 
general,  reconstructs  lower  absorption  coefficients  and 
therefore  lower  total  hemoglobin  concentrations  than  true 
values  [20].  This  can  be  seen  from  two  studies  of  the  same 
data  reported  in  Refs.  [7,13].  Authors  In  Ref.  [7]  reported  an 
average  of  35  ^imol  totaf  hemoglobin  concentration  of  a  2-cm 
ductal  carcinoma  in  situ  by  using  NIR  measurements  alone. 
After  reprocessing  the  same  NIR  measurement  data  using 
..  an  approximate  lesion  depth  obtained  from  a  separate 
ultrasound  image,  the  authors  reported  that  the  calculated 
average  total  hemoglobin  concentration  was  Increased  to  67 
pmol[13]. 

The  reported  small  early-stage  invasive  cancers  appear 
isolated  and  are  well  resolved  from  background  tissues  in 
optical  absorption  maps  as  well  as  in  total  hemoglobin 
distributions.  However,  the  combined  fibroadenoma  and 
^■fibrocystic  change  with  scattered  foci  of  lobular  neoplasia/ 
lobular  carcinoma  in  situ  case  showed  no  difference  In 
optical  absorption  and  hemoglobin  concentration  than  those 
obtained  from  fibroadenoma  cases.  This  suggests  that  opti¬ 
cal  tomography  may  not  be  sensitive  to  early-stage  mixed 
benign  changes  and  noninvasive  neoplasia/carcinoma  in  situ 
because  tumor  neovascularization  has  not  been  developed. 
This  result  is  consistent  with  magnetic  resonance  imaging 
(MRI)  findings  on  the  low  detection  sensitivity  of  carcinoma  in 
situ  [21-23].  However,  optical  imaging  may  be  a  valuable 
tool  for  monitoring  the  development  or  transition  of  lesions 
from  noninvasive  to  early  Invasive  stage.  Certainly,  more 
cases  are  needed  to  validate  the  observations  reported  here. 
For  larger  cancers,  highly  heterogeneous  wavelength-de¬ 
pendent  optical  absorption  distributions  and  total  hemoglobin 
distribution  have  been  observed  in  four  cases  [1 8],  and  these 
distributions  could  provide  valuable  information  for  monitor¬ 
ing  and  assessing  cancer  therapy  under  treatment.  Current¬ 
ly,  we  are  pursuing  research  along  this  line  and  more  cases 
will  be  reported  in  the  future. 

In  principle,  the  distribution  of  oxygenation  saturation  can 
be  estimated  as  S  =  oxyHb/(oxyHb+deoxyHb),  with  deox- 
yHb  and  oxyHb  distributions  calculated  from  absorption 
maps  at  the  two  wavelengths  of  780  and  830  nm.  However, 
because  background  tissues  mainly  consist  of  water  and 
lipid  and  these  two  chromophores  contribute  to  the  total 
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concern  that  the  lesion  could  be  a  carcinoma  because  of 
the  irregular  shape  and  the  posterior  shadowing  seen  by 
ultrasound.  An  ultrasound-guided  core  breast  biopsy 
revealed  that  the  lesion  was  simply  a  fibroadenoma. 

Optical  absorption  maps  are  shown  in  Figure  5,  b  and  c, 
as  well  as  the  total  hemoglobin  distribution  in  Figure  5d 
Compared  with  the  invasive  cancer  cases,  the  spatial  dis¬ 
tributions  of  the  absorption  coefficients  as  well  as  the  total 
hemoglobin  concentration  are  quite  diffused,  and  a  higher 
absorption  region  at  830  nm  corresponds  to  lesion  regions 
seen  by  ultrasound.  The  measured  maximum  lesion  ab¬ 
sorption  coefficients  at  780  and  830  nm  are  p|®°=0.10  cm”'* 
and  |j|®^=0.12  cm"\  respectively,  and  these  values  are  less 
than  half  of  those  small  invasive  cancer  cases.  The  calcu¬ 
lated  maximum  total  hemoglobin  concentration  is  52  pjnol 
and  the  background  hemoglobin  concentration  is  14  [imol 
The  calculated  average  maximum  total  hemoglobin  concen¬ 
tration  from  four  NIR  images  acquired  at  lesion  region  is  59 
\umo\  (±5.1  \imo\)  and  the  average  background  Is  16  \ixno\ 
(±1.0  ^mol). 

Table  1  lists  all  the  measured  parameters  of  the  1 9  cases 
[from  left  to  right:  biopsy  result,  lesion  size  measured  by 
ultrasound  in  xand  zdimensions  (because  lesions  are  small, 
the  lesion  size  In  y  is  similar  to  that  measured  in  the  x 
dimension),  region  of  Interest  (ROI)  used  for  fine-mesh  NIR 
imaging,  FWHM  measured  from  NIR  imaging,  maximum 
absorption  coefficients  measured  at  both  wavelengths, 
maximum  and  average  total  hemoglobin  concentrations 
(the  average  is  computed  within  FWHM)].  For  the  two 
invasive  cancer  cases,  NIR  parameters  are  given  for  two 
different  ROIs  and  the  results  show  that  the  choice  of  ROI 
has  negligible  effects  on  the  absorption  and  hemoglobin 
measurements. 

The  statistics  of  maximum  total  hemoglobin  concentration 
of  three  groups  obtained  from  the  19  cases  are  shown  in 
Figure  6.  The  benign  group  of  fibroadenoma  (15  cases)  and 
hyperplasia  (1  case)  has  an  average  of  67  ^mol  (±17.0 
pmol),  the  combined  fibroadenoma  and  fibrocystic  change 
with  noninvasive  neoplasia/carcinoma  in  situ  case  has  a 
maximum  of  48  pmol,  and  the  Invasive  cancer  group  of  two 
cases  has  shown  about  two-fold  greater  average  of  1 19  pmol 
(±1.6  |i,mol).  If  average  total  hemoglobin  concentration  is 
used,  the  values  are  46  (±1 1.3),  30,  and  86  (±2.12)  [jmol  for 
the  three  corresponding  groups,  respectively.  The  malignant 
group  also  presents  about  two-fold  greater  average  hemo¬ 
globin  concentration  than  that  of  the  benign  group.  The 
average  sizes  of  lesions  of  the  three  groups  measured  by 
ultrasound  are  1.05  (±0.3),  1.1,  and  0.9  (±0.07)  cm,  respec¬ 
tively.  The  lesion  size  is  the  geometric  mean  of  diameters 
measured  in  x  and  z  dimensions. 

Because  our  study  is  In  its  early  clinical  trial  stage,  we 
do  not  intend  to  provide  the  sensitivity  and  specificity  of  the 
combined  method  due  to  the  limited  sample  size  available. 
However,  our  Initial  findings  are  very  encouraging  and 
suggest  that  early-stage  small  invasive  cancers  have  much 
higher  optical  contrast  than  benign  lesions  and  could  be 
diagnosed  with  high  specificity  possibly  due  to  tumor 
angiogenesis. 


Discussion 

As  described  earlier,  due  to  intense  light  scattering,  optical 
tomography  alone  has  not  been  widely  used  in  clinical 
studies.  Data  in  the  published  literature  have  been  limited 
to  feasibility  studies  or  case  reports.  Pogue  et  al.  [9]  reported 
pilot  results  of  one  invasive  ductal  carcinoma  1  cm  in  size 
and  one  benign  fibroadenoma  3  cm  in  size.  The  reported 
maximum  total  hemoglobin  concentrations  were  68  iimol  for 
the  cancer  case  and  55  jimol  for  the  benign  case.  Although 
the  invasive  cancer  size  reported  by  this  group  is  comparable 
to  ours,  the  system  they  used  can  only  acquire  NIR  data  from 
a  ring  area  with  optica!  sources  and  detectors  deployed 
around  the  breast  (2-D  imager).  Therefore,  incomplete  infor¬ 
mation  could  lead  to  smaller  reconstructed  absorption  coef¬ 
ficients  and  total  hemoglobin  concentrations  than  those 
reported  here  using  our  3-D  NIR  Imager.  In  addition,  an 
investigation  with  phantoms  has  shown  that  NIR  alone,  In 
general,  reconstructs  lower  absorption  coefficients  and 
therefore  lower  total  hemoglobin  concentrations  than  true 
values  [20].  This  can  be  seen  from  two  studies  of  the  same 
data  reported  in  Refs.  [7,13].  Authors  in  Ref.  [7]  reported  an 
average  of  35  fxmol  totaf  hemoglobin  concentration  of  a  2-cm 
ductal  carcinoma  in  situ  by  using  NIR  measurements  alone. 
After  reprocessing  the  same  NIR  measurement  data  using 
an  approximate  lesion  depth  obtained  from  a  separate 
ultrasound  Image,  the  authors  reported  that  the  calculated 
average  total  hemoglobin  concentration  was  increased  to  67 
Mmol[13]. 

The  reported  small  early-stage  invasive  cancers  appear 
Isolated  and  are  well  resolved  from  background  tissues  In 
optical  absorption  maps  as  well  as  In  total  hemoglobin 
distributions.  However,  the  combined  fibroadenoma  and 
■  fibrocystic  change  with  scattered  foci  of  lobular  neoplasia/ 
lobular  carcinoma  in  situ  case  showed  no  difference  in 
optical  absorption  and  hemoglobin  concentration  than  those 
obtained  from  fibroadenoma  cases.  This  suggests  that  opti¬ 
cal  tomography  may  not  be  sensitive  to  early-stage  mixed 
benign  changes  and  noninvasive  neoplasia/carcinoma  in  situ 
because  tumor  neovascularization  has  not  been  developed. 
This  result  is  consistent  with  magnetic  resonance  imaging 
(MRI)  findings  on  the  low  detection  sensitivity  of  carcinoma  in 
situ  [21-23],  However,  optical  imaging  may  be  a  valuable 
tool  for  monitoring  the  development  or  transition  of  lesions 
from  noninvasive  to  early  invasive  stage.  Certainly,  more 
cases  are  needed  to  validate  the  observations  reported  here. 
For  larger  cancers,  highly  heterogeneous  wavelength-de- 
pendent  optical  absorption  distributions  and  total  hemoglobin 
distribution  have  been  observed  In  four  cases  [1 8],  and  these 
distributions  could  provide  valuable  information  for  monitor¬ 
ing  and  assessing  cancer  therapy  under  treatment.  Current¬ 
ly,  we  are  pursuing  research  along  this  line  and  more  cases 
will  be  reported  in  the  future. 

In  principle,  the  distribution  of  oxygenation  saturation  can 
be  estimated  as  S  =  oxyHb/(oxyHb4*deoxyHb),  with  deox- 
yHb  and  oxyHb  distributions  calculated  from  absorption 
maps  at  the  two  wavelengths  of  780  and  830  nm.  However, 
because  background  tissues  mainly  consist  of  water  and 
lipid  and  these  two  chromophores  contribute  to  the  total 
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absorption  estimate  as  well,  we  could  not  obtain  reasonable 
background  oxygenation  saturation  and  compare  it  with 
lesion  oxygenation  saturation.  Recently,  we  have  improved 
our  NIR  system  by  adding  another  wavelength  at  660  nm, 
which  may  allow  us  to  accurately  estimate  the  background 
oxygenation  saturation  and  to  compare  it  with  lesion  oxy¬ 
genation  saturation. 

If  we  reconstruct  the  lesion  area  only,  we  could  distribute 
partiai  perturbations  caused  by  background  to  lesion  and 
increase  calculated  lesion  absorption  and  therefore  hemo¬ 
globin  concentration.  However,  with  the  dual-mesh  scheme, 
we  reconstruct  the  entire  imaging  volume  instead  of  the 
lesion  area  only,  and  distribute  the  perturbations  to  both 
lesion  and  background.  We  did  phantom  experiments  using 
the  dual-mesh  scheme  and  obtained  absorption  coefficients, 
which  were  always  within  10%  of  the  true  values  depending 
on  the  phantom  contrasts  [24]. 

Two-dimensional  ultrasound  provides  fine  tissue  layer 
structures  in  depth  {2!)  direction  (Figure  2,  arrow  arrays), 
and  the  2-D  ultrasound  image  is  coregistered  with  optical 
data  in  z.  Therefore,  we  only  need  to  consider  possible 
angiogenesis  extensions  in  z  and  we  have  extended  the 
lesions  to  the  closest  normal  tissue  lines  in  z.  In  spatial 
dimensions,  we  have  some  uncertainty  in  another  dimension 
(3^  based  on  2-D  x-z  ultrasound  image,  and  we  need  to 
extend  the  lesion  region  to  a  larger  area  to  account  for  this  as ' 
well  as  for  possible  angiogenesis  extension. 

The  reported  optical  tomography  study  was  fised"  to  ; 
image  and  characterize  ultrasonically  detected  lesions.  How¬ 
ever,  in  one  benign  fibroadenoma  case,  the  lesion  was  not 
visible  in  ultrasound  but  was  seen  by  conventional  X-ray 
mammogram.  By  knowing  the  approximate  lesion  region- 
with  respect  to  the  nipple  location  from  the  patient’s ‘  mam-  : 
mogram,  we  used  a  fine  mesh  for  optical  reconstruction  in  a 
larger  region  of  9  x  9  x  1 .5  cm®,  and  identified  a  possible 
lesion  that  showed  slightly  higher  optical  absorption  than  that 
of  the  background.  We  believe  that  optical  tomography, 
assisted  by  conventional  mammography  and/or  ultrasound 
localization,  has  potential  as  a  screening  tool  to  identify  and 
characterize  malignancy. 

Our  initial  findings  provide  evidence  that  optical  tomogra¬ 
phy,  combined  with  ultrasound,  could  be  used  to  differentiate 
early-stage  small  invasive  breast  cancers  from  benign 
lesions.  Because  of  the  limited  sample  size  available,  we 
are  not  able  to  provide  sensitivity  and  specificity  results  in 
this  paper.  In  addition,  for  the  patients  studied,  four  fibroa¬ 
denoma  cases  and  one  intraductal  hyperplasia  were  exclud¬ 
ed  from  the  statistics  shown  in  Figure  6.  The  four 
fibroadenoma  cases  were  scanned  at  the  beginning  of  the 
study.  One  patient  was  scanned  after  core  biopsy  due  to 
scheduling  problem,  and  the  possible  blood  distribution 
change  due  to  biopsy  procedure  has  to  be  considered.  The 
other  three  young  patients  have  small  and  dense  breasts.' 
For  these  three  patients,  the  optical  amplitude  and  phase 
data  sets  were  highly  scattered  even  in  the  normal  contra¬ 
lateral  breast.  Therefore,  no  reliable  background  tissue 
absorption  and  scattering  coefficients  can  be  obtained. 
Three  possible  sources  for  these  findings  were  identified. 


First,  because  these  breasts  were  small  and  dense,  the  skin 
and  probe  contact  was  not  always  good.  Second,  a  thick  gel 
layer  used  for  coupling  the  ultrasound  transducer  with  the 
skin  can  serve  as  a  light  tunneling  medium  from  sources  to 
detectors.  This  portion  of  the  light  can  saturate  the  detectors 
and  give  false  readings.  Caution  was  taken  in  the  later 
studies  by  compressing  the  probe  harder  against  the  ex¬ 
amined  breast  to  ensure  good  probe-tissue  contact,  and  by 
placing  a  very  thin  gel  layer  underneath  the  ultrasound 
transducer  during  scanning.  The  third  source  are  tissue 
heterogeneities  of  the  young  dense  breasts.  However,  by 
carefully  removing  outliers  from  normal  breast  data,  we 
could  obtain  reasonably  good  reference  data  for  imaging. 
One  patient  with  intraductal  hyperplasia  data  was  excluded 
from  the  study.  This  patient  had  a  solid  lesion  located  at  the 
nipple  area  and  the  nipple  had  a  very  dark  color.  Four  NIR 
image  data  sets  obtained  at  the  lesion  area  had  the  same 
artifacts  and  the  artifact  consistently  showed  up  with  very 
high  absorptions  at  both  wavelengths.  When  the  lesion 
locations  were  changed  in  different  images  with  different 
probe  positions,  the  artifact  location  changed  very  little.  A 
similar  artifact  showed  up  in  the  two  reference  data  sets 
acquired  in  the  contralateral  normal  breast.  We  do  not  know 
the  exact  source  of  this  artifact  because  there  was  no  report 
that  this  patient  had  problems  at  the  contralateral  breast. 
Care  was  taken  in  the  later  studies  by  acquiring  at  least  four 
Teference  data  sets  at  the  normal  symmetric  region  of  the 
contralateral  breast,  as  well  as  at  normal  regions  of  the 
same  lesion  breast.  Reference  data  sets  were  always 
checked  for  normal  background  absorption  changes  before 
imaging  the  lesions. 


Summary 

Initial  clinical  results  from  the  use  of  optical  tomography 
combined  with  conventional  ultrasound  demonstrate  that 
there  is  a  huge  optical  contrast  between  early-stage  invasive 
cancers  and  benign  solid  lesions  due  to  angiogenesis.  An 
average  of  52  pmol  difference  was  obtained  between  two 
small  invasive  cancers  and  a  group  of  16  benign  solid 
lesions.  In  addition,  the  small  invasive  cancers  were  local¬ 
ized  well  in  absorption  maps  and  have  shown  wavelength- 
dependent  absorption  changes,  whereas  the  benign  lesions 
appeared  more  diffused  in  absorption  maps  and  have  shown 
relatively  wavelength-independent  absorption  changes.  The 
combined  fibrocystic  changes  and  noninvasive  neoplasia/ 
carcinoma  in  situ  case  did  not  show  a  significant  difference  in 
optical  absorption  and  total  hemoglobin  concentration  than 
those  of  benign  lesions.  This  suggests  that  the  angiogenesis 
may  not  be  developed  at  the  early  noninvasive  stage  of  the 
mixed  benign  changes  and  fibroadenoma  with  neoplasia/ 
carcinoma  in  situ.  However,  optical  tomography  may  have  a 
potential  role  in  monitoring  the  development  and/or  transition 
of  cancers  from  the  noninvasive  to  the  invasive  stage. 

Our  reported  initial  results  are  very  encouraging  and 
demonstrate  that  our  unique  approach,  which  combines 
optical  tomography  with  ultrasound,  has  great  potential  to 
detect  and  characterize  breast  lesions. 
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Abstract 

Purpose:  To  investigate  the  feasibility  of  optical  tomography  with  ultrasonography  (US)  localization 
in  differentiating  benign  from  malignant  breast  masses.  To  compare  optical  method  with  Doppler  US  in 
mapping  tumor  angiogenesis. 

Material  and  Methods:  Using  optical  tomography  with  US  localization,  we  have  overcome  the 
light  scattering  problem.  A  hand-held  probe  with  a  commercial  ultrasound  transducer  deployed  in  the 
middle  and  optical  sensors  distributed  at  the  periphery  is  used  for  simultaneously  acquiring  co-registered 
US  image  and  optical  measurements.  The  lesion  location  provided  by  US  is  used  to  segment  the 
imaging  volume  for  optical  imaging  reconstruction.  Two  optical  wavelengths  are  used  and  three- 
dimensional  (3D)  light  absorption  maps  at  these  wavelengths  are  obtained.  From  the  absorption  maps, 
lesion  total  hemoglobin  concentration  distribution,  which  is  a  direct  measure  of  tumor  angiogenesis,  can 
be  calculated. 

Results:  Initial  results  of  5  malignant  cancers  and  16  benign  lesions  have  shown  that  small  invasive 
cancers  can  be  well-resolved  in  3D  optical  absorption  maps.  The  average  maximum  total  hemoglobin 
concentrations  for  malignant  and  benign  groups  are  11 0.3 //moles  and  47.6  //moles,  respectively.  The 
malignant  group  presents  more  than  two-fold  greater  total  hemoglobin  concentration  than  that  of  the 
benign  group. 

Conclusion:  Optical  tomography  with  US  localization  has  a  great  potential  in  differentiating  benign 
from  malignant  breast  lesions.  The  optical  specificity  on  mapping  tumor  angiogenesis  is  much  higher 
than  Doppler  US. 

Index  terms:  US,  Optical  tomography.  Diagnosis,  Breast  Neoplasma 
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Introduction 

The  large  number  of  biopsies  performed  for  benign  breast  abnormalities  has  long  been  recognized  as  a 
serious  problem  (1).  The  use  of  US  as  an  adjunct  to  x-ray  mammography  decreased  the  number  of 
biopsies  by  enabling  reliable  identification  of  simple  cysts  (2-4)  from  solid  lesions.  However,  US 
features  were  not  yet  reliable  enough  to  determine  whether  biopsy  should  be  performed  on  a  solid  mass 
(5,6). 

Tumor  angiogenesis  or  neovascularity  is  known  to  be  critical  for  the  autonomous  growth  and  spread  of 
breast  cancers  (7,8).  Tumor  angeogenesis  is  a  complex  process  involving  both  the  incorporation  of 
existing  host  blood  vessels  into  the  tumor  and  the  creation  of  tumor  microvessels.  This  process  is 
moderated  by  tumor  angiogenesis  factor  (9).  In  principle,  the  altered  hemodynamics  that  accompany 
tumor  angiogenesis  provide  the  basis  for  discrimination  between  malignant  and  benign  masses  of  the 
breast  by  color  Doppler  sonography  (10).  However,  the  contribution  of  Doppler  US  to  differentiate 
benign  and  malignant  solid  masses  is  still  controversial  (11,12). 

Optical  tomography,  a  new  technique  that  employs  diffuse  light  in  near  infrared  spectrum,  has  recently 
been  under  intensive  investigations  in  pilot  clinical  studies  (13-17).  Optical  tomography  provides 
functional  images  of  tumor  angiogenesis  and  tumor  hypoxia.  If  a  single  optical  wavelength  is  used,  the 
optical  absorption  related  to  tumor  angiogenesis  can  be  measured.  If  appropriately  selected  multiple 
wavelengths  are  used,  the  optical  absorptions  at  these  wavelengths  can  be  measured  and  distributions  of 
oxyhemoglobin  (oxyHb)  and  deoxyhemoglobin  (deoxyHb)  of  tumors  can  be  deduced.  The  total 
hemoglobin  concentration  and  tumor  hypoxia  can  be  calculated  from  distributions  of  oxyHb  and 
deoxyHb  and  these  functional  parameters  are  highly  correlated  with  lesion  malignancy.  However, 
optical  tomography  when  used  alone  has  been  limited  to  feasibility  studies  in  the  past.  The  intense 
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scattering  caused  by  tissue  limits  the  imaging  resolution  as  well  as  location  certainty.  Novel 
combination  of  optical  tomography  with  other  modalities,  such  as  US  and  MRI,  has  shown  promising 
results  of  utilizing  full  advantages  of  light  to  characterize  malignancy  versus  benign  processes  (17-20). 
Flexible  light  guide  using  optical  fibers  makes  optical  imaging  compatible  with  many  other  imaging 
modalities  and  allows  for  simultaneous  imaging  under  identical  geometric  conditions.  Further  more,  the 
lesion  structure  information  provided  by  other  modalities  can  be  used  to  assist  optical  imaging 
reconstruction  and  therefore  reduce  the  location  uncertainty  and  improve  the  quantification  accuracy  of 
light. 

We  have  investigated  the  optical  tomography  using  a  priori  lesion  structural  information  provided  by  co¬ 
registered  US.  The  light  guide  using  optical  fibers  are  coupled  to  a  hand-held  probe  with  a  commercial 
ultrasound  transducer  deployed  in  the  middle.  Preliminary  results  of  19  solid  lesions  have  shown  that 
two  early  stage  invasive  cancers  have  shown  two-fold  high  total  hemoglobin  concentration  than  a  group 
of  benign  lesions  (19).  Preliminary  results  of  4  advanced  cancers  have  shown  heterogeneous  total 
hemoglobin  distributions  and  the  distributions  correlate  with  histological  micro-vessel  density  counts 
(20).  In  this  paper,  we  further  demonstrate  the  clinical  potential  of  optical  tomography  assisted  with  US 
localization  in  distinguishing  benign  from  malignant  lesions,  and  we  compare  optical  tomography  with 
Doppler  US  on  mapping  tumor  angiogenesis. 

Materials  and  Methods 

NIR  system 

A  hand-held  hybrid  probe  consists  of  a  commercial  ultrasound  transducer  located  at  the  middle  and 
near-infrared  source-detector  light  guides  (optical  fibers)  distributed  at  the  periphery  (see  Fig.l).  The 
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technical  aspects  of  the  NIR  imager  have  been  described  in  detail  previously  (21).  Briefly,  the  imager 
consists  of  12  pairs  of  dual  wavelength  (780nm  and  830nm)  laser  diodes,  which  are  used  as  light 
sources,  and  their  outputs  are  coupled  to  the  probe  through  optical  fibers.  On  the  receiving  side,  8 
photomultiplier  tubes  were  used  to  detect  diffusely  scattered  light  fi'om  the  tissue  and  8  optical  fibers 
were  used  to  couple  detected  light  to  the  PMTs.  The  laser  diodes’  outputs  were  amplitude  modulated  at 
140  MHz  and  the  detector  outputs  were  demodulated  to  20KHz.  Eight  detection  signals  and  one 
reference  were  amplified,  sampled  and  acquired  into  a  PC  simultaneously.  The  entire  data  acquisition 
took  about  3  to  4  seconds,  which  was  fast  enough  for  acquiring  data  from  patients. 

Clinical  protocol 

Clinical  studies  were  performed  at  the  Radiology  Department  of  Hartford  Hospital.  The  Hartford 
Hospital  IRB  committee  approved  the  human  subject  protocol.  Patients  who  were  scheduled  for 
ultrasound  guided  biopsy  were  recruited  to  the  study.  The  state-of-the  art  US  scanner,  Acuson  Sequoia, 
was  used  for  the  study.  Ultrasound  images  and  optical  measurements  were  acquired  simultaneously 
before  biopsy  procedures  at  multiple  locations  including  the  lesion  region,  a  normal  region  of  the  same 
breast,  and  a  normal  symmetric  region  of  the  contralateral  breast.  The  optical  data  acquired  at  normal 
region  were  used  as  reference  for  calculating  the  scattered  field  caused  by  lesions. 


Based  on  x-ray/US  or  MRI  images,  one  of  radiologists  (EC,  AC,  VH)  scored  each  biopsied  lesion  using 
BIRAD  as  possibly  benign,  suspicious,  and  highly  suspicious  for  malignancy.  For  Doppler  US 
measurements,  any  persistence  color  Doppler  signals  were  taken  to  represent  blood  vessels.  The  location 
of  the  vessels  with  respect  to  the  lesion  was  documented  as  peripheral  to  the  lesion,  within  the  lesion,  or 


both. 
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Optical  Imaging  Method 

The  details  of  our  dual-mesh  optical  imaging  reconstruction  algorithm  have  been  described  in  Ref  (19). 
Briefly,  the  NIR  reconstruction  takes  advantages  of  ultrasound  localization  of  lesions  and  segments  the 
imaging  volume  into  finer  grid  in  lesion  region  and  coarser  grid  in  non-lesion  regions.  A  modified  Bom 
approximation  is  used  to  relate  the  scattered  field  measured  at  the  source-detector  pairs  to  absorption 
variations  in  each  volume  element  of  two  regions  within  the  sample.  With  this  dual-mesh  scheme,  the 
inverse  optical  reconstruction  is  well-conditioned  and  converges  in  few  iterations. 

Results 

The  first  example  was  obtained  from  a  45-year  old  woman  who  had  a  suspicious  lesion  located  at  12  o’ 
clock  position  shown  on  MRI.  The  lesion  was  identified  sonographically  at  the  time  of  biopsy  as  ill- 
defined  slightly  heterogeneous,  mildly  hypoechoic  mass  of  7  mm  in  diameter.  Doppler  US  revealed 
large  blood  vessels  both  inside  and  at  peripheral  of  the  lesion  (Fig.2  (a)).  US  images  and  optical 
measurements  were  obtained  simultaneous  before  US  guided  biopsy  procedure.  Surgical  pathology 
reports  revealed  that  the  lesion  was  intraductal  and  invasive  mammary  duct  carcinoma  (nuclear  grade 
III,  histologic  grade  III). 

The  average  tissue  background  optical  absorption  coefficients  and  reduced  scattering  coefficients 
iu\  at  780  nm  and  830  run  were  measured  at  the  normal  contralateral  breast  as  p’*°=0.07  cm"',  = 

5.32cm'',  0.06  cm*'  and  p'f®=5.83  cm'',  respectively.  Optical  absorption  maps  of  lesion  region 

at  both  wavelengths  780nm  and  830  nm  were  reconstmcted  as  shown  in  Fig.  2(b)  and  (c).  The  first 
slice  is  the  spatial  x-y  image  of  0.5  cm  deep  from  the  skin  surface  and  the  last  slice  is  3.5  cm  deep 
toward  the  chest  wall.  The  spacing  between  the  slices  is  0.5  cm.  The  vertical  scale  is  the  absorption 
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coefficient  in  cm'’.  The  lesion  is  well  resolved  in  slice  #2.  The  total  hemoglobin  concentration  map  is 
shown  in  Fig  2.  (d)  and  the  vertical  scale  is  in  //M.  The  measured  maximum  value  at  lesion  area  is 
162.0  //M  with  the  average  value  of  121.1  // M  calculated  within  full  width  of  half  maximum  (FWHM) 
of  the  maximum  value.  The  average  background  total  hemoglobin  outside  of  the  FWHM  region  is  29.1 
H  M.  Since  our  hand-held  probe  can  be  easily  rotated  or  translated,  we  have  acquired  at  least  three  co¬ 
registered  ultrasound  and  NIR  data  sets  at  the  lesion  location  and  reconstructed  the  corresponding 
optical  absorption  maps  as  well  as  the  total  hemoglobin  concentration  distribution  under  the  co¬ 
registered  ultrasound  guidance.  The  average  maximum  total  hemoglobin  concentration  at  the  cancer 
region  is  149.91  p moles  (±39.23  pmoles)  and  the  average  background  total  hemoglobin  concentration 
is  29.03  p  moles  (±0.18  p  moles).  The  large  standard  deviation  obtained  from  slightly  different  probe 
positions  is  likely  due  to  the  probe  compression  at  different  positions  and  elastic  properties  of  the  blood 
vessels  in  the  cancer  region.  We  have  noticed  in  our  earlier  studies  that  one  of  the  two  early  stage 
invasive  cancers  showed  a  large  standard  deviation  on  total  hemoglobin  concentrations  at  the  cancer 
region  with  similar  scales  (19). 

Another  example  was  obtained  from  a  37-year  old  woman  who  had  a  known  infiltrating  lobular 
carcinoma  (nuclear  grade  II-III)  located  at  12  o’clock  position  of  the  right  breast  and  measuring  1  cm  in 
size  by  US  at  the  time  of  this  study  (Fig.  3(a)).  The  patient  had  family  history  of  breast  cancer.  Two 
new  lesions  located  at  7  o’clock  of  the  same  breast  and  10  o’clock  of  the  contralateral  breast  were  shown 
in  MRI.  The  new  lesions  were  identified  sono graphically  at  the  time  of  biopsy  measuring  6  and  7  mm 
in  diameter,  respectively.  Figure  4  (a)  is  the  US  image  of  the  7  o’clock  right  breast  lesion.  No  blood 
vessels  or  flow  were  identified  in  all  lesions  by  Doppler  US. 
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The  average  tissue  background  optical  absorption  coefficients  and  reduced  scattering  coefficients 
at  780  nm  and  830  nm  were  measured  at  contralateral  normal  breast  area  as  p™=0.03  cm'', 
^,780  _  g  08cm'',  pf°=  0.03  cm''  and  p'f°=7.76  cm'',  respectively.  The  absorption  maps  of  the  known 

cancer  at  both  wavelength  are  shown  in  Fig.3  (b)  and  (c)  and  the  computed  total  hemoglobin 
concentration  map  is  shown  in  Fig.3  (d).  The  lesion  is  resolved  in  slices  #3  and  #4  and  shown  110.4 
fj.  M  peak  and  77.8  /j.  M  average  total  hemoglobin  concentration,  respectively.  The  mean  maximum 
hemoglobin  concentration  of  three  NIR  scans  at  slightly  different  probe  locations  is  102.1  juM  and  the 
mean  within  FWHM  of  the  three  NIR  scans  is  68.7  //M  .  The  absorption  maps  of  the  new  lesion 
located  at  the  7  o’clock  position  of  the  same  breast  are  shown  in  Fig.  4(b)  and  (c)  and  the  total 
hemoglobin  concentration  is  shown  in  Fig.  4(d).  No  resolvable  lesion  was  found  in  absorption  maps 
and  the  maximum  total  hemoglobin  concentration  is  only  24.2  //  M.  The  biopsy  result  revealed  benign 
non-proliferative  fibroadipose  breast  tissue  with  focal  microcystic  alteration,  apocrine  metaplasia  and 
periductal  chronic  mastitis.  No  atypical  cells  were  found.  The  10  o’clock  suspicious  lesion  of  the 
contralateral  breast  had  similar  ultrasonic  and  optical  characteristics.  No  resolvable  lesion  was  found  in 
optical  absorption  maps  and  the  maximum  total  hemoglobin  concentration  is  only  27.5  jj.  M.  The  biopsy 
result  revealed  benign  proliferative  fibroadipose  breast  tissue  with  focal  aprocrine  papillary  epitheliosis. 
No  atypical  cells  were  found.  This  example  demonstrates  that  the  specificity  of  optical  contrast  is  much 
higher  than  that  of  Doppler  US  as  it  is  related  to  the  microvessel  development  of  tumor  angiogenesis 
process  in  contrast  to  the  larger  vessels  seen  by  Doppler  US. 

The  third  example  was  obtained  from  a  75-year  old  woman  who  had  a  suspicious  lesion  with  solid 
component  (pointed  by  the  arrow  in  Fig.  5  (a))  adjacent  to  a  cyst.  The  lesion  was  located  at  9  o’clock  of 
the  right  breast.  Doppler  US  has  shown  several  large  blood  vessels  located  at  the  peripheral  of  the 
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lesion.  Optical  absorption  maps  at  both  wavelengths  as  well  as  the  total  hemoglobin  concentration  map 
are  shown  in  Fig.  5(b),  (c)  and  (d),  respectively.  The  absorption  maps  of  the  lesion  have  showed  low 
light  absorption  and  low  hemoglobin  concentration.  The  absorption  distributions  are  quite  diffused  as 
compared  with  the  localized  distributions  of  malignant  cancer  cases.  The  measured  maximum  total 
hemoglobin  concentration  is  38.18  //M  and  the  average  is  26.31  juM.  The  biopsy  result  revealed  that 
the  lesion  was  intraductuctal  papilloma  with  no  evidence  of  atypical  cells  or  malignancy. 

The  last  example  was  obtained  from  a  47-year  old  woman  who  had  a  suspicious  tabulated  solid 
nodule  shown  in  Fig.  6(a).  No  blood  flow  was  observed  in  Doppler  US.  The  patient  had  family 
history  of  breast  cancer.  Optical  absorption  maps  as  well  as  the  total  hemoglobin  concentration 
map  are  shown  in  Fig.6  (b)-(c).  The  lesion  was  resolved  in  slice  2  and  showed  83.0  /u  M  maximum 
and  60.8  //M  average  total  hemoglobin  concentration.  The  biopsy  result  revealed  benign  breast 
tissue  with  minimally  proliferative  flbrocystic  changes  and  dense  fibrosis.  This  lesion  has  the 
highest  total  hemoglobin  concentration  among  all  the  benign  cases  reported  in  this  paper.  Since 
this  patent  had  family  history  of  breast  cancer,  6  months  follow  up  is  scheduled.  The  possibility  of 
core  biopsy  might  have  missed  the  suspicious  area  cannot  be  completely  ruled  out. 

Tables  I  and  II  list  results  of  5  malignant  cases  (5  malignant  cancers  and  3  benign  lesions)  and  10  benign 
cases  (13  benign  lesions)  studied  so  far.  Starting  from  column  2,  the  tables  provide  biopsy  results, 
measured  lesion  size  by  US,  combined  x-ray,  US  and/or  MRI  diagnosis,  Doppler  US  findings,  measured 
mean  maximum  value  of  total  hemoglobin  concentration  with  mean  taken  at  three  probe  positions, 
measured  mean  average  value  of  total  hemoglobin  concentration  with  mean  taken  at  three  probe 
positions,  and  measured  lesion  size  by  NIR  using  FWFM.  In  Table  I,  two  patients  (Ref  #52  and  Ref 
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#23)  have  multiple  benign  lesions  in  addition  to  malignant  ones  and  the  results  are  grouped  in  Table  I 
under  the  same  reference  number.  The  first  important  finding  is  that  the  invasive  cancers  have  shown 
much  higher  total  hemoglobin  concentrations  than  those  of  benign  cases.  Figure  7  provides  statistics  of 
these  two  groups.  The  mean  values  of  maximum  and  average  total  hemoglobin  concentration  of 
malignant  group  (5  lesions)  are  110.3  //M  (±21.2  //M )  and  76.7  /nM  (±  18.3 //M),  respectively,  and 
mean  values  of  benign  group  (16  lesions)  are  47.2  /^M  (±  18.5  //M  )  and  32.9  //M  (±  13.2 //M), 
respectively.  A  more  than  two-fold  greater  optical  contrast  is  obtained  from  malignant  group  and  it 
suggests  that  the  optical  differentiation  of  benign  and  malignant  breast  masses  is  feasible.  The  mean 
sizes  measured  by  US  for  these  two  groups  were  1.22  cm  and  0.7  cm,  respectively.  The  second 
important  finding  is  that  optical  tomography  provides  much  higher  specificity  than  Doppler  US. 
Although  Doppler  US  blood  flow  was  demonstrated  in  three  out  of  five  malignant  cases,  the  blood  flow 
was  also  present  in  many  benign  cases.  Since  the  diffused  light  probes  microvessel  density  rather  than 
relatively  large  blood  vessels  as  the  Doppler  US,  the  specificity  of  optical  contrast  is  expected  to  be 
much  higher  than  Doppler  US. 

Discussion 

With  the  increased  public  awareness  of  the  potential  benefits  of  early  detection  of  breast  cancer,  more 
women  are  now  practicing  self-examination  and  undergo  periodic  screening.  As  a  result,  more  lesions 
are  being  detected  by  palpation,  screening  mammography  and  screening  MRI  for  high  risk  or  dense 
breasts.  Consequently,  radiologists  and  surgeons  are  performing  an  even  increasing  number  of  breast 
biopsies,  with  the  vast  majority  of  these  revealing  benign  findings.  Therefore,  alternative,  less  invasive 
methods  of  distinguishing  malignant  from  benign  masses  are  needed  to  reduce  the  number  of 
unnecessary  biopsies.  US  is  frequently  used  as  an  adjunct  tool  to  x-ray  mammography  in  differentiating 
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simple  cysts  from  solid  lesions  and  also  plays  an  important  role  in  guiding  interventional  procedures 
such  as  needle  aspiration,  core-needle  biopsy,  and  prebiopsy  needle  localization  (4,  22,23).  Recently, 
screening  US  has  also  been  advocated  for  the  dense  breast  (24).  However,  the  specificity  of  Doppler  US 
and  gray  scale  imaging  in  differentiating  benign  from  malignant  lesions  is  low  as  many  solid  lesions  are 
diagnosed  as  suspicious.  Our  technique  of  using  optical  tomography  with  US  localization  has 
demonstrated  a  great  potential  for  non-invasively  distinguishing  malignant  and  benign  solid  masses  and 
therefore  a  potential  role  for  reducing  unnecessary  biopsies.  Earlier  results  obtained  from  two  invasive 
early  stage  cancers  and  17  benign  lesions  have  shown  that  malignant  cancers  of  1  cm  in  size  present  an 
average  of  1 19  //  M  maximum  total  hemoglobin  concentration  while  the  benign  group  has  an  average  of 
67  //  M.  A  nearly  two  hold  higher  contrast  has  been  obtained.  This  paper  further  demonstrates  the 
diagnostic  potential  of  this  technique  with  5  cancer  cases  and  a  group  of  16  biopsied  benign  solid 
lesions.  Again,  the  invasive  cancers  reveal  about  two-fold  greater  total  hemoglobin  concentration 
(mean  110.3  )  compared  with  benign  cases  (mean  47.3  //M).  Currently,  more  patients  are  being 

recruited  to  this  study  from  multiple  clinical  sites  and  the  specificity  of  the  optical  method  with  US 
localization  in  breast  cancer  diagnosis  will  be  reported  in  the  near  future.  Our  technique  may  also  be 
useful  for  monitoring  chemotherapy  responses  of  breast  cancer  treatment  and  assessing  treatment 
efficacy  (20). 

Due  to  intense  light  scattering,  optical  tomography  alone  has  not  been  widely  used  in  clinical  studies. 
The  data  in  the  published  literature  have  been  limited  to  case  reports.  Furthermore,  the  optical  systems 
and  measurement  configurations  vary  among  different  research  groups.  Pogue  et  al  reported  pilot 
results  of  one  invasive  ductal  carcinoma  and  one  benign  fibroadenoma  (13).  The  reported  maximum 
total  hemoglobin  concentrations  were  68  pmoles  for  the  cancer  case  and  55  pmoles  for  the  benign 
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case.  Since  the  system  they  used  can  only  acquire  NIR  data  from  a  ring  area  with  optical  source  and 
detectors  deployed  around  the  breast  (two-dimensional  imager).  Therefore,  incomplete  information 
could  lead  to  smaller  reconstructed  absorption  coefficients  and  total  hemoglobin  concentrations  than 
those  reported  here  using  our  three-dimensional  NIR  imager.  In  another  study  (25),  the  authors  used  the 
a  priori  lesion  depth  information  acquired  from  a  separate  US  image  and  obtained  an  average  of  67 
pmole  total  hemoglobin  concentration  from  a  ductual  carcinoma  in  situ.  This  value  is  close  to  the 
average  of  76.7  p, moles  obtained  from  the  malignant  group  reported  in  this  paper. 

In  principle,  the  distribution  of  oxygenation  saturation  can  be  estimated  as  S  =  oxyHb 
/(oxyHb+deoxyHb)  with  deoxyHb  and  oxyHb  distributions  calculated  from  absorption  maps  at  the  two 
wavelengths  of  780  nm  and  830  nm.  However,  since  background  tissue  mainly  consists  of  water  and 
lipid  and  these  two  chromophores  contribute  to  the  total  absorption  estimate  as  well,  we  could  not  obtain 
reasonable  background  oxygenation  saturation  and  compare  it  with  lesion  oxygenation  saturation. 
Currently,  we  are  in  the  process  of  adding  more  wavelengths  to  the  NIR  system.  Nevertheless,  initial 
experience  with  a  small  group  of  cancer  cases  suggests  that  the  deoxyHb  concentration  of  malignant 
lesions  may  not  be  as  significant  as  stressed  by  previous  investigators  but  the  total  hemoglobin 
concentration  (deoxyHb  +  oxyHb)  is  highly  specific  in  differentiating  small  invasive  malignant  cancers 
from  benign  lesions. 
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Figure  Captions 

Fig.  1.  (a)  Picture  of  the  eombined  probe  and  frequeney  domain  NIR  imager,  (b)  Sensor 

distribution  of  the  combined  probe.  Smaller  circles  in  (b)  are  optical  source  fibers  and  big  circles 
are  deteetor  fibers.  A  eommercial  ultrasound  probe  is  located  at  the  center  of  the  combined  probe 
and  the  optical  source  and  detector  fibers  are  distributed  at  the  periphery  of  the  ultrasound  probe. 

Fig.  2.  (DD41P2).  (a)  US  image  of  a  suspicious  lesion  located  at  3  o’clock  of  the  left  breast  of  a 
45-year  old  woman.  Doppler  US  revealed  large  blood  vessels  both  inside  and  at  peripheral  of  the 
lesion,  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  830  nm,  respectively.  The  first  slice 
is  the  spatial  x-y  image  of  0.5  cm  deep  from  the  skin  surface  and  the  last  slice  is  3.5  cm  toward  the 
chest  wall.  The  spacing  between  the  slices  is  0.5  cm.  The  lesion  is  well  resolved  in  slice  #2.  (d) 
The  total  hemoglobin  concentration  map  computed  from  absorption  maps  of  (b)  and  (c).  The 
vertical  scale  is  the  total  hemoglobin  concentration  in  //  M. 

Fig.  3.  (ID52P7).  (a)  US  image  of  an  infiltrating  lobular  carcinoma  located  at  12  o’clock  position 
of  the  right  breast  of  a  37-year  old  woman.  The  lesion  was  measured  1  cm  in  size  by  US  at  the 
time  of  this  study.  No  blood  vessels  or  flow  were  seen  by  Doppler  US.  (b)  and  (c)  are  optieal 
absorption  maps  at  780  nm  and  830  nm,  respectively,  (d)  The  total  hemoglobin  concentration  map 
computed  from  absorption  maps  of  (b)  and  (c).  The  lesion  is  resolved  in  slices  3  and  4. 

Fig.  4.  (ID52P17).  (a)  US  image  of  a  suspicious  lesion  located  at  7  o’elock  of  the  right  breast  of 
the  same  patient  with  an  infiltrating  lobular  carcinoma  shown  in  Fig.3.  No  blood  vessels  or  flow 
were  seen  by  Doppler  US.  (b)  and  (c)  are  optieal  absorption  maps  at  780  nm  and  830  nm, 
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respectively,  (d)  The  total  hemoglobin  concentration  map  computed  from  absorption  maps  of  (b) 
and  (c). 

Fig.  5.  (ID54P3_P17  REF),  (a)  US  image  of  a  suspicious  lesion  located  at  9  o’clock  of  the  right 
breast  of  a  75-year  old  woman.  Large  blood  vessels  were  seen  by  Doppler  US  at  the  peripheral  of 
the  lesion,  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  830  nm,  respectively,  (d)  The 
total  hemoglobin  concentration  map  computed  from  absorption  maps  of  (b)  and  (c). 

Fig.  6.  (ID58P1).  (a)  US  image  of  a  suspicious  lesion  located  at  7  o’clock  of  the  left  breast  of  a 
47-year  old  woman.  The  patient  had  family  history  of  breast  cancer.  No  blood  flow  was  seen  by  ■ 
Doppler  US.  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  830  nm,  respectively,  (d)  The 
total  hemoglobin  concentration  map  computed  from  absorption  maps  of  (b)  and  (c). 

Fig.  7.  Statistics  of  measured  total  hemoglobin  concentration  of  benign  and  malignant  groups. 

Table  Captions 

Table  I:  Biopsy,  US  and  NIR  Parameters  of  5  Cancer  Cases 
Table  II:  Biopsy,  US  and  NIR  Parameters  of  1 1  Benign  Cases 
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Fig.  2.  (ID41P2).  (a)  US  image  of  a  suspicious  lesion  located  at  3  o’clock  of  the  left  breast  of  a 
45-year  old  woman.  Doppler  US  revealed  large  blood  vessels  both  inside  and  at  peripheral  of 
the  lesion,  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  830  nm,  respectively.  The 
first  slice  is  the  spatial  x-y  image  of  0.5  cm  deep  from  the  skin  surface  and  the  last  slice  is  3.5  cm 
toward  the  chest  wall.  The  spacing  between  the  slices  is  0.5  cm.  The  lesion  is  well  resolved  in 
slice  #2.  (d)  The  total  hemoglobin  concentration  map  computed  from  absorption  maps  of  (b)  and 
(c).  The  vertical  scale  is  the  total  hemoglobin  concentration  in  jd  M. 
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(c)  (d) 

Fig.  4.  (ID52P17_P38  REF),  (a)  US  image  of  a  suspicious  lesion  located  at  7  o’clock  of  the 
right  breast  of  the  same  patient  with  an  infiltrating  lobular  carcinoma  shown  in  Fig.3.  No  blood 
vessels  or  flow  were  seen  by  Doppler  US.  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and 
830  nm,  respectively,  (d)  The  total  hemoglobin  concentration  map  computed  from  absorption 
maps  of  (b)  and  (c). 
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Fig.  3.  (ID52P7).  (a)  US  image  of  an  infiltrating  lobular  carcinoma  located  at  12  o’clock 

position  of  the  right  breast  of  a  37-year  old  woman.  The  lesion  was  measured  1  cm  in  size  by 
US  at  the  time  of  this  study.  No  blood  vessels  or  flow  were  seen  by  Doppler  US.  (b)  and  (c)  are 
optical  absorption  maps  at  780  nm  and  830  nm,  respectively,  (d)  The  total  hemoglobin 
concentration  map  computed  from  absorption  maps  of  (b)  and  (c).  The  lesion  is  resolved  in 
slices  3  and  4. 
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Fig.  5.  (ID54P3_P17  REF),  (a)  US  image  of  a  suspicious  lesion  located  at  9  o’clock  of  the  right 
breast  of  a  75-year  old  woman.  Large  blood  vessels  were  seen  by  Doppler  US  at  the  peripheral 
of  the  lesion,  (b)  and  (c)  are  optical  absorption  maps  at  780  nm  and  830  nm,  respectively,  (d) 
The  total  hemoglobin  concentration  map  computed  from  absorption  maps  of  (b)  and  (c). 
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Table  I:  Biopsy,  US  and  NIR  Parameters  of  5  Cancer  Cases  and  3  Benign  Lesions 


Ref 

BiopSy 

Size(US) 

US 

Doppler 

Total  Hb 

Total  Hb 

Size 

# 

Results 

cm 

Diagnosis 

US 

(max  /jM  ) 

(ave) 

(NIR) 

fjM 

cm 

#41 

intraductal  and 

' 

invasive 

O.SxO.S 

suspicious 

both 

147.9 

109.5 

2.67 

mammary  duct 

(±39.2) 

(±30.8) 

carcinoma 

#52 

#1:  infiltrating 
lobular 

0.91x0.92 

known  at 
biopsy 

.  no  flow 

102.1 

68.7 

3.91 

carcinoma 
#2:  non- 

(±8.1) 

24.2 
(±  12.8) 

(±8.6) 

17.3 

(±7.3) 

proliferative 
fibroadipose  breast 

0.62x0.33 

suspicious 

no  flow 

5.25 

tissue 

#3 :  proliferative 
fibroadipose  breast 
tissue 

0.66x0.33 

suspicious 

no  flow 

27.5 

21.1 

4.69 

#30 

invasive  mammary 

highly 

97.5 

67.7 

duct  carcinoma  & 
DCIS 

1. 1x1.1 

suspicious 

peripheral 

(±  12.5) 

(±9.4) 

4.62 

#46 

infiltrating  ductal 

1. 2x0.8 

suspicious 

no  flow 

100.0 

68.4 

3.71 

carcinoma  &  DCIS 

(±9.6) 

(±6.2) 

#23 

#1:  intraductal  and 

known  at 

Peripheral 

103.8 

69.2 

invasive  mammary 
duct  Carcinoma 

3. lx 

biopsy 

(±0.8) 

i 

(±0.9) 

48.5 

4.91 

#2:  non-proliferate 

0.72x0.43 

suspicious 

inside 

68.7 

4.58 

fibrocystic  changes 

(±2.34) 

(±4.9) 
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Table  II:  Biopsy,  US  and  NIR  Parameters  of  10  Benign  Cases  (13  lesions) 


Ref. 

# 

Biopsy 

Results 

Size(US) 

cm 

US 

Diagnosis 

Doppler 

US 

Total  Hb 
(max//M  ) 

Total  Hb 
(ave) 

/jM 

Size 

(NIR) 

cm 

#54 

Intraductal 

papilloma 

1.88x1.11 

suspicious 

peripheral 

38.3 

(±1.4) 

26.6 

(+1.1) 

5.24 

#45 

fibroadenoma 

1.92x1.02 

suspicious 

both 

28.7 

(±2.4) 

17.6 

(±2.2) 

8.0 

#28 

Apocrine  cyst 

0.6x0.35 

probably 

benign 

peripheral 

36.1 

(±4.8) 

24.6 

(±3.4) 

7.0 

#31  , 

Breast  tissue 
showing  nodular 
dense  fibrosis 
hyalinized 
fibroadenoma 

#l*;1.2x0.72 

#2*:0.53x0.58 

suspicious 

suspicious 

Peripheral 

no  flow 

37.8 

(±0.3) 

37.8 

(±0.3) 

26.2 

(±0.2) 

26.2 

(±0.2) 

4.89 

fibroadenoma 

#3:0.69x0.47 

suspicious 

no  flow 

70.9 

(±2.4) 

49.9 
(±  1.2) 

3.24 

#25 

Sclerotic 

fibroadenoma 

1.13x0.88 

suspicious 

no  flow 

48.7 
(±  2.9) 

32.2 

(±2.2) 

2.64 

#68 

Proliferative 

fibrocystic 

changes 

0.42x0.53 

suspicious 

peripheral 

25.8 
(±  1.8) 

17.9 
(±  2.5) 

7.14 

#61 

Fibrocystic 

changes 

0.47x0.48 

suspicious 

peripheral 

59.5 

(±7.2) 

40.7 

(±6.1) 

4.21 

#57 

fibroadenoma 

0.6  X  0.4 

0.6  X  0.4 

both  are 
suspicious 

no  flow 

43.0 

(±1.8) 

2.47 

2.91 

#58 

fibrocystic 

changes 

0.51x0.38 

suspicious 

no  flow 

56.9 
(±  10.7) 

2.81 

#21 

fibroadenoma 

0.61x0.38 

suspicious 

no  flow 

44.1 

(±0.4) 

2.16 

*these  two  lesions  were  located  at  2:30  and  3  o’clock  positions  and  were  functionally  un-resolvable. 
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A  novel  two-step  reconstruction  scheme  using  a  combined  near-infrared  and  ultrasound  technique  and  its  utility 
in  imaging  distributions  of  optical  absorption  and  hemoglobin  concentration  of  breast  lesions  are  demonstrated. 
In  the  first-step  image  reconstruction,  the  entire  tissue  volume  is  segmented  based  on  initial  coregistered  ultra¬ 
sound  measurements  into  lesion  and  background  regions.  Reconstruction  is  performed  by  use  of  a  finer  grid 
for  lesion  region  and  a  coarse  grid  for  the  background  tissue.  As  a  result,  the  total  number  . of  voxels  with  un¬ 
known  absorption  can  be  maintained  on  the  same  order  of  total  measurements,  and  the  matrix  with  unknown 
total  absorption  distribution  is  appropriately  scaled  for  inversion,  In  the  second  step,  image  reconstruction 
is  refined  by  optimization  of  lesion  parameters  measured  from  ultrasound  images.  It  is  shown  that  detailed 
distributions  of  wavelength-dependent  absorption  and  hemoglobin  concentration  of  breast  carcinoma  can  be 
obtained  with  the  new  reconstruction  scheme.  ©  2003  Optical  Society  of  America 
OCIS  codes:  170.0170,  170.3010,  170.5270,  170.7170,  170.3830. 


Tumor  blood  volume  and  microvascular  density  are 
parameters  that  are  anatomically  and  functionally 
associated  with  tumor  angiogenesis.  During  the  past 
decade,  modeling  of  light  propagation  in  the  near- 
infrared  (NIR)  region,  combined  with  advancements 
of  light  source  and  detectors,  has  improved  diffused 
light  measurements  and  made  possible  the  application 
of  tomographic  techniques  for  characterizing  and 
imaging  tumor  angiogenesis.^’^  However,  the  NIR 
technique  has  not  been  widely  used  in  clinics,  and 
the  fundamental  problem  of  intense  light  scattering 
remains.  As  a  result,  diffusive  light  probes  a  wide¬ 
spread  region  instead  of  providing  information  along 
a  straight  line,  and  tomographic  image  reconstruction 
is,  in  general,  underdetermined  and  ill-posed.  Zhu 
et  al.^  and  Chen  et  al^  demonstrated  a  combined 
imaging  technique,  using  a  priori  lesion  structure  in¬ 
formation  provided  by  coregistered  ultrasound  images 
to  assist  NIR  imaging  reconstruction  in  phantom  stud¬ 
ies.  As  a  result,  the  NIR  image  reconstruction  is 
well  defined  and  less  sensitive  to  noise.  In  this  Letter 
we  report  on  our  novel  two-step  image  reconstruction 
scheme  that  uses  the  combined  approach  and  demon¬ 
strate  its  utility  in  imaging  tumor  absorption  and 
hemoglobin  distributions.  It  is  shown  that  detailed 
heterogeneous  distributions  of  wavelength-dependent 
optical  absorption  and  hemoglobin  concentration  of  a 
breast  carcinoma  can  be  obtained.  To  the  best  of  our 
knowledge,  such  detailed  distributions  have  not  been 
reported  in  the  literature. 

A  picture  of  our  combined  hand-held  probe  used 
in  clinical  studies  is  shown  in  Fig.  1(a),  and  the 
probe  dimensions  and  optical  sensor  distributions 
are  shown  in  Fig.  Kb).  The  combined  probe  consists 
of  a  commercial  ultrasound  one-dimensional  array 
located  at  the  center  of  the  probe  and  optical  source 
and  detector  fibers  distributed  at  the  periphery 
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and  connected  to  the  NIR  imager.  The  NIR  imager 
consists  of  12  dual-wavelength  source  channels  and 
8  parallel  receiving  channels.^  In  the  transmission 
part,  12  pairs  of  dual-wavelength  (780-  and  830-nm) 
laser  diodes  are  amplitude  modulated  at  140  MHz. 
In  the  reception  part,  8  photomultiplier  tubes  de¬ 
tect  diffusely  reflected  light  from  the  tissue.  Both 
the  amplitude  and  phase  at  each  source -detector 
pair  are  obtained,  and  the  resulting  total  number 
of  measurements  is  12  X  8  X  2  —  192.  The  com¬ 
bined  probe  is  made  of  a  black  plastic  plate  10  cm 
in  diameter;  therefore,  a  semi-infinite  boundary 
condition  can  be  used  for  the  NIR  measurement 
geometry.  The  amplitude  and  phase  measured  from 
the  normal  side  of  the  breast  are  used  to  calculate  the 
background  absorption  and  reduced  scattering  co¬ 
efficient  In  our  two-step  image  reconstruction, 

we  first  segment  tissue  volume  into  two  regions,  L 
and  jB,  that  contain  a  lesion  as  measured  from  coreg¬ 
istered  ultrasound  images  and  background  tissue, 
respectively.  We  use  the  Born  approximation  to 


Fig.  1.  (a)  Hand-held  combined  probe  and  a  frequency- 
domain  NIR  imager,  (b)  Sensor  distribution  of  the  com¬ 
bined  probe  (diameter,  10  cm).  The  smaller  circles  are 
optical  source  fibers,  and  the  bigger  circles  are  detector 
fibers.  A  commercial  ultrasound  probe  is  located  at  the 
center  and  its  dimensions  are  5.6  cm  by  1  cm. 

©  2003  Optical  Society  of  America 


338  OPTICS  LETTERS  /  Vol.  28,  No,  5  /  March  1,  2003 


relate  the  scattered  field  Usc{rsi,rdi,o))  measured 
at  source -detector  pair  i  to  absorption  variations 
AjjLai^')  in  each  volume  element  of  two  regions  within 
the  sample,  where  rsi  and  rdi  are  the  source  and 
detector  positions,  respectively.  We  then  discretize 
the  lesion  volume  and  the  background  volume  with 
different  voxel  sizes  (a  finer  grid  for  lesion  volume  and 
a  coarse  grid  for  background).  The  scattered  field 
can  then  be  approximated  as 

Usc^rsiy  Tdi,  O))  ^ 

xTy 

Lj  J 

+  Y  Gir^k,rdi)Uinciruk.rsi)  f  A/Xa(rOd^r''^ ,  (1) 

Bk  * 

where  and  r^k  are  the  centers  of  voxels  j  and  k 
in  lesion  volume  L  and  background  volume  S,  respec¬ 
tively^,  and  Uinc{r\  rsi)  and  G(r^  rdi)  are  the  incident 
wave  and  the  Green  function  of  the  semi-infinite  ge¬ 
ometry,  respectively.  The  matrix  form  of  relation  (1) 
is  given  as 

[Usd]MXl  =  [Wl,  Wb]mXn[Mx„M5]^,  (2) 

where  ~  r ^siOJmxn^,  and  Wb  ~ 

[-yDGirykyrdi)Uincirvkyrsi)]Mxr^B  are  weight  matrices 
for  the  lesion  volume  and  the  background  volume,  re¬ 
spectively;  [Ml]  =  [  fi^  Afiair')d^r, A/iia{r')d^r^] 
and  [Mb]  =  are 

the  total  absorption  distributions  of  the  lesion  volume 
and  the  background  volume,  respectively. 

Instead  of  reconstructing  the  A/^a  distribution  di¬ 
rectly  as  the  standard  Born  approximation,  we  recon¬ 
struct  total  absorption  distribution  M  and  then  divide 
the  total  by  different  voxel  sizes  of  lesion  and  back¬ 
ground  tissue  to  obtain  A^ta  distribution.  By  choos¬ 
ing  a  finer  grid  for  the  lesion  and  a  coarse  grid  for  the 
background  tissue,  we  can  maintain  the  total  number 
of  voxels  with  unknown  absorption  on  the  same  scale 
of  the  total  measurements.  As  a  result,  the  inverse 
problem  is  less  underdetermined  and  ill-posed.  In  ad¬ 
dition,  since  the  lesion  absorption  coefficient  is  higher 
than  that  of  the  background  tissue,  in  general,  the  total 
absorption  of  the  lesion  over  a  smaller  voxel  is  on  the 
same  scale  of  total  absorption  of  the  background  over  a 
bigger  voxel,  and  therefore  the  matrix  [Ml,Mb]  is  ap¬ 
propriately  scaled  for  inversion.  The  reconstruction  is 
formulated  as  a  least-squares  problem.  The  unknown 
distribution  M  can  be  iteratively  calculated  with  the 
conjugate-gradient  search  method. 

The  lesion  location  and  volume  from  coregistered 
ultrasound  is  estimated  as  follows:  Since  the  com¬ 
mercial  one-dimensional  liltrasound  probe  that  we  use 
acquires  two-dimensional  ultrasound  images  in  the 
y-z  plane  {z  is  the  propagation  direction)  and  the  two- 
dimensional  NIR  probe  provides  three-dimensional 
images,  the  coregistration  is  limited  to  an  interception 
plane.  However,  if  we  approximate  a  lesion  as  an 
ellipsoid,  we  are  able  to  estimate  its  center  and  radii 
from  two  orthogonal  ultrasound  images  and  therefore 
obtain  the  lesion  volume.  Three  sources  of  error  may 


lead  to  inaccurate  estimation  of  the  lesion  center  and 
radii  and  therefore  to  cause  errors  in  reconstructed 
optical  properties.  First,  the  lesion  boundaries  may 
not  be  well  defined  in  ultrasound  images.  Second, 
two  separate  orthogonal  ultrasound  images  are  used 
to  estimate  the  radii  and  the  center,  and  these  pa¬ 
rameters  depend  on  the  ultrasound  probe  position  and 
compression  of  the  hand-held  probe.  Third,  the  target 
volumes  or  shapes  seen  by  different  modalities  may 
be  different  because  of  different  contrast  mechanisms. 
In  the  second  step,  we  refine  image  reconstruction 
by  perturbing  the  center  cq  and  then  the  radii  ro  and 
choosing  the  optimal  set  of  parameters  (copt,  ropt). 

Clinical  studies  were  performed  at  the  Health  Cen¬ 
ter  of  the  University  of  Connecticut,  and  the  human 
subject  protocol  was  approved  by  the  Health  Center 
IRB  committee.  Patients  with  palpable  and  nonpal- 
pable  masses  that  were  visible  on  clinical  ultrasound 
were  used  as  subjects.  These  subjects  were  scanned 
with  the  combined  probe,  and  ultrasound  images 
and  optical  measurements  were  acquired  at  multiple 
locations,  including  the  lesion  region  that  was  scanned 
at  two  orthogonal  positions,  and  a  normal  region  of 
the  contralateral  breast  scanned  at  two  orthogonal 
positions. 

An  example  is  given  in  this  Letter  to  demonstrate 
the  use  of  our  reconstruction  scheme.  Figure  2(a) 
shows  a  gray-scale  ultrasound  image  of  a  palpable 
lump  in  a  44-year-old  woman.  The  lesion  was  located 
at  the  6  to  8  o’clock  position  of  the  left  breast  at 
approximately  1,5-cm  depth.  Ultrasound  showed  an 
irregular  poorly  defined  hypoechoic  mass,  and  the 
lesion  was  considered  highly  suspicious  for  malig¬ 
nancy.  An  ultrasound  guided-core  needle  biopsy  was 
recommended.  Biopsy  results  (after  NIR  imaging) 
revealed  that  the  lesion  was  a  high-grade  in  situ 
ductal  carcinoma  with  necrosis. 

Multiple  optical  measurements  at  two  orthogonal 
positions  were  simultaneously  made  with  ultra¬ 
sound  images  at  the  lesion  location  as  well  as  at 
approximately  the  same  location  of  the  contralateral 
normal  breast.  The  fitted  average  tissue  background 
measured  on  the  normal  side  of  the  breast  at  both 
wavelengths  was  TT™  =  0.03  cm“\  =  0.05  cm"\ 
/if 80  =  9.22  cm“^,  and  =  7.58  cm^h  The 

perturbations  for  both  wavelengths  used  to  cal¬ 
culate  absorption  maps  were  normalized  as 
Use  {rsi,  rdi,  0))  ==  [Uiirsi,  rdi,  (o)  -  UN{rsi,  rdi,  o))] 

UB{rshrdi,a))/UN{rsiyrdi,(o),  where  UL{rsi,rdi,o))  and 
UN{rsi,rdi,o))  were  measurements  obtained  from  the 
lesion  region  and  the  contralateral  normal  region, 
respectively,  and  UB{rsi,rdi,o))  was  the  incident  field 
calculated  with  fitted  background  ~jl^  and  D  =  1/3 fXg'. 
This  procedure  cancels  unknown  system  gains  as¬ 
sociated  with  different  sources  and  detectors  as 
well  as  electronic  channels.  The  initial  estimates 
of  the  lesion  center  and  diameter  in  two  orthogonal 
ultrasound  images  were  (0,0.39  cm,  1.7  cm)  and 
3.44  cm  X  4.38  cm  X  1.76  cm,  respectively.  A  finer 
grid  of  0.5  cm  X  0.5  cm  X  0.5  cm  and  a  coarse  grid  of 
1.5  cm  X  1.5  cm  X  1.0  cm  were  chosen  for  the  lesion 
and  background  tissue,  respectively.  The  total  recon¬ 
struction  volume  was  chosen  to  be  9  cm  X  9  cm  X  4  cm, 
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Fig.  2.*  (a)  Gray-scale  ultrasound  image  of  a  palpable  lump  of  a  44-year-old  woman.  Ultrasound  showed  an  irregular 
poorly  defined  hypoechoic  mass,  and  the  lesion  was  considered  as  highly  suspicious  for  malignancy.  Reconstructed  optical 
absorption  maps  at  (b)  780  nm  and  (c)  830  nm.  The  vertical  color  bars  are  the  absorption  coefficient  [cm"^].  (d)  Total 
hemoglobin  concentration  map.  The  vertical  color  bars  are  ^tmol.  The  NIR  data  were  simultaneously  acquired  with 
the  ultrasound  image  shown  in  (a).  Each  image  consists  of  seven  slices  obtained  in  0.5-cm  spacing  from  0.5  to  3.5  cm  in 
depth.  The  vertical  and  horizontal  axes  correspond  to  x  and  y  dimensions  of  9  cm  by  9  cm. 


and  the  total  number  of  voxels  with  unknown  optical 
absorption  was  190,  which  was  of  the  same  order  as  the 
192  total  measurements.  Image  reconstruction  was 
performed  with  the  NIR  data  simultaneously  acquired 
with  the  ultrasound  image  shown  in  Fig.  2(a).  The 
second-step  refined  reconstruction  revealed  optimal  le¬ 
sion  centers  at  approximately  (~1.1  cm,  0.3  cm,  1.7  cm) 
for  780  nm  and  (--0.9  cm,  “0.7  cm,  1.7  cm)  for  830  nm 
and  optimal  diameters  of  4.28  cm  X  5.18  cm  X  1.96  cm. 
The  detailed  absorption  maps  with  high  absorption 
nonuniformly  distributed  around  the  lesion  bounda¬ 
ries  at  both  wavelengths  are  shown  in  Figs.  2(b)  and 
2(c).  By  assuming  that  the  major  chromophores  are 
oxygenated  (oxyHb)  and  deoxygenated  (deoxyHb)  he¬ 
moglobin  molecules  in  the  wavelength  range  studied, 
we  can  estimate  the  distribution  of  total  hemoglobin 
concentration  as  shown  in  Fig.  2(d).  The  measured 
average  cancer  and  background  total  hemoglobin 
concentrations  were  55.8  ^tmol  and  20.7  /^mol, 
respectively. 

It  is  interesting  to  note  that  the  absorption  dis¬ 
tributions  at  both  wavelengths  as  well  as  the  total 
hemoglobin  concentration  were  distributed  heteroge¬ 
neously  at  the  cancer  periphery.  To  the  best  of  our 
knowledge,  such  fine  distributions  have  not  been  ob¬ 
tained  by  use  of  NIR-only  reconstruction  techniques. 
However,  this  finding  agrees  with  the  published 
literature  showing  that  breast  cancers  have  higher 


blood  volumes  than  nonmalignant  tissue  because  of 
angiogenesis,  especially  at  the  cancer  periphery.  In 
addition,  the  carcinoma  reported  here  had  a  necrotic 
core,  which  could  lead  to  the  low  absorption  observed 
at  both  wavelengths  in  the  center  region. 
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Abstract 


In  this  paper,  we  report  high-resolution  identification  of  angiogenesis  distributions  of  advanced  breast 
cancers  imaged  by  a  novel  hybrid  optical  and  ultrasound  imaging  technique.  The  hybrid  technique 
deploys  near-infrared  (NIR)  optical  imaging  sensors  and  a  commercial  ultrasound  transducer  on  a  hand¬ 
held  probe.  The  ultrasound  transducer  is  used  for  lesion  localization  and  the  optical  sensors  are  used  for 
imaging  tumor  angiogenesis.  Preliminary  results  of  four  advanced  cancer  cases  have  shown  correlation 
between  cancer  angiogenesis  distributions  and  histological  microvessel  density  counts. 
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Introduction 

Noninvasive  in  vivo  breast  imaging  at  the  molecular  level  presents  a  unique  challenge  .  Optical 
tomography  using  near  infrared  diffused  light  has  a  great  potential  for  monitoring  tumor  angiogenesis 
development,  a  key  factor  required  for  tumor  growth  and  metastasis.  This  method  also  provides  insight 
into  tumor  metabolism  and  tumor  hypoxia,  an  important  indicator  of  tumor  response  to  various  forms  of 
therapy  (1-5).  These  functional  parameters  are  highly  felavant  for  assessing  tumor  response  to 
treatment  as  well  as  for  diagnosis.  However,  optical  tomography  when  used  alone  has  been  limited  to 
laboratory  tests  and  feasibility  studies  in  the  past.  The  fundamental  problem  remains  the  intense  light 
scattering  in  tissue,  which  makes  the  localization  of  lesion  difficult.  A  novel  ultrasound  guided  optical 
imaging  technique  has  shown  promising  results  of  overcoming  the  diffusive  light  localization  problem 
in  tissue  (6-8).  The  technique  is  implemented  by  simultaneously  deploying  optical  sensors  and  a 
commercial  ultrasound  transducer  on  a  hand-held  probe,  and  utilizing  co-registered  lesion  structure 
information  provided  by  ultrasound  to  improve  the  inverse  optical  imaging  reconstruction.  Preliminary 
results  have  shown  that  early  stage  invasive  breast  cancers  present  two-fold  high  total  hemoglobin 
concentration  than  a  group  of  fibroadenomas  and  benign  lesions  (9). 

Tumor  angiogenesis  distribution  of  advanced  cancers  is  highly  heterogeneous  (7).  This  is  related  to  the 
complex  angiogenesis  process  involving  both  the  incorporation  of  existing  host  blood  vessels  into  tumor 
and  the  creation  of  tumor  microvessels.  This  process  is  modeled  by  angiogenesis  factors  (10).  Blood 
flow  through  these  tumor  vessels  is  heterogeneous.  Some  areas  have  high  flow,  others  have  slower  flow 
and  develop  necrosis  (10).  Tumors  with  relatively  poor  blood  perfusion  may  not  receive  adequate 
delivery  of  systemic  therapy.  This  lack  of  perfusion  may  be  a  factor  in  poor  response  to  intravenous 
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chemotherapy  treatment  in  some  patients  (1 1).  Although  pilot  studies  of  imaging  tumor  angiogenesis 
have  been  conducted  with  optical  tomography,  no  one  has  obtained  high-resolution  heterogeneous 
angiogenesis  distributions  in  advanced  cancers.  In  this  paper,  we  report  what  we  believe  the  first 
heterogeneous  angiogenesis  distributions  of  large  cancers  imaged  by  optical  tomography.  We  also  show 
that  the  angiogenesis  distributions  correlate,  to  a  large  extent,  to  histological  microvessel  density  counts. 

2.  Materials  and  Methods 

A  hand-held  hybrid  probe  consists  of  a  commercial  ultrasound  transducer  located  at  the  middle  and 
near-infrared  source-detector  fibers  distributed  at  the  periphery.  The  picture  of  the  probe  and  its 
dimensions  are  shown  in  Fig.l.  The  technical  aspects  of  the  NIR  imager  have  been  described  in  detail 
previously  (8).  Briefly,  the  imager  consists  of  12  pairs  of  dual  wavelength  (780nm  and  830nm)  laser 
diodes,  which  are  used  as  light  sources,  and  their  outputs  are  coupled  to  the  probe  through  optical  fibers. 
On  the  receiving  side,  8  photomultiplier  tubes  (PMTs)  were  used  to  detect  diffusely  scattered  light  from 
the  tissue  and  8  optical  fibers  were  used  to  couple  detected  light  to  the  PMTs.  The  laser  diodes’  outputs 
were  amplitude  modulated  at  140  MHz  and  the  detector  outputs  were  demodulated  to  20KHz.  The 
demodulated  signals  were  further  amplified  and  band  pass  filtered  at  20  KHz.  A  reference  signal  of  20 
KHz  was  also  generated  by  directly  mixing  the  detected  RF  signals  with  the  RF  signal  generated  from 
the  oscillator.  The  reference  signal  was  necessary  for  retrieving  phase  shifts.  Eight  detection  signals 
and  one  reference  were  sampled  and  acquired  into  a  PC  simultaneously.  The  entire  data  acquisition  took 
about  3  to  4  seconds,  which  was  fast  enough  for  acquiring  data  from  patients. 

The  details  of  our  dual-mesh  optical  imaging  reconstruction  algorithm  have  been  described  in  Ref  (9). 
Briefly,  the  NIR  reconstruction  takes  advantages  of  ultrasound  localization  of  lesions  and  segments  the 
imaging  volume  into  finer  grid  in  lesion  region  and  coarser  grid  in  non-lesion  regions.  A  modified  Bom 
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approximation  is  used  to  relate  the  scattered  field  measured  at  the  source-detector  pairs  to  absorption 
variations  in  each  volume  element  of  two  regions  within  the  sample.  With  this  dual-mesh  scheme,  the 
inverse  optical  reconstruction  is  well-conditioned  and  converges  in  few  iterations. 

Clinical  studies  were  performed  at  the  University  of  Connecticut  Health  Center  (UCHC).  The  UCHC 
IRB  committee  approved  the  human  subject  protocol.  Patients  with  palpable  and  non-palpable  masses 
that  were  visible  on  clinical  ultrasoimd  and  who  were  scheduled  for  biopsy  were  enrolled  as  research 
subjects.  Ultrasound  images  and  optical  measurements  were  acquired  simultaneously  at  multiple 
locations  including  the  lesion  region  scanned  at  two  orthogonal  positions,  a  normal  region  of  the  same 
breast  if  the  breast  was  large,  and  a  normal  symmetric  region  of  the  contralateral  breast  scanned  at  two 
orthogonal  positions.  The  optical  data  acquired  at  normal  region  were  used  as  reference  for  calculating 
the  scattered  field  caused  by  lesions. 

Of  patients  studied  so  far,  four  cases  presented  with  advanced  breast  cancers;  tumor  sizes  ranging  from 
T2  to  T4.  One  patient  received  chemotherapy  for  three  months  before  surgery,  and  the  other  three  have 
received  definitive  surgery.  To  correlate  the  imaged  angiogenesis  distribution  with  histology 
microvessels  density,  we  have  performed  microvessel  density  counts.  Samples  obtained  at  surgical 
biopsy  or  breast-conservation  surgery  were  used  for  counting.  For  each  sample,  sections  3  to  5 
micrometer  thick  were  stained  on  an  immunohistochemistry  slide  staining  system  (DAKO  autostainer) 
with  factor  8/86  mouse  monoclonal  antibody  (anti-human  von  Willebrand  factor,  DAKO  Corp, 
Carpinteria,  Calif)  at  1:100  dilution  digested  by  proteinase  K  for  3  minutes,  by  labeled  polymers 
(DAKO  Envision  plus)  immunoperoxidase  method.  The  microvessel  density  counts  were  performed  in 
ten  consecutive  fields  with  the  use  of  an  ocular  grid  at  X200  magnification.  The  first  field  chosen  was  a 
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hotspot  (area  of  maximum  vascular  density  either  within  the  infiltrating  tumor  mass  or  at  the  tumor- 
stromal  interface).  Since  the  cancers  were  large,  two  to  three  separate  sample  blocks  were  selected  for 
microvessel  counts  as  per  the  specimen  orientation  in  the  surgical  pathology  report. 

3.  Results 

During  the  diagnostic  imaging  study,  we  scanned  a  patient  who  was  undergoing  chemotherapy.  This 
44-year-old  woman  had  a  large  4  cm  x  4  cm  x  1.5  cm  palpable  mass  (Fig.2  (a))  located  at  the  6  to  8 
o’clock  position  of  the  left  breast  that  was  considered  as  highly  suspicious  for  malignancy.  The  lesion 
center  was  approximately  1.5  cm  in  depth  relative  to  the  skin.  An  ultrasound  guided  needle  biopsy 
revealed  that  the  lesion  was  a  high-grade  invasive  carcinoma  with  necrosis.  Optical  absorption  maps  of 
both  wavelengths  are  shown  in  Fig.  2(b)  and  (c)  and  the  distributions  are  highly  heterogeneous  with  high 
absorption  at  the  cancer  periphery.  Slice  1  is  the  spatial  x-y  image  of  9  em  x  9  em  obtained  at  0.5  cm 
deep  from  the  skin  surface.  Slice  7  is  3.5  cm  deep  toward  the  ehest  wall  and  the  spacing  between  slices 
is  0.5  cm.  The  total  hemoglobin  eoncentration  map  is  shown  in  Fig.2  (d)  and  the  maximum  and  average 
hemoglobin  concentrations  are  92.1  ?  moles  and  26.2  '!  moles.  The  measured  maximum  absorption 
coefficients  at  780  nm  and  830  nm  are  0. 17  cm*'  and  0.22  cm"',  respectively.  Since  this  cancer  was  too 
large  for  breast  conserving  surgery,  the  patient  was  treated  with  chemotherapy  in  the  neo-adjuvant 
setting  for  three  months.  At  the  time  the  patient  eompleted  the  chemotherapy,  we  imaged  her  again  with 
the  combined  US  and  NIR  probe.  Figure  3  (a)  is  the  ultrasound  image  of  the  eancer  three  month  later. 
The  eancer  contrast  is  poor  and  cancer  boundaries  are  completely  unclear,  probably  due  to  treatment. 
Figure  3  (b)  and  (c)  are  optical  absorption  maps  at  both  wavelengths  and  (d)  is  the  total  hemoglobin 
distribution.  The  maximum  and  average  hemoglobin  concentrations  of  the  lesion  are  79.8  ?  moles  and 
24.9  ?  moles,  respectively.  The  measured  maximum  absorption  coefficients  at  780  nm  and  830  nm  are 
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0.15  cm'  and  0.19  cm'',  respectively.  Compared  with  the  images  acquired  before  treatment,  the  spatial 
extension  of  the  angiogenesis  pattern  is  much  smaller  and  more  confined  to  the  core  area.  The 
maximum  total  hemoglobin  concentration  is  reduced  by  about  10  ?  moles  and  the  average  is  about  the 
same  as  before.  This  example  demonstrates  the  feasibility  of  monitoring  the  treatment  using  the 
combined  technique. 

To  correlate  the  NIR  images  with  microvessel  densities,  w.e  selected  three  block  samples  obtained  at 
breast-conservation  surgery  marked  as  anterior,  lateral  and  posterior  (see  Fig.  3  insertion)  for 
microvessel  counting.  The  total  number  of  microvessels  were  196  (lateral),  114  (anterior)  and  48 
(posterior  and  inferior)  per  10  consecutive  fields  at  200  magnification,  respectively  (see  table  I:  ID#  19). 
The  high  coimts  obtained  at  anterior  and  lateral  correlate  with  the  high  optical  absorption  and  total 
hemoglobin  concentration  shown  in  slice  3  of  Fig.3  (b)-(d).  A  representative  section  demonstrating 
high  microvessel  density  is  shown  in  Fig.3  (e).  The  low  counts  obtained  at  posterior  and  inferior  block 
sample  correlate  with  the  low  optical  absorption  as  well  as  the  low  hemoglobin  concentration  seen  in 
deeper  slices  4  and  5  of  Fig.  3(a)-(c). 

The  second  imaging  example  was  obtained  from  a  47-year-old  woman  who  had  a  3  cm  x  3  cm  x  2  cm 
dominant  mass  at  the  2  o’clock  position  in  her  left  breast.  The  lesion  center  was  about  2.3  cm  in  depth 
relative  to  the  skin.  Ultrasound  showed  hypoechoic  mass  with  irregular  margins  (see  Fig.  4(a))  and  the 
lesion  was  considered  as  highly  suspicious  for  malignancy.  Figure  4  (b)  and  (c)  are  optical  absorption 
maps  obtained  at  780  nm  and  830  nm  and  the  vertical  scale  of  (b)  and  (c)  is  half  of  that  used  in  Fig.3  (b) 
and  (c).  Figure  4(d)  is  the  total  hemoglobin  concehtration  distribution  and  the  vertical  scale  is  the  same 
as  that  used  in  Fig.  3(d).  The  light  absorption  at  both  wavelengths  is  much  lower  than  the  previous 
example  but  the  distributions  are  highly  heterogeneous.  The  measured  maximum  absorption 
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coefficients  at  780  nm  and  830  nm  are  0.08  cni'  and  0.10  cm'',  respectively.  The  measured  maximum 
total  hemoglobin  concentration  of  the  tumor  and  average  of  the  lesion  are  40.6  ?  moles  and  17.2 
?  moles,  respectively.  Surgical  pathology  report  revealed  that  the  mass  was  an  infiltrating  carcinoma 
(histological  grade  II,  nuclear  grade  II)  with  low  mitotic  activity.  The  total  counts  of  microvessels 
obtained  from  anterior  and  posterior  core  biopsy  samples  were  61  and  40  per  10  consecutive  fields  at 
200  magnification,  respectively  (see  Table  I:  ID#5).  The  total  counts  measured  from  anterior  and 
posterior  tumor  samples  were  52  and  29,  respectively.  These  low  counts  correlate  well  with  the  low 
optical  absorption  shown  in  Fig.  4  (b)-(c)  and  indicate  that  the  tumor  was  poorly  perfused.  A 
representative  section  demonstrating  low  microvessel  density  is  shown  in  Fig.4  (e). 


Since  our  hand-held  probe  can  be  easily  rotated  or  translated,  we  have  acquired  at  least  three  co¬ 
registered  ultrasound  and  NIR  data  sets  at  the  lesion  location  for  each  patient  and  reconstructed 
corresponding  optical  absorption  maps  as  well  as  the  total  hemoglobin  concentration  distribution  under 
the  co-registered  ultrasound  guidance.  The  measurements  shown  in  Table  II  are  average  values  taken 
over  three  images.  The  maximum  absorption  coefficient,  the  average  absorption  coefficient  within  the 
hot  area  defined  as  within  6  dB  of  the  maximum  value,  the  ratio  of  average  absorption  coefficient  inside 
the  lesion  over  the  hot  area,  the  maximum  total  hemoglobin  concentration,  the  average  total  hemoglobin 
concentration  inside  the  lesion,  are  all  given  in  Table  II.  The  ratio  of  average  absorption  coefficient 
inside  the  lesion  over  the  hot  area  partially  reflects  the  angiogenesis  heterogeneity. 

The  third  patient  (Table  I:  ID#23)  was  a  33-year-old  pregnant  woman  who  had  palpable  left  breast  lump 
at  12’clock  position  measuring  3  cm  x  3cm  x  1.5  cm.  The  ultrasound  image  showed  that  the  lesion  had 
discrete  nodularity  (see  Fig.  7(a)).  The  center  of  the  lump  was  about  2.5  cm  deep  from  the  skin.  An 
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ultrasound  guided  core  biopsy  was  obtained  and  revealed  that  the  lesion  was  both  in  situ  and  invasive 
ductal  carcinoma  as  seen  on  multiple  cores  (histologic  grade  III,  nuclear  grade  III).  Optical  absorption 
maps  as  well  as  total  hemoglobin  concentration  map  were  obtained  8  days  after  core  biopsy  at  the  time 
of  the  patient  visit  (see  Fig.5  (b)-(c)).  The  measured  average  maximum  absorption  coefficients  at  780 
nm  and  830  nm  are  0.15  cm'’  and  0.33  cm"’,  respectively  (Table  II:  ID#23).  The  maximum  total 
hemoglobin  concentration  of  the  tumor  is  1 13.0  ?  moles  and  the  average  is  24.5  ?  moles.  As  seen  from 
Table  II,  the  standard  deviation  of  maximum  total  hemoglobin  concentration  measured  from  different 
probe  positions  is  21  ?  moles,  which  is  much  larger  than  those  obtained  from  other  cases.  Since  the 
optical  imaging  was  performed  eight  days  after  core  biopsy  and  possible  hemoglobin  changes  due  to  a 
host  response  at  biopsy  sites  could  not  be  ruled  out. 

Histological  microvessel  counts  of  three  sample  blocks  were  60  (anterior/lateral),  8 8 (anterior/medial), 
and  152  (posterior  and  inferior)  per  10  consecutive  fields  at  200  magnification,  respectively.  The  larger 
variation  in  total  counts  at  different  tumor  locations  is  partially  related  to  the  inherent  heterogeneity  of 
breast  tumors  and  the  resulting  misdistribution  of  angiogenesis  in  the  viable  and  schirrous  regions. 
Nevertheless,  the  relatively  higher  counts  obtained  at  anterior/lateral,  anterior/medial  correlate  to  some 
extent  to  the  high  optical  absorption  and  total  hemoglobin  distribution  seen  in  slice  4  of  Fig.  5(b)-(d). 
The  high  counts  obtained  in  the  posterior  and  inferior  sample  do  not  correlate  with  the  low  optical 
absorption  distribution  seen  by  optical  imaging  in  deeper  slices.  The  reason  for  this  discrepancy  is 
explained  in  the  following  section. 

The  last  example  was  obtained  from  a  5  3 -year  old  woman  who  had  a  palpable  mass,  but  a  normal  - 
mammogram.  An  ultrasound  obtained  revealed  an  irregularly  shaped  lesion  of  2  cm  x  2  cm  x  1.3  cm 
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and  ultrasound  guided  surgical  biopsy  confirmed  an  invasive  and  in  situ  ductal  carcinoma  (histological 
grade  II,  nuclear  grade  II).  The  lesion  center  was  about  2.5  cm  in  depth  relative  to  the  skin.  Optical 
absorption  maps  as  well  as  the  total  hemoglobin  concentration  distribution  were  obtained.  The 
measured  maximum  absorption  coefficients  at  780  nm  and  830  nm  are  0.20  cm'  and  0.29  cm*  (Table 
II:  ID#21),  respectively.  The  maximum  total  hemoglobin  concentration  of  the  tumor  is  114.8  ?  moles 
and  the  average  is  20.9  ?  moles.  The  total  number  of  microvessels  were  83  (anterior),  121  (posterior 
and  lateral),  124  (posterior)  per  10  consecutive  fields  at  200  magnification,  respectively.  The  higher 
anterior  and  lateral  counts  correlate  with  the  high  optical  absorption  and  high  total  hemoglobin 
concentration.  Similar  to  the  third  example,  the  higher  posterior  counts  do  not  correlate  with  the  low 
light  absorption  and  low  hemoglobin  concentration  seen  in  deeper  slices. 

4.  Discussion  and  Summary: 

Reflection  geometry  is  used  for  optical  tomography  because  optical  source  and  detector  fibers  can  be 
easily  integrated  with  a  hand-held  ultrasound  probe.  For  deeply  located,  highly  absorbing  tumors,  such 
as  the  third  and  fourth  examples,  the  diffusely  reflected  photon  density  waves  from  the  bottom  of  the 
tumor  are  weak  when  they  reach  the  detectors.  Therefore,  the  perturbations  from  the  deeper  part  of  the 
tumor  are  much  weaker  compared  with  perturbations  from  the  top  part  of  the  tumor.  As  a  result,  the 
reconstructed  images  show  higher  absorption  at  top  and  lower  absorption  at  the  bottom.  For  the  poorly 
perfused  case  (Fig.4),  the  light  absorption  of  the  cancer  was  not  high  and  the  lesion  was  imaged  more 
uniformly  from  top  to  bottom.  This  depth  dependent  distribution  imaged  by  dffused  wave  may  be 
minimized  by  increasing  the  detection  sensitivity  and  appropriately  scaling  the  weight  matrix  for 
imaging  reconstruction.  Currently,  we  are  pursuing  investigations  along  this  direction. 
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For  the  large  cancers  presented  here,  ultrasound  can  accurately  diagnose  the  malignancy.  The 
value  of  the  optical  tomography  is  to  estimate  tumor  perfusion  and  evaluate  treatment  efficacy  by 
providing  angiogenesis  maps  as  well  as  tumor  hypoxia  during  the  chemotherapy.  This  type  of 
information  could  prove  invaluable  in  determining  whether  or  not  the  tumors  are  chemosensitive. 
Tumor  hypoxia  which  can  also  be  measured  is  related  to  tumor  growth  rate,  drug  response  etc  (12). 

In  principle,  the  distribution  of  oxygenation  saturation  can  be  estimated  from  absorption  maps  at 
the  two  wavelengths  of  780  nm  and  830  nm.  However,  since  background  tissue  mainly  consists  of 
water  and  lipid  and  these  two  chromophores  contribute  to  the  total  absorption  estimate  as  well,  we 
.could  not  obtain  reasonable  background  oxygenation  saturation  and  compare  it  with  lesion 
oxygenation  saturation.  One  or  two  more  wavelengths  will  allow  us  to  accurately  estimate  the 
background  oxygenation  saturation  and  to  compare  it  with  lesion  oxygenation  saturation. 

In  this  paper,  we  have  reported  initial  results  of  heterogeneous  angiogenesis  distributions  of  four  average 
size  cancer  cases  imaged  by  optical  tomography  with  ultrasound  localization.  The  angiogenesis  maps 
correlate,  to  a  large  extent,  to  the  histological  microvessel  density  counts.  These  initial  results  are 
exciting  and  suggest  that  imaging  with  optical  tomography  at  molecular  level  is  feasible. 
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Figure  Captions: 

Figure  1.  (a)  The  hand-held  probe  with  ultrasound  transducer  located  in  the  middle  and  optical 
sources  and  detectors  distributed  at  the  periphery,  (b)  Optical  source  and  detector  distribution. 

Figure  2.  (a)  Ultrasound  image  of  a  44-year-old  woman  who  had  a  large  4  cm  by  4cm  palpable  left 
breast  mass,  (b)  Optical  absorption  map  at  780  nm.  (c)  Optical  absorption  map  at  830  nm.  (d) 
Total  hemoglobin  concentration  of  the  lesion. 


Figure  3  (ID#19P10).  (a)  Ultrasound  image  of  the  cancer  treated  with  chemotherapy  for  three 
months,  (b)  Optical  absorption  map  at  780  nm.  (c)  Optical  absorption  map  at  830  nm.  (d)  Total 
hemoglobin  concentration  of  the  lesion,  (e)  Representative  section  of  the  infiltrating  ductal 
carcinoma  of  breast  (ID#  19)  demonstrating  high  microvessel  density  (original  magnification  X200; 
anti-human  von  Willebrand  factor,  DAKO  EnVision  plus). 

Figure  4(ID#5P4).  (a)  Ultrasound  image  of  a  47-year-old  woman  with  a  highly  suspicious  lesion  located 
at  2’oclock  position  of  her  left  breast,  (b)  Optical  absorption  map  at  780  nm.  (c)  Optical  absorption  map  at 
830  nm.  Note  the  scales  for  Fig.  4(b)  and  (c)  is  half  of  that  used  in  Fig.  3(b)  and(c).  (d)  Total  hemoglobin 
concentration  of  the  lesion.  The  scale  is  the  same  as  the  Fig.3  (d).  (e)  Representative  section  of 
infiltrating  ductal  carcinoma  of  breast  (ID#5,  Fig.7)  demonstrating  low  microvessel  density 
(original  magnification  X200;  anti -human  von  Willebrand  factor,  DAKO  EnVision  plus). 

Figure  5  (ID#23P7).  (a)  Ultrasound  image  of  a  33-year-old  woman  with  a  highly  suspicious  lesion 
at  12’clock  position  measuring  3  cm  x  3cm  x  1.5  cm,  and  the  ultrasound  showed  that  the  lesion  had 
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discrete  nodularity.  Ultrasound  image  cannot  clearly  identify  the  lesion  sue  the  acoustic 
heterogeneous  of  tissue  at  pregnancy,  (b)  Optical  absorption  map  at  780  nm.  (c)  Optical 
absorption  map  at  830  nm.  (d)  Total  hemoglobin  concentration  of  the  lesion. 
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Table  I.  Histological  microvessel  density  counts 


ID# 

core  bx  # 

block  # 

location 

mvd/10200x 

block# 

mvd/10  200x 

NA 

NA 

NA 

S02-7823 

B2 

m^mm 

196 

5SS3H 

NA 

NA 

NA 

S02-7823 

B6 

114 

NONE 

NA 

NA 

S02-7823 

B8 

48 

SOl-5578 

61 

S02-392 

mm 

ANT 

52 

SOl-5578 

mM 

InB 

SQHHHI 

S02-392 

A2 

POST 

29 

ID#23 

NONE 

NA 

NA 

NA 

S02-4024 

A3 

152 

NONE 

NA 

NA 

NA 

S02-4024 

A4 

NA 

NA 

NA 

S02-4024 

mm 

ANT/MED 

88 

mUSM 

mm 

NA 

mmm 

S02-4299 

A3 

121 

5SS3Hi 

NA 

NA 

A2 

POST 

124 

NA 

NA 

A1 

ANT 

83 

Table  II.  Optical  measurements 


ID# 

max  ?a  cm"' 
(780  nm) 

ave.6dB 
?.  cm'' 

(780  nm) 

aveL/ave.6d 

B. 

(780  nm) 

max  ?a  cm'' 
(830  nm) 

ave.6dB 
?a  cm''  (830 
nm) 

aveL 
/ave.6dB 
(830  nm) 

max  total  Hb 
?  moles 

ave.  Hb 
?  moles 

ID#  19 

0.16(0.01) 

0.11(0.01) 

33.7% 

0.18(0.004) 

0.13(0.006) 

50.2% 

80.6(2.4) 

24.1(1.1) 

56.7% 

0.06  (0.01) 

67.4% 

17.2(1.1) 

ID#23 

SQSEIEHi 

mmm 

iiWKilAM 

SSEIESESH 
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We  propose  a  novel  method  for  measuring  time-dependent  optical  quantities.  A  train  of  excitation  pulses 
modulated  by  a  pseudorandom  bit  sequence  is  used  as  the  light  source,  and  a  cross-correlation  scheme  is 
used  to  retrieve  the  impulse  response.  Simulation  results  of  the  temporal  point-spread  function  of  a  diffu¬ 
sive  wave  are  provided,  as  well  as  experimental  results  of  a  fluorescence  decay  profile.  It  is  demonstrated 
that  our  new  time-resolved  technique  can  lead  to  high  signal-to-noise  ratios  and  short  data  acquisition  times, 
A  fluorescence-time-dependent  suppression  process  was  also  been  discovered.  ©  2002  Optical  Society  of 
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Time-resolved  techniques  are  important  in  optical 
instrumentation.  Measurements  of  time-dependent 
transmittance,  reflectance,  and  fluorescence  in  re¬ 
sponse  to  illumination  by  an  impulse  of  light  always 
contain  rich  information  that  can  be  used  to  retrieve 
the  dynamics  and  other  properties  of  the  sample 
under  investigation.  Time-resolved  photon-migration 
imaging  and  time-resolved  fluorescence  spectroscopy 
and  imaging  are  good  examples  of  such  techniques. 

Photon-migration  imaging  has  been  emerging  as  a 
potential  clinical  diagnostic  modality.^  Conventional 
time-resolved  photon-naigration  imaging  measures 
temporal  point-spread  functions  (TPSFs),  which  can  be 
used  for  quantitative  reconstruction  of  distributions 
of  scattering  and  absorption  coefficients.^  Accurate 
measurement  of  the  rising  edge  of  a  TPSF,  which 
is  related  to  that  quasi-straightforward  propagating 
component,  is  desirable  for  improving  spatial  resolu¬ 
tion.  Time-resolved  fluorescence  spectroscopy®  and 
imaging®  have  been  used  in  a  wide  range  of  research 
areas  such  as  biosciences,  chemistry,  and  clinical  di¬ 
agnosis.  The  fluorescence  lifetime  of  a  chromophore 
depends  on  its  characteristic  internal  structure  as 
well  as  on  its  physical  and  chemical  environments. 
Discrimination  of  signals  related  to  different  fluo¬ 
rescence  lifetimes  has  been  exploited  to  suppress 
background  autofluorescence,  improve  signal-to-noise 
ratio,  achieve  high  sensitivity  and  signal  specificity, 
and  measure  local  environmental  parameters  such  as 
pH,  temperature,  and  ion  concentrations.  Fluores¬ 
cence-lifetime  imaging  microscopy  makes  it  possible  to 
reveal  the  biomedical  processes  and  reactions  within 
a  cell  with  high  spatial  resolution.  Determination 
of  fluorescence  lifetime  is  generally  made  in  the  time 
domain,  in  which  short  excitation  pulses  are  used. 

Currently,  the  temporal  profile  of  the  response  to  an 
ultrashort  light  pulse  can  be  measured  with  either  a 
streak  camera  or  a  time-correlated  single-photon  count¬ 
ing  (TCSPC)  system.  TCSPC  is  preferred  because  its 
dynamic  range  and  temporal  linearity  are  better  than 
those  of  streak  cameras.  A  few  time-resolved  systems 
for  medical  optical  tomography  were  reported  recently, 
all  of  which  use  the  TCSPC  technique.*^’®  Nonetheless, 
the  principle  of  single-photon  counting  requires  that 
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no  more  than  one  photon  be  detected  in  each  cycle. 
This  restriction  causes  a  big  problem  in  data  acquisi¬ 
tion.  The  maximal  count  rate  is  limited  by  the  repe¬ 
tition  rate  of  laser  pulses,  the  processing  speed  of 
electronic  devices,  and  the  time  span  of  the  impulse 
response.  A  typical  count  rate  is  100,000  counts 
per  second,  which  means  that  approximately  15-30 
are  required  for  acquisition  of  10^-10®  photons  for 
one  temporal  profile.  In  some  applications,  multiple 
detection  channels  are  necessary  to  keep  the  overall 
data-acquisition  time  under  control.  For  example, 
Schmidt  et  al.  reported  on  a  32-channel  time-resolved 
instrument  that  needs  10-20  min  for  a  complete 
scan.®  It  is  worthy  of  mention  that  far  more  photon 
counts  are  needed  for  accurate  measurement  of  the 
intensity  of  early-arriving  photons. 

We  propose  in  this  Letter  a  coded  excitation  and 
correlation  detection  mechanism  that  borrows  ideas 
from  spread-spectrum  communications.  It  is  well 
known  that  a  spread  spectrum  communication  system 
possesses  many  desirable  properties,  such  as  selective 
addressing  capability,  low  error  rate,  and  interference 
rejection.®  A  broadband  pseudorandom  code  has 
only  a  weak  cross  correlation  with  other  codes  and 
an  autocorrelation  function  that  is  analogous  to  a 
delta  function.  Consequently,  a  receiving  system  can 
pick  up  the  correct  code  sequence,  addressing  it  from 
environmental  noise  and  interference,  and  is  able  to 
distinguish  the  same  sequence  that  arrives  at  different 
times  from  multiple  paths.  This  property  is  utilized 
in  our  time-resolved  optical  system,  in  which  photons 
detected  at  different  time  delays  need  to  be  resolved. 
Better  signal-to-noise  ratio,  shorter  data-acquisition 
time,  and  low  system  cost  are  the  expected  advantages 
of  our  new  technique.  Our  computer  simulations 
have  demonstrated  the  feasibility  of  a  time-resolved 
photon-migration  imaging  system  based  on  spread 
spectrum  excitation.  We  have  also  built  a  primitive 
time-resolved  f  luorometer  and  have  conducted  a  series 
of  experiments  with  it.  The  exponential  decay  profile 
of  a  long-lifetime  dye  was  successfully  retrieved. 

The  principle  of  our  method  is  simple  and  straight¬ 
forward.  We  denote  by  7(^)  the  time-dependent  re¬ 
sponse  of  a  sample  to  the  excitation  of  an  ultrashort 
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pulse.  A  light  source  continuously  modulated  with  a 
pseudorandom  bit  sequence  is  used  to  illuminate  the 
sample.  So  the  ac  component  in  the  detected  signal  is 
proportional  to  a  convolution  of  the  impulse  response 
and  the  excitation: 

R(t)^AI{t)*P{t),  (1) 

where  A  is  the  modulation  depth  and  P{t)  is  an 
N-bit-long  pseudorandom  bit  sequence.  P{t)  is  a 
maximal  length  sequence,  which  has  a  circular  auto¬ 
correlation  function  similar  to  a  delta  function: 

gi'r)  =  {p{t)p{t  -  t)> 

1  r/To  =  0 

-1/N  t/To  =  1,2,...,N/2;, 

^  -1,-2,...,  -N/2 

(2) 

in  which  r  is  the  time  delay  and  To  is  the  bit  pe¬ 
riod.  g{r)  increases  linearly  when  r  is  within  [-To,  0] 
and  decreases  linearly  in  [0,  To].  In  regions  out  of 
[-To,  To],  the  autocorrelation  values  are  essentially 
zero  when  N  is  big  enough.  By  correlating  R(t)  with 
P(t)  we  have 

f{r)^{R{t)P{t-r))^I{r)^g{T).  (3) 

Equation  (3)  is  valid  when  the  time  span  of  the  TPSF  is 
less  than  NTq.  If  To  is  small  enough,  fir)  at  a  specific 
delay  time  t  is  well  approximated  by  I{r).  When  To 
becomes  comparable  to  the  time  scale  of  7(t),  gir)  acts 
as  an  equivalent  temporal  gating  window.  So  the  im¬ 
pulse  response  or  its  integral  over  a  time  window  can 
be  retrieved  by  the  correlation  method. 

Time-resolved  measurement  of  diffusive  photon  den¬ 
sity  waves  by  the  spread  spectrum  approach  was  simu¬ 
lated  on  a  computer.  A  simulated  photon  migration 
imaging  system  consisted  virtually  of  a  four-channel 
12.5-Gbit/s  pattern  generator  (MP1775A,  Anritsu), 
a  10-Gbit/s  850-nm  vertical-cavity  surface-emitting 
laser  transmitter  (1780,  New  Focus),  a  12-GHz  re¬ 
ceiver  (1580-A,  New  Focus),  a  light-collecting  device 
that  effectively  increased  the  detection  aperture  to 
1  mm^,  and  electronic  components  such  as  mixers, 
filters,  and  amplifiers.  A  4095-bit  long  pseuodran- 
dom  bit  sequence  pit)  (To  =  100  ps)  generated  by  the 
pattern  generator  was  used  to  directly  modulate  the 
transmitter  to  generate  an  optical  pulse  sequence  with 
a  1-mW  peak-to-valley  value.  The  optical  fibers  of 
the  transmitter  and  the  receiver  were  embedded  in  an 
infinite  homogeneous  turbid  medium  and  were  sepa¬ 
rated  by  5  cm.  The  reduced  scattering  coefficient  of 
the  medium  was  6  cm"^  and  the  absorption  coefficient 
was  chosen  to  be  0.02  cm“^  The  ideal  response  of 
the  receiver  to  a  1-nJ  pulse  from  the  transmitter 
was  denoted  y(0  and  is  proportional  to  the  TPSF 
that  can  be  computed  with  the  diffusion  equation. 
The  receiver's  response  to  the  pseudorandom  normal 
sequence  stimulation  was  given  by 


M{t)  =  AVit)  *  [Pit)  +  UsRiit)]  +  ndR2it) ,  (4) 

where  A  —  10“^  and  where  Riit)  and  R2it)  were 
random  sequences  of  normal  distribution  with  a  unit 
standard  deviation.  Noise  levels  rig  =  0.022  and 
rid  =  2  were  derived  from  specifications  of  the 
transmitter  and  the  receiver.  The  correlation  can 
be  performed  by  a  mixer  (e.g.,  ZMX-IOG  from  Mini- 
Circuit)  and  a  low-pass  filter.  The  reconstructed 
Vit)  is  compared  in  Fig.  1  with  the  original  curve. 
The  noise  level  in  the  reconstructed  temporal  profile 
was  S.13  X  10“^,  or  3.2%  of  the  peak  value  of  y(^). 
Note  that  fewer  than  41  photons  could  be  detected  by 
TCSPC  within  the  same  measurement  time. 

To  validate  our  spread  spectrum  approach  ex¬ 
perimentally,  we  built  a  time-resolved  fluoroscopy 
system  and  conducted  experiments.  The  experimen¬ 
tal  system  consists  of  an  ultraviolet  LED  (NSHU550, 
Nichia  America  Corporation)  whose  output  spectrum 
is  centered  at  375  nm  and  a  photomultiplier  tube 
(R928,  Hamamatsu)  that  is  used  to  detect  fluorescent 
emissions.  A  personal  computer  is  responsible  for 
generating  the  same  pseudorandom  code  as  men¬ 
tioned  above,  synchronizing  input  and  output,  and 
processing  digital  signals.  A  train  of  pseudorandom 
voltage  signals  superimposed  upon  a  dc  signal  is 
generated  from  the  analog  output  port  of  a  mul¬ 
tifunction  data- acquisition  card  (PCI-MIO-16E-1, 
National  Instruments)  and  then  is  power  amplified  by 
a  current  driver  to  feed  the  LED.  The  update  rate 
for  the  bit  sequence  is  10®  bits/s,  and  the  sequence 
is  generated  repeatedly.  The  coded  excitation  light 
is  directed  to  a  fluorescent  sample  through  a  light 
guide.  The  fluorescent  photons  are  converted  by  the 
photomultiplier  tube  into  electrical  signals,  which  are 
amplified  by  a  40-dB  amplifier  before  being  digitized 
by  the  data-acquisition  card.  Digital  cross  correlation 
is  performed  on  the  computer  to  retrieve  the  impulse 
response,  i.e.,  the  temporal  fluorescence  decay  profile. 
A  solution  of  a  kind  of  europium  chelate  (Quantum 
Dye,  Research  Organics)  was  used  as  the  fluorescence 


Fig.  1.  Impulse  response  representing  the  TPSF  of  a  diffu¬ 
sive  photon  density  wave.  Thin  curve,  theoretical  predic¬ 
tion;  thick  curve,  reconstruction  from  noisy  environments. 
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Fig.  2.  Florescence  decay  temporal  profile  of  the  solution 
of  eropium  chelate,  Quantum  Dye. 

sample.  Quantum  Dye  has  a  long  lifetime  (more 
than  300  ^6s),  which  leads  to  moderate  requirements 
for  system-response  speed.  The  specified  excitation 
wavelength  ranges  from  360  to  390  nm,  and  the  emis¬ 
sion  has  a  peak  near  620  nm.  An  optical  long-pass 
filter  was  used  to  separate  the  emission  from  the 
excitation. 

An  acquired  decay  profile  is  plotted  in  Fig.  2.  Ini¬ 
tially,  the  fluorescence  intensity  decays  almost  expo¬ 
nentially.  However,  the  intensity  crosses  zero  near 
-^700  fjis  and  then  remains  negative  for  a  few  millisec¬ 
onds.  The  best  fit  of  profile  7(^)  with  two  exponential 
decay  terms  is  given  by 

I{t)  -  1.33  exp(-^/Ti)  -  0.389  exp(-^/T2) ,  (5) 

in  which  Ti  —  424  fxs  and  T2  =  1745  /xs.  The  second 
term  on  the  right-hand  side  of  Eq.  (5)  represents 
an  interesting  phenomenon,  which  we  call  time- 
dependent  suppression.  To  the  best  of  our  knowl¬ 
edge,  this  time-dependent  suppression  process  had 
not  been  discovered  previously.  A  simple  theoretical 
explanation  of  it  is  as  follows:  An  exciting  pulse  will 
bring  more  molecules  to  higher  energy  levels  and  thus 
temporarily  reduce  the  population  of  the  ground  state. 
Consequently,  the  excitation  rate  is  suppressed  before 
the  system  returns  to  prepulse  status.  From  the 
additional  time  constant  T2,  one  can  learn  more  about 
radiative  and  nonradiative  transitions  involved. 

The  temporal  profile  in  Fig.  2  was  an  average  over 
500  repeated  excitations,  which  took  20.5  s  in  total. 
The  acquisition  time  can  be  reduced  to  ^^4  s  if  a  bit 
rate  of  5  X  10®  bit/s  is  used.  The  standard  devia¬ 
tion  of  the  intensity  is  essentially  independent  of  time 


and  averages  0.022,  —2.3%  of  the  intensity  at  the  ori¬ 
gin  of  time.  Suppose  that  we  measure  a  fluorescence 
process  described  by  the  first  term  on  the  right-hand 
side  of  Eq.  (5)  with  an  ideal  TCSPC  system  and  that 
the  photon-counting  rate  is  set  to  be  200  counts/s.  If 
the  probability  that  a  photon  is  detected  in  each  circle 
(5  ms)  is  1,  the  relative  statistical  error  at  the  origin  of 
time  can  be  estimated  by 


where  C  is  the  total  number  of  photon  counts  and 
A  =  10  /xs  is  the  temporal  resolution.  For  example, 
we  have  err  =  2.1%  when  C  =  10®,  In  such  a  case 
the  total  data-acquisition  time  would  be  500  s.  Our 
method  can  reduce  the  data-acquisition  time  by  more 
than  tenfold,  or  significantly  improve  the  signal- 
to-noise  ratio  if  the  same  time  is  spent.  Of  course, 
no  one  will  really  use  single-photon  counting  for  such 
a  slow  process.  However,  the  above  comparison 
indicates  the  great  potential  of  the  spread  spectrum 
approach  in  its  future  applications  in  faster  optical 
processes,  such  as  fluorescence  decays  of  nanosecond 
lifetimes. 

To  conclude,  we  have  proposed  a  spread  spectrum 
approach  to  measuring  time-dependent  optical  pro¬ 
cesses.  Our  simulation  and  experimental  results 
have  demonstrated  that  this  new  method  can  remark¬ 
ably  reduce  data-acquisition  time  and  (or)  improve 
overall  signal-to-noise  ratio.  In  some  situations,  a 
much  lower  system  cost  can  also  be  expected. 

This  research  was  supported  in  part  by  U.S.  Depart¬ 
ment  of  Defense  Army  Breast  Cancer  Research  Pro¬ 
grams  (grants  DAMD17-00- 1-0217  and  DAMD17-01~ 
1-0216).  N.  G.  Chen’s  e-mail  address  is  chenng@engr. 
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Appendix  F 


Simultaneous  reconstruction  of  absorption  and 
scattering  maps  with  ultrasound  localization: 
feasibility  study  using  transmission  geometry 


Minming  Huang,  Tuqiang  Xie,  Nan  Guang  Chen,  and  Quing  Zhu 


We  report  the  experimental  results  of  the  simultaneous  reconstruction  of  absorption  and  scattering 
coefficient  maps  with  ultrasound  localization.  Near-infrared  (NIR)  data  were  obtained  from  frequency 
domain  and  dc  systems  with  source  and  detector  fibers  configured  in  transmission  geometry.  High-  or 
low-contrast  targets  located  close  to  either  the  boundary  or  the  center  of  the  turbid  medium  were 
reconstructed  by  using  NIR  data  only  and  NIR  data  with  ultrasound  localization.  Results  show  that  the 
mean  reconstructed  absorption  coefficient  and  the  spatial  distribution  of  the  absorption  map  have  been 
improved  significantly  with  ultrasound  localization.  The  improvements  in  the  mean  scattering  coeffi¬ 
cient  and  the  spatial  distribution  of  the  scattering  coefficient  are  moderate.  When  both  the  absorption 
and  the  scattering  coefficients  are  reconstructed  the  performance  of  the  frequency-domain  system  is  much 
better  than  that  of  the  dc  system.  ©  2003  Optical  Society  of  America 
OCIS  codes:  170.0170,  170.3010,  170.5270,  170.7170,  170.3830. 


1.  Introduction 

Functional  imaging  with  near-infrared  (NIR)  dif¬ 
fused  light  has  found  potential  applications  in  many 
areas.  However,  one  fundamental  problem  that 
we  have  to  overcome  is  the  intense  scattering  of  light. 
As  a  result,  diffused  light  probes  a  widespread  region 
instead  of  providing  information  along  a  straight  line. 
Multiple  measurements  are  always  correlated  as  a 
result  of  the  overlapping  of  probed  regions.  There¬ 
fore  increasing  the  total  number  of  measurements 
does  not  necessarily  provide  more  independent  infor¬ 
mation  for  image  reconstruction.  In  general  the  in¬ 
verse  image  reconstruction  is  underdetermined*  and 
ill-posed.  The  behavior  of  reconstruction  algorithms 
is  affected  by  many  factors,  such  as  the  system  signal- 
to-noise  ratio,  the  probe  configuration,  and  regulation 
schemes  used  in  image  reconstructions. 

The  current  approaches  to  image-reconstruction  al¬ 
gorithms  are  (1)  simple  backprojection  methods, (2) 
perturbation  methods, ^nd  (3)  finite-element 
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methods  (FEMs).^3-i5  ^he  backprojection  method 
provides  a  real-time  estimation  of  the  coarse  optical 
properties  of  lesions.  However,  the  reconstructed- 
image  resolution  is  low  and  the  lesions  that  appear  in 
the  images  are  often  displaced  from  true  locations. 
Perturbation  methods  are  in  general  based  on  linear 
approximations  to  the  heterogeneous  functions. 
Bom  and  Rytov  approximations  are  examples. 
Measurements  between  the  background  and  the  het¬ 
erogeneous  medium  are  used  to  relate  the  optical 
signals  at  the  measurement  surface  to  absorption 
and  scattering  variations  in  each  volume  element 
within  the  sample.  The  least-squares  method  is  in 
general  used  to  formulate  the  inverse  problem.  It¬ 
erative  search  methods,  such  as  conjugate-gradient 
techniques,  are  employed  to  iteratively  solve  the  in¬ 
verse  problem.  However,  an  accurate  estimation  of 
the  target  optical  properties  depends  on  an  accurate 
estimation  of  the  background  optical  properties, 
which  are  in  general  not  easy  to  obtain  in  breast 
tissues.  In  addition,  when  the  absorption  or  the 
scattering  coefficients  of  the  lesions  are  significantly 
higher  than  the  background,  the  linear  perturbation 
methods  cannot  give  accurate  optical  properties. 
Reconstructions  based  on  FEMs  provide  higher-order 
estimations  to  heterogeneous  functions.  However, 
as  the  scattering  and  absorption  coefficients  are  ex¬ 
panded  over  local  basis  functions,  the  number  of  un¬ 
knowns  is  increased  considerably.  Nonetheless  the 
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independent  information  from  diffused  light  is  lim¬ 
ited,  and  it  essentially  stops  increasing  when  the 
source- detector  pairs  reach  a  certain  number.  In 
this  case  the  inverse  problem  becomes  more  ill-posed 
and  the  behavior  of  the  reconstruction  algorithm  may 
be  unpredictable. 

Reconstructions  with  the  aid  of  a  priori  target  ge¬ 
ometry  information  provided  by  coregistered  ultra¬ 
sound  have  shown  promising  results  in  improving  the 
accuracy  of  reconstructed  optical  properties  and  the 
localization  of  targets.  ^^-20  In  this  method  an  ultra¬ 
sound  probe  and  NIR  source  and  detector  fibers  have 
been  deployed  on  a  handheld  probe  and  configured  in 
the  reflection  geometry.  The  a  priori  tissue  type  and 
lesion  location  as  well  as  shape  provided  by  coregis¬ 
tered  ultrasound  can  guide  the  NIR  image  recon¬ 
struction  to  localized  target  regions.  Therefore  the 
total  number  of  voxels  with  unknown  optical  proper¬ 
ties  can  be  reduced  significantly  and  the  inversion 
becomes  well  defined.  In  general  the  solution  is 
unique  and  the  iterative  search  algorithms  converge 
very  fast. 

In  this  paper  we  explore  experimentally  the  utility 
of  the  combined  imaging  with  NIR  sources  and  de¬ 
tectors  configured  in  transmission  geometry.  The 
standard  pulse-echo  technique  is  used  to  obtain 
coregistered  ultrasound  images.  Compared  with 
NIR  reflection  geometry,  transmission  geometry  may 
have  the  advantages  of  a  lower  requirement  on  the 
system  dynamic  range  and  higher  sensitivities.  In 
the  reported  experiments  the  FEM  forward  solver  is 
used  to  generate  a  Jacobian  weight  matrix  for  simul¬ 
taneous  reconstruction  of  the  absorption  and  the 
scattering  coefficient  maps,  and  the  inversion  is  per¬ 
formed  by  using  NIR  data  only  and  NIR  data  with  a 
priori  target  geometry  provided  by  coregistered  ul¬ 
trasound. 

Frequency-domain  and  dc  systems  are  widely  used 
by  many  research  groups.  The  advantage  of  using  a 
frequency-domain  system  is  that  both  amplitude  and 
phase  information  can  be  obtained,  and  therefore  si¬ 
multaneous  reconstruction  of  absorption  and  scatter¬ 
ing  coefficients  is  feasible.  The  disadvantage  is  the 
complexity  in  constructing  high-frequency  transmis¬ 
sion  and  reception  circuitries,  particularly  when 
many  parallel  channels  are  needed  for  fast  data  ac¬ 
quisition.  This  situation  is  more  pronounced  when 
nearly  100  source  and  detector  positions  are  needed 
for  three-dimensional  (3D)  imaging  in  transmission 
geometry.  In  addition  the  system  cost  is  high. 
Direct-current  systems  have  the  advantages  of  sim¬ 
plicity  in  circuit  construction  and  low  cost.  The 
drawbacks  are  that  the  phase  information  is  not 
available  and  cross  talk  between  the  reconstructed 
absorption  and  scattering  coefficients  may  be  un¬ 
avoidable.  In  this  paper  we  compare  the  perfor¬ 
mance  of  both  systems  in  simultaneously 
constructing  absorption  and  diffusion  coefficient 
maps  in  conditions  of  NIR  reconstruction  only  and 
ultrasound  guided  NIR  reconstruction. 


2.  Forward  Model 

A  commercial  FEM  package  femlab  has  been  used  as 
a  forward  solver  to  generate  the  Jacobian  weight  ma¬ 
trix.  This  software  has  an  interactive  environment 
for  modeling  and  simulating  scientific  and  engineer¬ 
ing  problems  based  on  partial  differential  equations. 
Users  can  utilize  various  existing  models  by  defining 
relevant  parameters  and  can  also  design  specific 
models  based  on  equations  under  consideration. 

In  FEMLAB  the  scalar-governing  equation  and 
boundary  condition  are  given  as 

rfa^-V-[c-Vf/(r)  +  a[/(rh] 
ot 

+  aU{f)  +  3Vt/(P)  =  f{f)  in  a,  (1) 

n‘[cVC/(r)  +  aU(r)7]  +  gf7(r) 

=  g{f)  on  boimdary  afl.  (2) 

n  is  the  domain  of  interest  where  Eq.  (1)  is  satisfied, 
and  an  is  the  boundary  of  fl  where  Eq.  (2),  represent¬ 
ing  a  generalized  Neumann  boundary  condition,  is 
satisfied.  f/(r)  is  the  unknown  scalar  to  be  solved,  f 
is  the  source  term  in  fl,  c  is  the  diffusion  coefficient, 
a  is  the  absorption  coefficient,  a  is  the  conservative 
flux  convection  coefficient,  d^,  is  the  mass  coefficient, 
P  is  the  convection  coefficient,  y  is  the  conservative 
flux  source  term,  g(r)  is  the  source  at  the  boundary,  n 
is  the  outward  unit  normal  vector  on  and  9  is  a 
positive  number  related  to  the  internal  reflection  at 
the  boundary.  The  time-domain  photon-diffusion 
equation  can  be  obtained  from  Eq,  (1)  as 

1  a4>(^  _  ^ ^  ^ 

V  dt 

(3) 

where  U  is  replaced  with  <l)(r,  t),  the  photon  density  at 
location  r;  f  with  t)y  the  point  source  term;  a  with 

|x^,  the  optical  absorption  coefficient;  c  with  D,  the 
diffusion  coefficient  [D  =  l/(3p,5')j  where  is  the 
reduced  scattering  coefficient];  and  with  1/v, 
where  v  is  the  speed  of  light  in  the  medium.  All 
other  coefficients  are  chosen  to  be  zero. 

By  letting  g(r)  equal  zero,  we  can  obtain  the 
Neumann-type  boundary  condition  as 

n  •D(r)  •  V4)(r,  0  +  0  =  0-  (4) 

By  setting  d^  to  zero,  we  can  obtain  the  steady-state 
diffusion  equation  as 

V  •  D(r)  •  V4)(r)  -  ^Mr)  =  Soir).  (5) 

Furthermore,  by  replacing  with  jx^  -  m/v,  where 
CO  is  the  modulation  frequency,  we  can  obtain  a 
frequency-domain  diffusion  equation  as 

(i(o\ 

=  -Soir,  (o).  (6) 
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Assume  that  the  changes  due  to  the  absorption  and 
the  scattering  inhomogeneities  are  small.  The  Ja¬ 
cobian  matrix 


W  = 


[Am.,/  ADj\ 


that  relates  the  photon-density  perturbation  at 
boundary  cell  i  and  imaging  voxel  j  with  absorption 
coefficient 


TO 


change 

Ajjl  and 

I  difiusion 

coefficient 

)y,  is  given  as 

A<j)ii 

A^ii 

A(|)ii 

A<|)u, 

Ajtal 

A(J.o£, 

ADi 

AOi, 

A<|)21 

A4)2£, 

A4)21 

A<|)2I, 

A|a„i 

A(iat 

A^i 

ADt 

1 

a4mi 

Ac^mz, 

A<^mi 

A<^mz. 

.A|x„i 

A|x„t 

ADi 

adJ 

(7) 


mesh.  As  a  result  the  total  number  of  unknown  vox¬ 
els  is  400,  which  is  comparable  with  the  total  number 
of  measurements. 

In  the  Jacobian  matrix  calculation  we  evaluated 
the  linearity  of  A(I)jy/AiXa^-  versus  Ap-^y-  with 
changing  from  0.02  to  0.3  cm”^  (the  typical  range  of 
the  absorption  coefficient  for  breast  lesions)  at 
many  cell  locations.  For  dc  component  calcula¬ 
tions  the  error  due  to  the  linear  approximation  is 
always  less  than  0.4%.  For  frequency-domain  cal¬ 
culations  the  real  and  imaginary  parts  of  4)y-  of  sev¬ 
eral  cells  are  1  or  2  orders  smaller,  and  they  are 
very  sensitive  to  small  perturbations.  However, 
even  for  the  worst  case,  the  error  due  to  the  use  of 
the  linear  approximation  is  also  less  than  0.5%  for 
24  out  of  28  detector  locations.  Therefore  we  can 
simplify  the  Jacobian  matrix  calculation  by  comput¬ 
ing  at  two  points  I^oh)  f^at  are  be¬ 

tween  0.02  and  0.3  cm'*^.  The  selection  of  these 
two  points  is  based  on  the  following  considerations. 
The  linear  slope  A(t)y/AM.aj  is  best  estimated  when  a 
smaller  perturbation  Apu^y  is  used.  However,  a 
ApLq;  that  is  too  small  is  not  suitable  for  computing 
As  a  good  compromise  we  chose  =  0.02 
cm“^  and  cm”^  for  all  the  calculations. 

We  also  evaluated  the  linearity  of  A^ij/ADj  versus 
ADj  with  Dj  changing  from  0.0556  (the  correspond¬ 
ing  M's'  is  6  cm”^)  to  0.0306  cm  (the  corresponding 
jXg'  is  10.9  cm“^)  (the  typical  diffusion  coefficient 
range  for  breast  lesions)  at  many  cell  locations. 
For  both  the  dc  component  and  frequency-domain 
calculations  the  error  due  to  the  use  of  a  linear 
approximation  is  always  less  than  0.2%.  There¬ 
fore  we  compute  A(j>y/AD^-  by  using  two  points  at 
0.0556  and  0.0356  cm. 

The  experimental  data  included  many  unknown 
systematic  factors  (unknown  source,  detector  gains, 
etc.).  To  remove  those  unknown  factors,  we  used  the 
normalized  difference  method  in  Ref.  22  as 


4>m(t)  4>mr(t) 


4>cr(t) 


=  [Wy]-[Ato  AD,,],  (8) 


where  M  is  the  total  number  of  boundary  cells  and  L 
is  the  total  number  of  imaging  voxels. 

3.  Computation  Procedures 

We  have  generated  a  two-dimensional  (2D)  mesh 
with  4050  triangle  elements  and  2088  nodes  (Fig.  1). 
The  radius  of  the  physical  boundary  is  chosen  to  be 
45.26  cm  and  the  extrapolation  distance  =  0.7  • 

[1^^  =  l/(p-s'  +  P'J]  is  calculated  to  be  0.1169  cm,2i 
where  is  the  background  value  and  is  chosen 

to  be  6  cm^^  (a  typical  value  for  normal  breast  tissue). 
To  speed  up  the  Jacobian  matrix  calculation,  we 
chose  a  coarse  mesh  of  5  cm  X  5  cm  in  the  middle  of 
the  2D  mesh  by  assuming  that  there  are  no  pertur¬ 
bations  outside  this  area.  The  imaging  pixel  size  is 
chosen  to  be  0.25  cm  X  0.25  cm  inside  the  coarse 


where  (|)„i(,)  is  the  measured  heterogeneous  data  as¬ 
sociated  with  the  source- detector  pair  i  with  the  tar¬ 
get  in  the  homogeneous  medium,  ^^nrn)  is  the 
measured  homogeneous  data,  <^cr(i)  is  the  calculated 
forward  homogeneous  data  used  for  the  weight  ma¬ 
trix  calculation.  The  (}>  is  a  complex  value  in  the 
frequency-domain  case.  The  unknown  system  fac¬ 
tors  present  in  both  sets  of  measured  heterogeneous 
and  homogeneous  data  are  canceled  by  taking 
the  ratio  of  the  perturbation  (heterogeneous- 
homogeneous)  to  the  reference  (homogeneous)  mea¬ 
surements. 

To  alleviate  the  cross  talk  between  the  absorption 
coefficient  and  the  diffusion  coefficient,  we  also  nor¬ 
malized  each  column  and  row  of  the  weight  matrix  by 
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their  own  mean  value  to  improve  the  matrix  condi¬ 
tion.  Thus  the  inverse  problem  is  modified  as 


[Fu]' 


^cr(0 


-[Fu]-[W,]'[Gj,^-[AXjil  (9) 


where  and 


2U0^ 


I  =  1,  2,  . . .  ,  iW  (Fy  =  0  when  i  i^j), 


/  1  "  \ 


>  =  1,  2,  .  . .  ,  L  (Gy  =  0  when  i  ¥=j). 


Finally  the  conjugate  gradient  method  is  used  for 
solving  Eq.  (9).  Since  the  error  function 

F- W-G- AX 

reduces  consistently  and  becomes  flat  after  certain 
iterations  (around  30  for  reconstruction  with  the  NIR 
data  only  and  3  for  the  NIR  with  ultrasound  guid¬ 
ance),  the  stopping  criteria  are  chosen  to  be  100  and 
3  iterations  for  the  NIR  only  and  the  NIR  with  ultra¬ 
sound  localization,  respectively.  The  corresponding 
computation  time  is  around  2  min  for  NIR  data  only 
and  30  s  for  NIR*  with  ultrasound  information. 

In  the  experiments,  data  are  obtained  from  a  3D 
model  (the  finite  cylindrical  medium  and  the  finite 
length  target).  However,  the  Jacobian  matrix  com¬ 
putation  is  based  on  a  2D  model  (the  infinite  cylin¬ 
drical  medium  and  the  infinite-length  target). 
According  to  Refs.  23  and  24  the  2D/3D  difference 
between  intensity  is  reasonably  independent  of 
source- detection  separation  and  it  is  reasonable  to 
correct  it  by  simply  multiplying  or  dividing  a  con¬ 
stant.  Therefore  the  normalized  difference  method 
not  only  cancels  the  unknown  source- detector  gains 
but  also  partially  eliminates  the  2D/3D  model  mis¬ 
match  by  taking  the  ratio  of  the  perturbation  and  the 
reference  measurements. 


4.  Experimental  Systems 

To  compare  the  performance  of  the  frequency-domain 
and  the  dc  systems  with  and  without  ultrasound  lo¬ 
calization,  we  constructed  both  systems.  Details  of 
our  multisource  and  multidetector  frequency-domain 
system  can  be  found  in  Ref  18.  Briefly,  it  consists  of 
12  pairs  of  dual-wavelength  sources  [780  nm  (LT024) 
and  830  nm  (LT015),  Sharp  Corporation]  modulated 
at  140  MHz  and  8  parallel  detectors  (R928, 
Hamamatsu  Corporation).  These  sources  and  de¬ 
tectors  are  coupled  to  the  medium  through  optical 
fibers. 

Our  new  dc  system  consists  of  six  pairs  of  dual¬ 
wavelength  laser  diodes  (780  and  830  nm)  amplitude 
modulated  at  20  kHz  to  avoid  dc  fluctuations  (Fig.  2). 


Fig.  2.  Schematic  of  our  new  dc  system:  AO,  analog  output;  AI, 
analog  input;  DO,  digital  output;  DI,  digital  input. 


Twelve  detector  fibers  are  sequentially  coupled  to  a 
photomultiplier  tube  (PMT)  (R928,  Hamamatsu  Cor¬ 
poration)  through  a  mechanical  multiplexer  that  uses 
a  stepper  motor  to  collimate  accurately  a  coupling 
light  guide  on  a  rotation  plate  to  one  of  the  detection 
fibers.  Each  light  source  has  two  output  levels  of  30 
dB  in  difference  to  avoid  PMT  saturation  and  to  im¬ 
prove  the  system  dynamic  range.  A  National  In¬ 
strumentation  data  acquisition  (DAQ)  card  is  used  to 
generate  sinusoidal  waveforms  on  top  of  a  constant 
current  to  drive  one  laser  diode  at  a  time.  The  pho¬ 
ton  density  wave  detected  by  the  PMT  is  amplified 
and  then  digitized  by  the  same  DAQ  card.  The  cor¬ 
responding  amplitude  is  retrieved  from  the  acquired 
waveform  on  a  PC.  The  digital  ports  of  the  DAQ 
card  are  used  to  switch  laser  diodes  sequentially, 
control  the  output  level,  and  input  a  feedback  signal 
for  closed-loop  stepper  motor  control.  The  total 
DAQ  time  for  a  complete  data  set  is  around  1  min. 

Since  our  dc  system  is  modulated  at  20  kHz,  the 
measured  amplitude  is  an  ac  signal  instead  of  the  dc 
component  that  we  used  for  forward  Jacobian  matrix 
calculation.  We  estimated  the  difference  by  using  a 
spherical  wave  of  light  energy  density  [exp(j^p)/Dp] 
of  two  wave  numbers.  ki  =  (-jx^/Z?)^^^  and  ^2  ” 
[(-|x^  +  j  X  (o/v)/Z)]^^^,  where  w  =  (277)2  X  10^  Hz 
and  i;  and  D  are  the  speed  of  light  and  the  diffusion 
coefficient,  respectively.  The  normalized  difference 
[expC/^iP)  -  exp(j^2p)]/exp(j/2ip)  can  be  approxi¬ 
mated  as 

-■’-p 

2v 

and  it  is  larger  when  the  wave  propagates  deeper  into 
the  tissue  and  the  background  absorption  is  smaller. 
For  the  worst  case  of  =  0.01  cm""^  and  p  =  10  cm 
the  difference  is  less  than  3.295”“®.  Therefore  we 
directly  use  the  Jacobian  matrix  calculated  from  the 
steady-state  diffusion  equation  [Eq.  (5)]  to  relate  the 
measured  amplitude  data  for  the  inversion. 

The  experimental  setup  is  shown  in  Fig.  3.  Dur¬ 
ing  experiments,  source  and  detector  fibers  were 
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Fig.  3.  Experimental  setup.  A  commercial  ultrasound  probe  is 
located  at  the  top  of  a  water  tank,  and  the  NIR  source  and  detector 
fibers  are  deployed  armmd  the  tank  and  configured  in  transmission 
geometry. 


placed  on  a  circular  plane  (Fig.  4).  A  commercial 
ultrasound  probe  of  3.5-MHz  central  frequency 
(HP77020A  ultrasound  imaging  system)  was  placed 
on  top  of  a  water  tank  filled  with  Intralipid.  This 
one-dimensional  commercial  probe  provides  a  cross- 
section  image  in  the  y-z  imaging  plane  (called  the 
B-scan),  where  y  is  the  lateral  direction  and  z  is  the 
propagation  direction.  By  translating  the  probe  me¬ 
chanically  in  the  x  direction,  we  acquired  3D  volu¬ 
metric  image  data.  Windowing  the  3D  data  in  the  z 
direction  at  a  particular  depth  provides  2D  target 
spatial  images  at  the  x-y  plane  (called  the  C-scan). 
C-scans  are  coregistered  with  the  NIR  images.  An 
Intralipid  solution  of  0.6%  was  used  as  a  homoge¬ 
neous  background.  The  fitted  and  of  the  In¬ 
tralipid  at  780  nm  were  0.02  and  6  cm“^,  respectively. 


Since  we  have  a  limited  number  of  source  and  de¬ 
tector  positions,  we  placed  8  source  fibers  instead  of 
12  on  a  quarter  circle  and  8  detector  fibers  on  the 
opposite  quarter  for  frequency-domain  measure¬ 
ments  [Fig.  4(a)].  For  dc  measurements  we  placed  6 
source  fibers  on  a  quarter  circle  and  12  detector  fibers 
on  the  opposite  quarter  [see  Fig.  4(b)].  After  obtain¬ 
ing  one  set  of  heterogeneous  NIR  data  with  a  target 
in  the  Intralipid  (8x8  measurements  for  the  fre¬ 
quency  system,  6  X  12  measurements  for  the  dc  sys¬ 
tem),  we  simply  rotated  the  target  by  90°  for  the 
frequency-domain  system  and  77.14°  for  the  dc  sys¬ 
tem  to  get  another  set  of  data.  Figure  4(c)  shows  the 
rotation  scheme  used  during  frequency-domain  ex¬ 
periments.  Three  rotations  were  performed  to  ob¬ 
tain  one  complete  data  set  for  reconstruction,  which 
covered  360°.  For  dc  measurements  a  similar  rota¬ 
tion  scheme  was  used  and  four  rotations  were  per¬ 
formed  to  obtain  one  complete  data  set.  Three- 
dimensional  ultrasound  data  were  obtained 
simultaneously  with  NIR  data  acquisition. 

5.  Results 

For  NIR  reconstruction  with  transmission  geometry 
the  perturbations  in  caused  by  the  target  are  re¬ 
lated  to  the  location  of  the  target  and  the  contrast 
between  the  target  and  the  background.  It  is  much 
harder  to  reconstruct  a  low-contrast  target  located  in 
the  middle  of  the  medium.  In  our  experiments  a 
cylindrical  target  (high/low  contrast)  1.0  cm  in  diam¬ 
eter  was  placed  at  two  typical  locations  (at  the  center 
and  close  to  the  boundary).  The  target  was  made  of 
polyester  resin  with  calibrated  (jl^  and  Based 

on  the  calibration  results  from  Dartmouth  College, 
the  absorption  coefficient  of  the  target  was  from  0.20 
to  0.3  cm”^  for  the  high-contrast  target  case  and  was 
0.07  cm"^  for  the  low-contrast  target;  of  both 
targets  was  9.6  cm“^.  The  repeatability  of  the  cali¬ 
bration  was  within  ±5%. 

A.  Reconstruction  when  the  Frequency-Domain  System 
is  Used 

For  the  high-contrast  target  case,  improvements  in 
the  reconstructed  pL^^  are  dramatic  when  ultrasound 
localization  is  used.  Figures  5(a)  and  5(b)  show  the 
simultaneously  reconstructed  absorption  and  diffu- 


(a)  (b)  (c) 

Fig.  4.  Configurations  of  NIR  sources  and  detectors  used  in  the  reported  experiments:  (a)  frequency-domain  experiments,  (b)  dc 
experiments,  (c)  top  view  of  target  rotation  scheme  used  in  frequency-domain  measurements. 
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Fig.  5.  Comparisons  of  simultaneously  reconstructed  absorption  and  diffusion  coefficient  maps  of  a  high-contrast  target  located  at  the 
center  of  the  turbid  medium:  (a),  (b)  reconstructed  |x„  and  D  distributions  when  NIR  data  only  are  used;  (c),  (d),  reconstructed  and  D 
distributions  with  ultrasound  localization;  (e)  coregistered  ultrasound  C-scan  image  used  to  guide  the  NIR  reconstruction.  X  and  Y  are 
the  spatial  dimensions  in  centemeters. 


sion  maps  where  NIR  data  only  is  used  when  the 
target  is  located  at  the  center  of  the  medium,  while 
Figs.  5(c)  and  5(d)  show  the  reconstructed  absorption 
and  diffusion  maps  with  ultrasound  localization. 
Figure  5(e)  is  the  C-scan  ultrasound  image  obtained 


at  the  source- detector  plane.  A  threshold  was  ap¬ 
plied  to  the  ultrasound  image,  which  was  chosen  as 
the  FWHM  from  the  peak  value.  The  pixels  with 
gray-scale  values  greater  than  the  threshold  were 
mapped  out  as  a  target  region,  while  the  pixels  with 
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Fig.  6.  Comparisons  of  simultaneously  reconstructed  absorption  and  diffusion  coefficient  maps  of  a  high-contrast  target  located  closer  to 
the  boundary  of  the  turbid  medium:  (a),  (b)  reconstructed  and  D  distributions  when  NIR  data  only  are  used;  (c),  (d),  reconstructed 
and  D  distributions  with  ultrasound  localization;  (e)  coregistered  ultrasound  C-scan  image  used  to  guide  the  NIR  reconstruction. 


values  less  than  the  threshold  were  identified  as  a 
nontarget  region.  The  inverse  reconstruction  was 
localized  to  the  target  region.  With  ultrasound  the 
mean  reconstructed  absorption  coefficient  was  im¬ 
proved  from  32.4%  to  64.4%  of  the  calibrated  mean 
value  and  the  FWHM  was  improved  from  140%  to 


100%  of  the  true  target  size.  The  mean  value  was 
calculated  from  the  pixels  with  absorption  coeffi¬ 
cients  greater  than  the  FWHM  in  the  target  mass 
region.  As  for  the  diffusion  coefficient  the  mean  re¬ 
constructed  values  were  within  5%  of  the  calibrated 
value  both  with  and  without  ultrasound  localization. 
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(c)  (d) 

Fig.  7.  Comparisons  of  simultaneously  reconstructed  absorption  and  diffusion  coefficient  maps  of  a  low-contrast  target  located  closer  to 
the  boundary:  (a),  (b)  simultaneously  reconstructed  and  D  distributions  when  NIR  data  only  are  used;  (c),  (d),  simultaneously 
reconstructed  and  D  distributions  with  ultrasound  localization. 


However,  the  edge  artifacts  seen  in  Fig.  5(b)  were 
eliminated  by  using  ultrasound  guidance. 

Figure  6  shows  the  experimental  results  of  the 
same  high-contrast  target  located  close  to  the 
boundary.  Figure  6(a)  is  the  distribution  recon¬ 
structed  by  using  NIR  data  only.  The  recon¬ 
structed  mean  absorption  coefficient  was  —45.2%  of 
the  calibrated  mean  value,  and  the  FWHM  was 
around  1.3  cm,  which  was  130%  of  the  true  target 
size.  Figure  6(b)  is  the  D  distribution  simulta¬ 
neously  reconstructed  with  the  map  by  using 
NIR  data  only.  Artifacts  in  the  boundary  of  D  dis¬ 
tribution  are  visible.  The  reconstructed  mean 
value  was  -65.2%  of  the  calibrated  value,  and  the 
FWHM  was  around  1.2  cm.  Figures  6(c)  and  6(d) 
are  the  simultaneously  reconstructed  and  D  dis¬ 
tributions  with  coregistered  ultrasound  localiza¬ 
tion.  Figure  6(e)  is  the  coregistered  C-scan 
ultrasound  image  obtained  in  the  source- detector 
plane.  With  ultrasound  localization  the  mean 


was  improved  to  65.2%  of  the  calibrated  mean  value 
and  the  FWHM  was  improved  to  1  cm.  The  mean 
D  was  improved  to  102%  of  the  calibrated  value. 
Artifacts  in  the  D  distribution  were  eliminated  with 
ultrasound  localization. 

Figure  7  shows  the  experimental  results  of  the  low- 
contrast  target  located  close  to  the  boundary.  Fig¬ 
ure  7(a)  is  the  distribution  reconstructed  by  using 
NIR  data  only.  The  reconstructed  mean  absorption 
coefficient  was  —93%  of  the  calibrated  value,  and  the 
FWHM  was  around  1.3  cm.  Figure  7(b)  is  the  D 
distribution  simultaneously  reconstructed  with  the 
|x^  map  by  using  NIR  data  only.  The  reconstructed 
mean  value  was  -81%  of  the  calibrated  value,  and 
the  FWHM  was  around  1.2  cm.  Figures  7(c)  and 
7(d)  are  the  simultaneously  reconstructed  |x^  and  D 
distributions  with  coregistered  ultrasound  localiza¬ 
tion.  With  ultrasound  localization  the  mean  fx^  was 
within  11%  of  the  calibrated  value  and  the  FWHM 
was  improved  to  1  cm.  The  mean  D  was  improved  to 
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Fig.  8.  Comparisons  of  simultaneously  reconstructed  absorption  and  diffusion  maps  of  the  same  low-contrast  target  located  at  the  center 
of  the  turbid  medium:  (a),  (b),  reconstructed  and  D  distributions  when  NIR  data  only  are  used;  (c),  (d),  reconstructed  and  D 
distributions  with  ultrasound  localization. 


102%  of  the  calibrated  value.  The  edge  artifacts 
seen  in  the  D  distribution  were  eliminated  with  ul¬ 
trasound  localization. 

Figure  8(a)  and  8(b)  show  the  simultaneously  re¬ 
constructed  fjL^  and  D  distributions  of  the  same  low- 
contrast  target  located  at  the  center.  As  one  can 
see  the  spatial  distributions  of  |ul^  and  D  are  worse 
than  those  shown  in  Figs.  7(a)  and  7(b).  Artifacts 
in  the  boundary  of  the  D  distribution  are  visible. 


The  reconstructed  mean  jjL^  was  ~94%  of  the  cali¬ 
brated  value,  and  the  reconstructed  mean  D  was 
—92%  of  the  calibrated  value.  Figures  8(c)  and 
8(d)  are  the  simultaneously  reconstructed  absorp¬ 
tion  and  diffusion  coefficient  maps  with  ultrasound 
localization.  The  reconstructed  mean  is  within 
13%  of  the  calibrated  value  and  the  reconstructed 
mean  D  is  within  6%  of  the  calibrated  value.  The 
FWHM  of  the  map  has  been  improved  by  30%, 


Table  1.  Experimental  Results  from  the  Frequency-Domain  System 


NIR  Only 

NIR  +  US 

Calibrated  Value  from 
Dartmouth  College 

Target  Location 

Target 

Contrast 

Mean 

Mean  D 
(cm) 

Mean 

(cm"') 

Mean  D 
(cm) 

Mean  p.^ 
(cm"') 

Mean  D 
(cm) 

Close  to  the  boundary 

At  the  center 

High  contrast 
Low  contrast 
High  contrast 
Low  contrast 

0.113  (45.2%“) 
0.065  (93%) 
0.081  (32.4%®) 
0.067  (94%) 

0.030  (87%) 
0.028  (81%) 
0.0353  (101%) 
0.032  (92%) 

0.163  (65.2%“) 
0.078  (111%) 
0.161  (64.4%“) 
0.0793  (113%) 

0.0355  (102%) 
0.0355  (102%) 
0.0363  (105%) 
0.037  (106%) 

0.2-0.3 
0.07  ±  5% 
0.2-0.3 
0.07  ±  5% 

0.0347  ±  5% 
0.0347  ±  5% 
0.0347  ±  5% 
0.0347  ±  5% 

“Tlie  percentage  is  calculated  based  on  the  mean  value  of  0.2- 0.3  cm 
US,  ultrasound 
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Fig.  9.  Comparisons  of  simultaneously  reconstructed  absorption  and  diffusion  maps  of  a  high-contrast  target  located  at  the  center  of  the 
medium.  The  dc  system  of  the  configuration  in  Fig.  4(b)  was  used  for  the  experiments:  (a),  (b),  reconstructed  (r„  and  D  distributions 
when  NIR  data  only  are  used;  (c),  (d),  reconstructed  p,„  and  D  distributions  with  ultrasound  localization;  (e)  coregistered  C-scan  ultrasound 
image  used  to  guide  the  NIR  reconstruction. 


and  artifacts  in  the  D  distribution  have  been  elim¬ 
inated. 

In  Table  1  are  summaries  of  the  experimental  re¬ 
sults  of  the  two  cases  in  which  the  target  was  located 


at  different  positions.  In  general,  when  the  target 
was  located  at  the  center,  the  results  with  NIR  data 
only  were  poorer  than  those  when  the  target  was 
located  closer  to  the  boundary.  With  ultrasound  lo- 
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(c)  (d) 

Fig.  10.  Comparisons  of  simultaneously  reconstructed  absorption  and  diffusion  coefficient  maps  of  a  low-contrast  target  located  at  the 
center  of  the  medium:  (a),  (b),  reconstructed  and  D  distributions,  respectively,  when  NIR  data  only  are  used;  (c)  (d),  reconstructed 
and  D  distributions,  respectively,  with  ultrasound  localization. 


calization,  improvements  in  terms  of  spatial  distribu¬ 
tions  were  higher  than  those  when  the  target  was 
located  at  the  boundary. 

B.  Reconstruction  when  the  Direct-Current  System  is 
Used 

Figures  9(a)  and  9(b)  show  the  simultaneously  recon¬ 
structed  |jl,q  and  D  distributions  of  the  high-contrast 
target  located  at  the  center.  The  reconstructed 
mean  was  only  15.6%  of  the  calibrated  value,  and 
the  reconstructed  mean  D  was  —132%  of  the  cali¬ 
brated  value.  Figures  9(c)  and  9(d)  are  the  simulta¬ 
neously  reconstructed  absorption  and  diffusion  maps 
with  ultrasound  localization.  Figure  9(e)  shows  the 
C-scan  ultrasound  image  obtained  at  the  source- 
detector  plane  and  used  for  the  localization.  The 
reconstructed  mean  has  been  improved  to  22%  of 
the  calibrated  value  and  the  reconstructed  mean  D 
has  been  improved  to  121%  of  the  calibrated  value. 
The  FWHM’s  of  and  the  D  maps  have  been  im¬ 
proved  by  50%  and  30%,  respectively. 

Figure  10  shows  the  experimental  results  of  the 
low-contrast  target  located  at  the  center  of  the  me- 
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dium.  Figure  10(a)  is  the  jx^^  distribution  recon¬ 
structed  by  using  NIR  data  only.  The  spatial 
distribution  of  ix^  is  much  worse  than  the  high- 
contrast  case.  The  reconstructed  mean  absorption 
coefficient  is  only  47%  of  the  calibrated  value,  and  the 
FWHM  is  around  1.6  cm.  Figure  10(b)  shows  the  D 
distribution  simultaneously  reconstructed  with  jx^  by 
using  NIR  data  only.  The  reconstructed  mean  value 
was  —140%  of  the  calibrated  value,  and  the  FWHM 
was  around  1.4  cm.  Figures  8(c)  and  8(d)  are  the 
simultaneously  reconstructed  |x^  and  D  distributions 
with  coregistered  ultrasound  localization.  With  ul¬ 
trasound  localization  the  accuracy  of  reconstructed 
mean  fx^  was  improved  to  73%  of  the  calibrated  value, 
and  the  FWHM  was  improved  to  1  cm,  while  the 
accuracy  of  reconstructed  mean  D  was  126%  of  the 
calibrated  value  and  the  FWHM  was  improved  by 
40%. 

In  Table  2  are  summaries  of  the  experimental  re¬ 
sults  of  high-  and  low-contrast  cases  when  the  target 
is  located  at  different  positions.  Generally  the  re¬ 
constructed  absorption  coefficients  are  much  lower 
than  those  obtained  from  the  frequency-domain  sys- 


BEST  AVASLAB.LE  COPY 


Table  2.  Experimental  Results  from  the  dc  System 


NIR  Only 

NIR  +  US 

Calibrated  Value  from 
Dartmouth  College 

Target  Location 

Target 

Contrast 

Mean  Mua 

(cm“b 

Mean 

D  (cm) 

Mean  Mua 

(cm"b 

Mean 

D  (cm) 

Mean 

Mua 

(cm"h 

Mean  D 
(cm) 

Close  to  the  boundary 

At  the  center 

High  contrast 
Low  contrast 
High  contrast 
Low  contrast 

0.0485  (19.4%°) 
0.047  (67%) 
0.039  (15.6%°) 
0.0328  (47%) 

0.043  (123%) 
0.044  (126%) 
0.046  (132%) 
0.0486  (140%) 

0.059  (23.6%°) 
0.049  (70%) 
0.055  (22%°) 
0.051  (73%) 

0.043  (123%) 
0.046  (132%) 
0.042  (121%) 
0.044  (126%) 

0.2-0.3 
0.07  ±  5% 
0.2-0.3 
0.07  ±  5% 

0.0347  ±  5% 
0.0347  ±  5% 
0.0347  ±  5% 
0.0347  ±  5% 

"The  percentage  is  calculated  on  the  basis  of  the  mean  value  of  0.2-0. 3  cm  ^ 
US,  ultrasound 


tern,  and  the  reconstructed  diffusion  coefficients  are 
much  higher  than  those  from  the  frequency  domain. 
The  improvements  with  ultrasound  localization  are 
much  less  than  that  obtained  from  the  frequency- 
domain  system. 

Comparing  the  results  from  frequency-domain  and 
dc  systems,  we  find  that  the  reconstructed  results 
from  the  frequency-domain  system  are  much  better 
in  terms  of  accuracy.  With  ultrasound  localization 
the  mean  absorption  coefficients  obtained  from  the 
frequency-domain  data  are  approximately  111-113% 
of  the  calibrated  value  for  the  low-contrast  target  and 
around  64-65%  of  the  calibrated  mean  value  for  the 
high-contrast  target,  while  they  are  only  70-73%  and 
22-23.6%  for  the  same  cases  in  dc.  The  mean  diffu¬ 
sion  coefficients  obtained  from  the  frequency-domain 
data  are  within  102-108%  of  the  calibrated  value, 
while  they  are  approximately  121-132%  for  the  same 
cases  in  the  dc.  The  FWHMs  obtained  from  the 
frequency-domain  reconstructions  are  always  within 
90-130%  of  the  calibrated  value,  while  the  FWHMs 
are  approximately  110-160%  for  the  same  cases  in 
the  dc. 

6.  Discussion 

The  reconstructed  absorption  and  scattering  distri¬ 
butions  can  be  partially  improved  by  iteratively  up¬ 
dating  the  Jacobian  matrix  to  account  for  higher- 
order  terms  of  the  weight  matrix.  However,  this 
procedure  is  time-consuming  and  prohibits  near-real- 
time  NIR  image  processing.  The  iterative  updating 
of  a  localized  Jacobian  matrix  within  the  ultrasound 
specified  region  could  further  improve  the  accuracy  of 
reconstructed  absorption  and  diffusion  coefficients 
with  reduced  computation  load.  Currently  we  are 
pursuing  this  research,  and  results  will  be  reported  in 
the  near  future. 

In  the  reported  phantom  studies  we  assume  that 
the  lesions  are  isolated  and  embedded  in  a  homoge¬ 
neous  background.  Therefore  we  localize  the  image 
reconstruction  by  using  the  target  geometry  obtained 
from  ultrasound.  In  the  clinical  breast-imaging 
cases  the  background  tissues  also  scatter  and  absorb 
the  diffusive  light,  which  are  also  accounted  for  in  the 
dual-mesh  reconstruction  scheme  reported  in  Ref.  20. 


7.  Summary 

We  have  reported  experimental  results  from  the 
simultaneous  reconstruction  of  absorption  and  scat¬ 
tering  maps  with  ultrasound  localization.  Trans¬ 
mission  geometry  was  used  for  NIR-imaging 
reconstruction,  and  standard  pulse-echo  ultrasound 
was  used  for  obtaining  coregistered  ultrasound  im¬ 
ages.  When  a  frequency-domain  system  was  used 
the  ultrasound  localization  results  showed  that  the 
reconstructed  mean  absorption  coefficients  of  high- 
contrast  phantom  targets  improved  more  than  20%, 
while  improvements  to  the  reconstructed  mean  dif¬ 
fusion  coefficients  were  moderate.  For  low-contrast 
targets  no  significant  improvement  has  been  shown 
in  the  mean  reconstructed  values.  However,  im¬ 
provements  in  the  target  spatial  distributions  are 
significant.  When  a  dc  system  was  used,  improve¬ 
ments  in  using  ultrasound  localization  were  small  in 
terms  of  mean  reconstructed  values.  Owing  to  the 
lack  of  phase  information  and  possible  cross  talk  be¬ 
tween  absorption  and  scattering  coefficients,  the  dc 
system  performance  is  much  worse  than  that  of  the 
frequency-domain  system.  With  ultrasound  local¬ 
ization  the  reconstruction  speed  has  improved  by  a 
factor  of  10  and  near-real-time  optical  imaging  be¬ 
comes  feasible. 
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development  of  the  dc  system  and  the  mechanical 
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toms.  The  authors  thank  the  following  for  funding 
support:  Donaghue  Foundation,  U.S.  Department 
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00-1-0217,  DAMD17-01-1-0216),  and  the  National  In¬ 
stitutes  of  Health  (1R01CA94044-01A1).  The 
authors  greatly  appreciate  the  donation  of  a  proto¬ 
type  dc  system  by  Multi-Dimensional  Technology  Inc. 
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